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Abstract 



The enormous impact of computers on everyday life is an incredible fact. The current 
technology of the computers is based on th technology of transistors, which is an appre- 
ciated synergy between computer science and quantum physics. Smaller transistors can 
work with less power, can be highly densed and can be switched on and off very fast. 
While the miniaturization is necessary to increase the computation power, this minia- 
turization also provides us with an intuitive way of understanding why, in a near future, 
quantum laws will become important for computation. 

Gordon Moore (0) in 1965 codified his law, known as Moore's law, on the growth of 
computers. It states that the power of computers will double for constant cost roughly 
once every two years. This statement amazingly has come true for decades, since the 
1960s. This growth is achieved through miniaturizing the size of the computers, by 
compacting the number of transistors. There is of course a limit for the number of 
transistors per chip and at present time it is almost 10 18 for a typical size of circuit 
components to be of the order of 100 nanometers. Then this is a fact that around 2020 
we will reach the atomic size for storing a single bit of information, hence quantum effects 
will become incredibly dominant. 

So far, we have considered the increase of the computational power of the current 
computers, so called classical computers. These nowadays commonly being used com- 
puters go back into not a well-known origin in their very old history. Up to now, the 
technology of computers has been improved incredibly. However, Alan Turing, a great 
mathematician, pioneered the concept of what is called now programmable computer (0) 
and this concept yields to the strong Church- Turing thesis (J3|)(l4|). 

"Any model of computation can be simulated on a probabilistic Turing machine with 
at most a polynomial increase in the number of elementary operations required. " 

Even although the polynomial difference in speed might be still significant, this state- 
ment verifies the equivalency of all the computers, from the most pioneer ones, to the 
current supercomputers. Whereas, intuitively, if a quantum computer works truly based 
on quantum postulates namely superposition of states, then an exponential enhancement 
would be achieved by massive operations on a total quantum state. 

The situation on our desire for a quantum computer should be clear in the sense 
that it's not just our intuition for a far future but it is somehow our current demand, 
too. Particularly, for some special cases, the atomic scales have already been achieved 
(0). Then classical physics cannot explain the evolution of the corresponding system any 
more. 

One considerable step through our desire for a quantum computer has been established 
by Feynman, 1980's (0). He suggested that a quantum computer based on quantum logic 
would be ideal for simulating quantum-mechanical system by making a hint on how to 
get involved to these problems. 

"So, as we go down and fiddle around with the atoms down there, we are working 
with different laws, and we can expect to do different things. We can manufacture in 
different ways. We can use, not just circuits, but some system involving the quantized 
energy levels, or the interactions of quantized spins, etc. " 
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Apart from the serious requirement on quantum algorithms, one of the most important 
points that has to be considered for any quantum physical system is that the system should 
be stable, while the operations and processing are going on. This statement explains the 
requirement of the maintaining coherent superposition of quantum states. 

Generally speaking, decoherence is one of the ultimate obstacle to the practical re- 
alization of a quantum computer. For any quantum computer, which has inevitable 
interaction with environment, the quantum information is decaying and this means that 
the quantum system is suffering from decoherences. Then, errors are provided by de- 
coherence. Hence, there are also another sources of errors, such as the imperfections in 
the quantum computer hardware. By the way, it is inevitable to have an appropriate 
quantum error correction code, which can eliminate the errors and make the quantum 
processing stable. 

For a while it has been assumed that quantum error correction codes would not be 
accessible due to some quantum fundamental difficulties in interchanging the classical 
error correction code to the quantum counterpart. However, quantum error correction 
code has been invented by Peter Shor (0) and Andrew Stean (p), provided that the 
probability of error per computational step is not large. This field has been extended 
much after these pioneer works and we will revisit some of the following works in this 
thesis. Therefore, employing quantum error correction codes, the problem of the random 
error due to decoherence is assumed to be solved out. Nevertheless, the threshold on the 
error probability is still so strict that it is not trivially accessible. 

In this work, the objects of quantum algorithms and quantum error corrections will 
be reviewed in separate sections. 

Chapter 1 Classical computers, as much as being advanced, still can give only polyno- 
mial enhancement compared to the previous not advanced versions. A quantum 
computer can be much more faster, exponentially, as it involves superposition of 
states. Then, any process is applied to the whole massive quantum states. The 
result of the processing is however hidden in between the superposed states. Then 
a large number of outputs would be only potentially possible and totally useless. 
Nevertheless, it is the task of quantum algorithms to exploit the inherent quantum 
parallelism of quantum mechanics and to highlight the desired output. Therefore, 
quantum computers, in order to be useful, require the development of appropri- 
ate quantum software, that is of efficient quantum algorithms. As we will see in 
the following sections, entanglement is introduced as an essential prerequisite for 
some of these quantum algorithms specially those with exponential enhancement 
over classical processing. Furthermore, quantum information is decaying, which 
means that quantum system is suffering from decoherence. Errors are provided by 
decoherence. Hence, there are also another sources of errors, such as the imper- 
fections in the quantum computer hardwares. By the way, it is inevitable to have 
an appropriate quantum error correction code which can eliminate errors and make 
the quantum processing stable. As we will see, in the following sections, quantum 
teleportation or entanglement plays a crucial role in a fault tolerant quantum com- 
puter. Therefore, in order to realize a quantum computer in its real meaning as to 
be quantum digital computer and to be different from the classical analog comput- 
ers, furthermore in order to get the advantages of quantum computer by running the 
quantum algorithms, entanglement should be provided among the quantum states. 
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Chapter 2 Therefore, we conclude that entanglement might be a resource for a quan- 
tum computations's power through the theoretical point of view. Entanglement is 
important for a number of quantum algorithms, also it is an essential requirement 
for fault tolerant quantum computing. Then, one challenge would be understanding 
the quantum entanglement in more details. Enormous studies on the concept of 
quantum entanglement exist. We study the problem in its corresponding chapter. 
We will emphasize parts, which will be used for our introduced new measure of 
entanglement, applicable for the particular experimental system. 

Chapter 3 After all, the main goal of this study would be reasonably focused on gen- 
erating or even manipulating quantum entanglement. Nuclear magnetic resonance 
(NMR) based quantum information processing seems to be a good candidate for 
this idea as there are variety of seemingly successful implementations of quantum 
algorithms with NMR. We first explain NMR system and NMR quantum infor- 
mation processing and quantum computation, as much as we need for this work. 
Then, we will prove, by studying a particular quantum algorithm, that NMR at 
high spin temperature does not give a true quantum processing . The experimental 
conditions in order to fulfill all the requirements for generalizing an entangled state 
will be shown to be very strict, almost out of reach. However, there are still some 
indirect approach to solve this problem. This will be discussed extremely in its 
corresponding part. 

Chapter 4 In order to generalize a true quantum entanglement, while getting benefits 
of the inevitable advantages of the previously studied NMR systems, we move to 
a rather different physical system, so called Electron Nuclear DOuble Resonance 
(ENDOR). We study ENDOR system with several experiments and finally we go 
to our final goal of realization of a true quantum entanglement. This will be shown 
by several experiments on introducing ENDOR system, which is a rather more 
complicated experimental scheme as compared to NMR because that ENDOR is 
a double resonance experiment and needs to have simultaneously controls of radio 
frequency and microwave frequency. After all, we will show ENDOR for quantum 
computing experiments under different conditions, magnetic fields from 9.5 GHz up 
to 95 GHz and temperature from ~ 300 K down to ~ 3 K. Results demonstrate 
truly establishment of (pseudo-)entangled state. This will be discussed to be par- 
ticularly important because it demonstrates our ability to manipulate electron and 
nuclear spins on our arbitrary desired values and conditions. Also, we will introduce 
an almighty molecular sample, which will be discussed to be a proper sample for 
quantum computing even for a larger number of qubits or under very different ex- 
perimental conditions of magnetic field and temperature. Finally, we will go much 
closer to the realization of quantum entanglement as it is based on experimental 
conditions just to be in hand. 
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Preface 



Once in 1850, William Gladstone, the British minister of finance asked Faraday, "If 
electricity had any practical value". He replied, "One day, sir, you may tax it". 

Soon or later, this is a dream though, we may be asked to pay monthly also, for purified 
entanglement supplied by Entanglement Purification Centre of our city. Entanglement 
would be purified there and sent to our houses via teleportation similar way that can 
exchange the used entangled states with refresh and pure entangled ones. We may use it 
then for running the quantum part of our computers for computing also for information 
processing. Let's go further in such a sweet dream, those days people may know all about 
the dynamics of entanglement, and how entanglement plays role for macroscopic worlds 
of physical properties, down to its explanation for gauge field theory. I am but still not 
so brave to say that it will be a key for unification in physics but I'd like to say it will 
work for small systems of course and even large systems as then we know how to use it 
in place for Partition functions of large number of molecules. Such a world would be a 
wonderful place and I wish I could choose my life time to enjoy the smell of knowledge 
in such an era. 
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Chapter 1 



Classical Computing/Quantum 
Computing 

Let's start with this question. What do we want from a quantum computer? In principle, 
all the calculations that can be performed in a classical computer can also be performed in 
a quantum computer, simply by replacing the irreversible gates of the classical computer 
with their reversible counterparts. The new circuit, which includes only the reversible 
gates, representing the unitary operations, can be implemented in a quantum computer. 
Then, the question remains that what is the advantage if we construct such a quantum 
computer only by interchanging from a classical computation to a quantum computer. 

The appeal of a quantum computer is the ability of quantum algorithms faster than 
classical ones. This is the core statement for any attempt regarding improvement in 
quantum computation. 

The most welcome advantage of a quantum computing compared to the classical one 
is for the case of having a polynomial complexity in solving a particular problem with a 
quantum computer versus an exponential complexity of solving the same problem with 
a classical computer. Shor's algorithm of integer factorization is a specific case of this 
class. 

In this chapter, we review the concept of quantum computing very briefly in order to 
find out the inherent advantages to the classical computation. In this regard, we review 
quantum algorithms to verify the advantages of each case compared to the classical ones. 
Quantum algorithms must use quantum features, which are not available in classical 
computers, in order to improve the computational complexities. The quantum algorithms, 
which have been defined up to now, can be separated to two main classes of the database 
search algorithms and the algorithms for finding the generators of a normal subgroup of 
a given group. We will see that quantum superposition is the critical quantum property 
for the former class including Grover's algorithm, which gives quadratic enhancement for 
an unsorted search problem. While for the latter class, which includes Shor's factoring 
algorithm, entanglement plays a crucial role. 

However, most of the hard problems are NPC, NP complete, and it is most likely 
that a quantum computer also cannot solve them in polynomial time. Hence, building a 
very expensive quantum computer, with which only very rare quantum algorithms give 
exponential enhancement, seems not to be so wise. As a matter of fact, polynomial com- 
plexity class versus exponential complexity class is mostly a mathematical convenience. 
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Sometimes, a polynomial case is even worse than exponential type of complexity and it 
is of course the nature of the problem itself which identifies a good or a bad class of 
complexity. 

In order to look for more elaborated advantages of a quantum computer, we might 
look for more subtle realistic distinctions specially inside polynomial type of complexities. 
This concept is discussed in the part of measurement based quantum computer. We will 
see that entanglement plays the conceptual role, in this regard too. 

Then the straight conclusion would be that quantum features are reasonably cru- 
cial for quantum computational advantages over classical computation. While quantum 
superposition, coherence, is an important feature, entanglement also might be very im- 
portant in this sense. 

In this chapter our aiming would be on giving enough motivations for the study of 
quantum entanglement through emphasizing its critical role in quantum algorithms, par- 
allelism and even more importantly its unique role for a complete fault tolerant quantum 
computing. Quantum algorithms, measurement based quantum computing and fault 
tolerant quantum computing are explained briefly in this regard only to emphasize the 
critical role of entanglement rather than to be the concrete explanation of each topic. 



1.1 Classical Computing 

In a classical computer data are stored in a sequence of bits that are represented with 
a string of O's and l's. Computational processing is applied on this sequence in discrete 
steps such as Boolean gates, (0). 

A classical computational task is not a single task such as finding the primarity of 
a particular number but it is a family of similar tasks, such as "given an integer N, is 
N prime?" Input size is defined to be the space in computer that is needed to store the 
particular input N. 

Now, we want to know how hard is the particular task to be done on the computer. 
The computational effort or complexity is then given by finding out that how the compu- 
tational effort grows as the input size grows. 



1.1.1 Polynomial Versus Exponential Complexity 

Growth in the number of computational steps as the input size grows is mainly divided 
to two classes of polynomial and exponential complexities. The former case is known as 
problems, which are feasible computation on a classical computer while the latter one 
is known as unfeasible computation. Very wellknown example is the problem of finding 
the prime factor for a given number. Up to now, there is no classical algorithm, which 
can solve this problem efficiently, meaning that the number of steps that are required 
to solve the problem is not bounded by any polynomial function. Then, the problem of 
factoring is known as a classical hard problem, whereas quantum Shor's algorithm gives 
an efficient way to solve this problem. 
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1.1.2 Query Complexity 



The complexity of a classical task would also be explained in a query complexity. The 
assumption is that a "black-box" is given, which computes a function / : x — > y. During 
the processing, each call of the query is counted as one step. 

For example, for / : n bits \— > 1 bit, the problem is that whether / is balance or 
constant. The promise is that / is either balance or constant. For a classical computer 
the query complexity of this particular example is 2 n_1 + l, in a worst case. However, with 
a quantum computer, this problem can be solved with a single query to the black-box. 
The point is that the property, which has been asked is a global property. Therefore, 
there is no need to examine all the elements and compare them to each other, in order to 
find the answer. Quantum superposition makes it possible to manipulate all the elements 
of initial states, simultaneously and extract the desired global property with only one 
query to the black-box. 



1.2 Quantum Computing 

Generalization of the above concepts to the quantum computing, gives the corresponding 
quantum complexities. In quantum computation, bits are represented by qubits, two- 
level system with chosen computational bases. The computational steps are quantum 
operations on qubits and are represented by unitary gates or quantum measurement. 

There are various models possible for structure of a quantum computer. Each case 
is evaluated by a particular computational complexity, and entanglement is particularly 
an important quantum property for each class to represent the desired advantages of the 
quantum complexities compared to the classical counterparts, (jjj. 



1.2.1 Gate Array Model 

This is the standard model for quantum computation. It starts with row of qubits in a 
simple state, generally initialized state. The computational steps are unitary gates and 
measurement that is performed at the end of the computation in {|0), |1)} basis. 

Quantum algorithms prescribe the quantum processing and quantum computational 
steps. Quantum algorithms use the properties of quantum physics in order to provide new 
modes of computation and information processing, which are not available to classical 
computers. Because of the fundamental postulates of quantum physics, namely superpo- 
sition, a quantum computer, which works based on the laws of quantum physics is able to 
store and process large volumes of information. All the information is represented in a(n) 
(entangled) quantum state. However, generally the access to the whole information is 
strictly restricted by quantum measurement theory. Then, only a relatively small amount 
of the whole stored information can be read out. 

However, the extracted information maybe of a global nature and then gives an intu- 
itive information on the whole quantum physical system. On the other hand, this process 
may be impossible for a classical computer to be performed efficiently to give a global 
feature of the total system. Because a classical computer may still need to examine each 
individual state in order to get the total information. 
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Quantum algorithms are actually to give a routine on a way how to extract the 
desired information from the total information stored in superposed quantum states, in 
a rather more efficient way compared to the classical approaches. Particularly, the most 
celebrated quantum algorithm is Shor's algorithm for integer factorization. While the 
best known classical algorithm for this problem runs in superpolynomial time of order 
exp(n 1//3 (logn) 2//3 ), Shor's algorithm provides a method for factoring an integer of n digits 
in the order of 0(n 3 ). 

The significance of all quantum algorithms lies in the issue of computational com- 
plexity theory. Then, we classify our review based on the complexity classes for quantum 
algorithms. 

The computational step for the gate array model of quantum computation is counted 
as the required unitary gates with dimension equal to or less than two qubits, e.g. CNOT, 
Hadamard and so on, in addition to the read-out step, which is the measurement at the 
end of the quantum circuit and only is performed in the computational basis. 



Polynomial Versus Exponential Complexity; Shor's Algorithm 

Suppose TCn to be the Hilbert space with dimension N and with basis |0), |1), \N — 1), 
being labeled by integers module N. For this space, quantum Fourier transformation 
QFT^ is defined as unitary transformation 7i N 7i N as follows, (0), ([left. (11), fl2h . 



■ ■ M . ■ 



QFT w |x) = 7 =^e^|y), (1.1) 

ViV j=Q 

for x, y = 0, 1, N — 1. The matrix elements are represented as 



[QFT) xy = A=W*», (1.2) 



where 



2iti 



W = e^r. (1.3) 



For example, QFT 2 is just the Hadamard gate. 

In a classical scheme with dimension N = 2 n , Fourier transformation is an iV x iV 
matrix. So, calculation of the matrix needs 0(N 2 ) steps. The number of steps can be 
decreased to 0(N log N), i.e. still exponential function of n. 

However, it is proved that QFT W can be implemented in 0((logiV) 2 ) = 0(n 2 ) com- 
putational steps. If iV ^ 2 n then QFTjy cannot be exactly implemented in poly(log) 
steps, but QFTjy can be still be closely approximated (to any accuracy e) by QFT^, 
where M > N, M = 2 m and M = 0{N log N). i.e. we need only O (log log N) additional 
qubits. 

This is a remarkable fact that QFTjy can be implemented polynomially and this fact 
is fundamental step in Shor's algorithm. 

The problem of factoring can be reduced to the problem of periodicity finding. This 
connection is entirely mathematical and does not have any quantum physical aspects. 
Therefore, if we could have an efficient classical algorithm for solving the periodicity 
problem then factoring problem also could be solved entirely classically. 
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Suppose that the integer number N is given. Pick an integer N 2 < q < 2N 2 . This step 
is clearly always possible. Also, pick a random number integer x, such that gcd(a;, N) — 1. 
This step can be run using Euclid's algorithm for calculating the greatest common divisor, 
b = gcd(x, N). If b > 1 then b is a factor of iV and actually we are done. Otherwise keep 
the choice of x. 

Now, create two quantum registers. Register I must contain enough qubits to rep- 
resent numbers as large as q — 1. Register II must contain enough qubits to represent 
numbers as large as N — 1. Then, load quantum register I with an equally weighted 
superposition of all the integers from to q — 1. Load register II with zeros. The total 
state at this point will be 

4 £ !«.<»■ 

Now, apply the transformation x a modN to each number stored in register I and store 
it in register II. The total state is changed into 

1 9-1 

— ^2\a,x a modN). (1.5) 
a =o 

Measure the state of the second register. The measurement must yield one value of 
x a modiV, let's say V. The measurement on the register II has the effect of collapsing 
the register I into an equal superposition of each state between and q — 1, such that 
a^modiV = V, because these are the only states, which could give that results. 

By Euler's theorem, there is the lowest power r such that 

x r = lmodiV, (1.6) 

where r is the periodicity of f(x) = x a modN. 

Up to now, QFT^y is explained to be implementable efficiently Also, factoring problem 
is shown to be a periodicity problem. Then, the remaining step in Shor's algorithm is to 
compute the QFT^ on the collapse register I. 

If the periodicity of QFT^ is r, then the QFT^ output has the periodicity q/r. With 
a classical computer the r period is evaluated. Now, the factor of N can be determined 
by taking the greatest common devisor, gcd(x r//2 + 1, N) and gcd(x r//2 — l,N). If this 
does not work, the same process should be tried for a new number x. 

The important point supporting the efficiency of Shor's factoring algorithm is the fact 
that QFT can be implemented efficiently, then exponentially faster than the classical 
case. The following example makes it clear how Shor's algorithm works. 

Problem is to find the prime factors for iV = 15. Firstly, pick an integer q such that 
N 2 <q< 2N 2 . 

(N 2 = (15) 2 = 225) < (q = 2 8 = 256) < (2iV 2 = 2(15) 2 = 450). (1.7) 

This implies that 8 qubits are sufficient for register I. The second register must contain 
enough qubits to represent numbers up to 15 — 1 = 14. Let register II holds numbers up 
to 2 4 = 16. 

Therefore, we are working with two registers, I, and, II. 
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Register I : given as 



100000000) + 100000001) + 100000010) + ... (1.8) 

or equivalent ly 
|0) + |1) + |2) + ... 

Register II : all zeros 

|0000) = |0). (1.9) 
The total state would be then as follows 

1 [|0) + |1) + ...]|0>. (1.10) 



v/256 

Now, apply the transformation x a modN to each number stored in register I and 
store it in register II. The number x is picked at random here to be x = 7, satisfying 
gcd(x, N) = 1. 

An important result from number theory is shown as follows. 
f(a) = x a modN is a periodic function, (x and a and iV are all integers; the greatest 
common divisor of x and N is 1.) 

Definition: BmodC is the remainder after division of B by C. 

This is known that f(a) = x a mod N has a period P. If P is an even number then 
(x p l 2 — l)(x p l 2 + 1) is an integer multiple of N. 

Hence, we have 

First register 7 a mod 15 Second register (1-11) 

|0) 7° mod 15 1 

|1) 7 1 mod 15 7 

|2) 7 2 modl5 4 

|3) 7 3 modl5 13 

|4) 7 4 modl5 1 

|5) 7 5 modl5 7 

|6) 7 6 modl5 4 

1 7) 7 7 modl5 13 

1 8) 7 s mod 15 1 

Now, measure the second register. For this example, the possible outcomes are |1), |4), |7), |13). 
E.g. if the measurement on the second register gives V = |1), then the first register must 
collapse into |0) + |4) + |8) + |12) + |16) + |20) which has the periodicity of 4. Then 
apply the quantum Fourier transformation (QFT). 

We know that the discrete classical Fourier transformation is as follows. 
Input: a series of iV complex numbers, x ,x 1 ,x 2 , ...,Xn- 
Output: a second series of iV complex numbers, yo,yi,y2, ■■■jVn, where 

JV-l 
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#of outcomes 



64 128 192 



Figure 1.1: The result of the measurement. 



For the input series, if it has a period P, the discrete Fourier transform (DFT) output 
will repeat with a period N/P. 



Input 
10000000 
10001000 
10101010 
11111111 



Period(p) DFT 
8 

4 => 
2 => 
1 => 



Output 
11111111 
10101010 
10001000 
10000000 



Period(N/P) 
1 
2 
4 



;i.i3) 



The quantum Fourier transform (QFT) is almost the same as the discrete Fourier trans- 
form, however it operates on the complex numbers which represent phase and amplitude 
of the various contributions of a superposition of states. For example consider three-qubit 
case. 



= aooo|000) + aooi|001) + aoio|010) + a ii|011)+ (1.14) 

|100) + aioi|101) + QfnolllO) + a m |lll) 
= aooo|0) + «ooi|l) + «oio 1 2) + «on|3) + «ioo|4) + aioi|5) + a uo |6) + a m |7). 

For this general case, suppose that the amplitudes are equal to each other. Then the 
Fourier transform would be written as follows. 

Input Period(p) QFT Output Period(N/P) 

|0> 8 =► {|0) + |1) + |2) + |3)+ 1 

|4) + |5) + |6) + |7)} 
|0) + |4) 4 => |1) + | 3 ) + |5) + |7) 2 

|0) + |2) + |4) + |6) 2 => |0) + |4) 4 

{|0) + |1) + |2) + |3)+ 1 => |0> 8 

|4) + |5) + |6> + |7» 

(1.15) 

As with the discrete Fourier transform, if the period of the N input is P, then the 
period for the output of the quantum Fourier transform is N/P. Then in the above 
explained example for N = 15, we expect the periodicity of 256/4 in the output. The 
new state can be written as |0) + |64) + 1 128) + 1 192) for q/r = 256/4 = 64. With a 
high probability, the measurement yields an integer which is an integer of g/4. If the 
measurement were to be repeated many times, the output then with high probability 
could be plotted as Figure fTTTl 

Now, it turns out that gcd(7 2 + 1, 15) and gcd(7 2 — 1, 15) or 5 and 3 are the factors. 
If this does not work, the process would be repeated for another random number x. 
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Query Complexity; Grover Algorithm 

The black-box is a unitary operator, which can be applied as a single gate. Then the 
complexity is counted with the number of query to the box. 

Grover's algorithm is for searching in an unstructured database for a single item, (0), 
ffl6h . flit), flih . The particular problem is among N database, to find the special x 
value for which f(x = xo) = l, where / is a black-box defined as 

/: {1,2,..., A} -{0,1}. (1.16) 

Unstructured search simply means that if the answer for the trial is "NO" then we 
get absolutely no idea where the answer could be. Instead, for an structured search 
the situation is completely different. For example, suppose an ordered search problem 
with a given N = 2™ number in ascending order. Suppose that one of them (xq) is in 
hand. Then try to find Xq in the list. Classically, log 2 N = n queries are necessary and 
sufficient to solve the problem. With quantum computation, however it is different, hence 
representing the quantum benefits. 

We assume an unstructured database search problem. Classically, O(N) queries are 
required to solve this problem. But, in a quantum case, 0(yN) queries are necessary. 
This is given by Grover algorithm and it is proved that "no quantum process can do this 
with fewer steps", fjlih . 

The important point to be discussed here before starting with Grover algorithm is to 
make sure about the benefits of quantum searching. We try to go around this question 
with a very simple example. Assume that we are given a database including N values, for 
which f{x = Xq) = 1 and otherwise f(x) = 0. The task is to find x = Xq. The important 
point is that we do not have different compositions of all the states. Instead, we have 
only an equal superposition of items. Then making an equal superposition of the items 
need no exponential effort. However, if we originally have all the different cases, then 
to make an equal superposition of all the items requires exponential steps and hence the 
benefits of Grover algorithm is lost, (0). 

We start with a black box / : nbits — ► lbit with f(x ^ Xq) = and f(x = Xq) = 1. 
However, it turns out that it is more convenient to use equivalent box, phase box as 
follows 

'J*> = f |x > if f ^ X0 ' d-17) 

I — \X) if X = Xq. 

The equivalency becomes clear if we try ot for ^75 (|0) — |1)). In other words, one use 
of I Xo needs just one query to Uf. The Grover iteration is introduced as follow 

Q = —H n I H n I Xo , (1.18) 

where 

H n = H ® ...®H (1.19) 

with H to be the Hadamard operation. The starting state is an equal superposition of all 
possible x-values as follows 

|^>=#n|0...0> = -L5>>. (1.20) 
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Then (n + l)th state is written based on the previous nth state as follows 

= Wn>- (1.21) 

For any state in space, we have 

I W = I-2|V>)M, (1.22) 

indicating that I\m is a reflection in H±(ip), orthogonal subspace. The proof is clear 
because that for any chosen state |£) = a\ip) + b\£'), where |£')_L|^>), then = 
-a\%>)-b£). 

I\m preserves the 2D subspace spanned by \ip) and Also, for any unitary operation 
U, we have 

UI w U^ = I um , (1.23) 

because 

U(I -2\ip)(ip\)rf = I -2\Uip)(Uip\ (1.24) 
= Iu\i>)- 

If all the above properties are applied on the iteration operator Q and also recall that 
H n = Hi, then 

Q = -H n I H n I X0 (1.25) 
= -Ih„\o)Ix 

In the same manner, both and I XQ preserve the 2D plane, spanned by |^o) and \xq). 

Now, we want to show, that —I\^ ) is actually inversion in the average. Recall j^o) = 
-j= Y2i\x). For | a) = YlH=o a x\x), with average amplitude a = ^j^ 2 -, then 



-I m \a) = -(I-2|Vo)(V'o|)|«) (1.26) 

2 

N 



= -yXax - 2a)\x) 

X 

= -y]a x - 2(a x -a)\x) 

X 

X 

where a' x = a x — 2(a x — a) and (a x — a) is the distance of a x from average, see Figure 

The total steps of the Grover algorithm can be represented as follows 
. Start with a uniform superposition, Figure 11.31 
. Apply I xo , hence the inverse of x , Figure IP1 
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average a 



Figure 1.2: a! x is the inverse of a x as to the average a. 



1/V2 



average 



Figure 1.3: The initial states are in a uniform superposition. 



1/V2 



new average 



Figure 1.4: Applying I XQ gives the inverse of x and makes the average less than the 
previous one. 
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previous 
average 



Figure 1.5: — 1\^ ) gives the inverse in average. 

. Apply — I\ip Q ), which gives the inverse in average, Figure [T3J 

Amplitude for x = xq much increased, while amplitude for i ^ xo has been slightly 
decreased. Now, repeat the iteration operation. 

It is proved that in 0(^/N) the value Xq can be extracted and no other quantum algo- 
rithm can solve this problem faster. Clearly, the scheme is based on quantum properties 
such as superposition of states. 

Discussion 

Our review on quantum algorithms up to now has been mainly on the quantum compu- 
tational algorithms with the main idea of emphasizing the quantum benefits compared 
to the classical counterparts. Shor's algorithm is introduced as a quantum algorithm, 
which represents polynomial complexity compared to any classical known algorithms, 
generally with exponential complexities. Entanglement is an essential quantum property, 
which supports exponential speed-up for Shor's algorithm. Grover algorithm represents 
quadratic query complexity compared to classical algorithms. Original algorithm rep- 
resented by Grover, involved entanglement in order to represent quantum enhancement, 
however quantum coherence or superposition seems to be the required quantum property, 
(0- 

There are some other quantum algorithms, but we are not going to get involved for 
their explanations here. Through our study on quantum entanglement, we emphasize 
that quantum enhancement is available if entanglement exists between the states even 
though the vice explanation is not correct, generally. Both quantum superposition and 
quantum entanglement play crucial role for quantum algorithms. Then, serious study on 
the concept of entanglement still seems to be reasonable. 

We should maybe go around explanation on quantum information processing, too. 
Namely quantum teleportation and superdense coding. But, here our main idea of study- 
ing quantum algorithms has been finding out the complexity advantages due to quantum 
properties such as superposition and entanglement, compared to classical counterpart 
algorithms for solving the similar problems in quantum gate array model. Quantum tele- 
portation is an essential part of the measurement based model. Then, we postpone it 
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to the next section, when speaking about measurement based model and corresponding 
complexities. Superdense coding turns out to be very important for introducing a new 
measure of entanglement, as we proposed for ensemble quantum computer. Then we will 
go around superdense coding in the following chapter when we discuss on entanglement 
properties and different measures of entanglement. Superdense coding will be explained 
precisely in order to introduce a new entanglement measure. 



1.2.2 Measurement Based Model 

Measurement based models are particularly important for several reasons. Firstly, they 
do not have any evident classical analogous. This means they can be regarded to be 
more intrinsically quantum. Secondly and more importantly for us, measurement based 
models offer new perspective on the role of entanglement in computation. Finally, they 
might be shown to be more interesting for experimental implementations, then to be in 
our interest. Therefore, we study this model of quantum computation. 

In gate array model, measurement is the last part of quantum circuit for converting 
quantum information into classical information in order to read out classical answers. In 
contrast to unitary gates, measurements are very destructive. In other words, we carry 
large amount of information, all during the calculation in quantum circuit whereas in 
the last step we spoil most part of the whole information. It can be then surprisingly 
beneficial to perform measurement at the most beginning part of calculation or even use 
several measurements for steps of calculation. 



In this formalism of measurement based quantum computation, (|9|), (|20[) . (|21l). (J2 



()23f ). (|24f). (|25l). (|26T). fixed entangled state of many qubits are used to start up the 
computation. Therefore, entanglement plays a crucial role. Two very general schemes 
exist, teleportation quantum computing TQC and Cluster model or one-way quantum 
computer (1WQC). The latter one is based on multi-partite entanglement among all 
the cluster state, while the former one is based on bipartite entanglement. We will 
discuss TQC and give up 1WQC for future because that these two models are shown to 
be equivalent. Though the role of entanglement is so clear in TQC through quantum 
teleportation, as a non-local quantum algorithm. 

After all, we will conjecture that the model can basically give an exponential en- 
hancement, however on polynomial classical algorithms, by parallelism and due to the 
manipulation of entangled states. 



Teleportation Quantum Computation 



Quantum teleportation, (|27l ). is a non-local quantum algorithm for transferring quan- 
tum state between two parts. This algorithm is based on the entangled state preshared 
between the involved parties. 

Suppose that Alice posses a quantum state \ip) and wants to send it to Bob. Each 
of the parties, Alice and Bob, posses a part of an entangled state |/3oo}> where \(3 ZX ) is 
maximally entangled Bell state, defined as follows 

l/SU) - IM> + ™ g> . (1.27) 
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i2<A*#>i|Ax)>23 = ^x d z c ^) 3 - (1-30) 



Totally, Alice has two qubits in her disposal, one to be in the state \ip) and the other 
that is entangled to the qubit in Bob's disposal. Alice applies Bell measurement on 
qubits and gets the experimental outcomes c and d. The state is then teleported to Bob, 
as X d Z c \ip). She sends the classical information, outcomes of the measurement to Bob. 
He applies corresponding unitary operations on the state and corrects it into 

The rotated Bell basis is defined as 

\P(v,) ed ) = tf®I\p ed ), (1.28) 

for any 1-qubit U. Then, consider the situation, where we want to apply the unitary 
transformation U on a state One way is as follows 

Mu)M)M23 =12 (ArfWil/Waa = ^X d Z c |^) 3 . (1.29) 

The last equation is written in a similar manner as the standard teleportation for 
as follows 

1 
2 

The interpretation of the above equations is that given a resource of |/3oo)> we can 
apply any 1-qubit unitary transformation (up to known Pauli operators) to any state 
by measurement alone. Generalization for 2-qubit is done with more elaborated 
teleportation, which we do not mention here but may be found in some review papers, 

a. 

Generalization of teleportation for dimension d of space gives 

12 (0| a) i | </>)23 = da)3> ( L31 ) 

where is teleported into \a)^. Here, the maximally entangled state in d dimension is 
shown by 

1 d-l 

In a same manner, for rotated maximally entangled state, we have 

|0(«)>i 2 = ^®/ 2 |0) 12 , (1.33) 

and hence 

12«>(tl)l<*}l|<«23 = 12«>|lMl«23 (1-34) 

= il/|a>,. 

The unitary operation U is applied on the state by only measurement. 

The unitary transformations (Uj) are chosen such that |0(uj))'s are orthonormal basis 
in 12-space. In the standard teleportation, \<j>) is the Bell basis and U^s are Pauli operators 
I, X, Y, Z and the output state is (Pauli) {0)3. For a general 1-qubit gate Vi then LVs are 

Pauli's and the output state is (Pauli) V\ a) 3. For a general 2-qubit gate V, just apply 
the same but with d = 4 and use 

|0) = 1/2(|00) + |ll))i 2 (|H) + 111))!'*, (1-35) 
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which is the maximally entangled state in four dimensions. 

For a universal computation, we should know, how to perform arbitrary sequence of 
gates from a universal set, by measurement. One set of universal gates is controlled- Z, X- 
rotation and Z-rotation. Hence, we know how to implement the required quantum gates. 
Then, we are almost done if only we can get rid of unpleasant extra Pauli operations 
premium the measurement outcomes. The reason for these extra Pauli operations is the 
randomness of the measurement. 

So, in order to get rid of the extra Pauli operators, we use unrandom measurement 
or adaptive measurements. This means that the result of a measurement determines 
and modifies the basis for the following measurement. Simple computational relationship 
between the Pauli operators and the rotational gates are useful in this regard. 

Suppose we want to do 

Ul-.U^) = X™'Z T l ll ...X™ 2 Z% 2 X™ 1 Z T l 1 (U l ...U 1 )\^}, (1.36) 

where ni and m 8 are the outcomes of ith. measurement. So, we find out the Pauli correc- 
tions for the sequence of the unitary transformations. 

After all, measurement is done in Z-basis for all qubits. The Z Pauli operators do not 
have any effect on Z-measurement outcomes. For X Pauli operators also the situation is 
not complicated. Recall that U\ip} = a\0) +b\l), then XU\ip) = a\l) + 6|0). Therefore, 
probability of outcome Si — 0,1 for Zj measurement on U\ip) is equal to the probability 
of outcome S[ = Si + for Zj measurement on X™*[/|^). The discussion above can be 
easily generalized for controlled Z operation supporting the universality of the gates in 
this model. 

Up to here, we gave a very brief sketch on the TQC. The idea of this studying is to 
show that how this model, which is based on entanglement can give quantum advantages 
over classical ones. The answer is not cleared up yet. But in some cases we can get the 
desired advantages as described in the followings. 

In fact in TQC, we require the "adaptive measurement ". However, if at least for some 
cases the measurement basis choices are not adaptive then we can do all the measurements 
in parallel as they all commute as quantum operations. The special operators, in this 
regard, are known as Clifford operations on n qubits. We recall definitions as follows 

. Pauli group P n on n qubits is generated by ±1, ±i, X, Y, Z on each qubit. 

. C is a Clifford group operation on n qubits if 

CP n C- l = P n i.e. CP n = P n C. (1.37) 

The Clifford group on n qubits is generated by X, Z, H, CX, p = 

Any n-qubit Clifford operator C can be represented by an O(logn) depth circuit. For 
any given general Poly(n) sized gate, there is generally Poly(n) depth so as a special 
parallelisability property of Clifford array are appropriate. 

Therefore, in measurement based model, to apply sequence of Clifford operations on a 
state, it is equivalently correct to apply all measurements, simultaneously and in parallel, 
i.e. Clifford operations are done in depth 1 process. However, Pauli exponents m; and 
rii require further 0{n) classical computation. Adding up the steps, then we conclude 
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Figure 1.6: The depth of calculation is shown for a particular example for a sequential 
adding, where the elements are added step by step to each other. 




Figure 1.7: For the particular example as shown on Figure [OH Adding one by one at 
each level decreases the depth of calculation. 



that any Clifford operation on n qubits can be represented as constant depth quantum 
0(log(n)) depth computation followed by purely classical 0{n) computation. 

We try to re-organize the quantum algorithm into a quantum part and a classical 
part, then minimize the quantum part in its depth. In this sense, entanglement inherent 
in measurement based model has been shown to be useful to come up with quantum 
advantage. 

Suppose the problem for adding up 8 digits. This simple problem may be solved by 
sequential addition as shown in Figure lL6l The sequential addition requires 0(n) depth. 
However, if we add up one by one at each level then the same problem would be solved 
in log(n) depth, which verifies the logarithmic enhancement, see Figure ITT7I 

We note that the quantum gate array model is equivalent to the measurement based 
model. But the equivalency is not generalized to parallelisibility. In that sense, mea- 
surement based model seems more adequate to find out whether quantum algorithms 
exhibit better parallelisibility properties. There are very interesting research outcomes 
for parallelisability for special circuits with special gates, however for the general case of 
problems, still this is an open problem. 

The conjecture is that any Poly(n) sized quantum gate array can be done with 
Poly ( log (n)) layer, and Poly(n) classical computation. This is a remarkable result if can 
be proved. It means that any polynomial time quantum computation can have its quan- 
tum part squeezed into logarithmic depth with supplementary polynomial time classical 
computation along the way. Clearly, the advantage is based on the quantum bipartite en- 
tanglement in TQC and multipartite entanglement on the whole state in 1WQC. This is 
then another motivation for us to study entanglement and to try to realize true entangled 
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state in the experiment. 



1.3 Fault Tolerant; Classical Computing 

Ideally the attempt is that to implement algorithms, classical or quantum algorithms, 
with computers far from any unpleasant attack of noise. However, this is only possible in 
a theoretical study and in reality even in a very clean experiment in a laboratory, this idea 
does not come true. Always there are some noises, through environment or the interaction 
inside the system or what so ever, which make the system unstable. Therefore, in order to 
get benefits of theory, there might be some possibility to overcome the undesired effects 
of the noise. 

Old classical computers encountered problems due to the existence of noise. The 
computational elements, on those old days, had been unstable against noise. Then the 
efforts had been made on making the fault tolerant computations. Nowadays, the classical 
computation is approximately well stable, however the situation for quantum computation 
is still unstable. 

The final aiming of the study of quantum information theory and quantum computa- 
tion is of course the possibility of realizing a true quantum computer, which can realize 
the enhancement of quantum algorithms. Still nobody knows about the most appropriate 
hardware for constructing a quantum computer. Nevertheless, we are quite sure that any 
practical trial for making a quantum computer will be so fragile against the noise. Then 
it is a formidable task to study the noise and the ways to overcome the effects of the 
noise, error correction, and to make a fault tolerant computation. 

In this regards, the experience through the classical computation can be useful, only 
to some extent. In this section, we study the fault tolerant quantum computation after 
a brief overview on the fault tolerant classical computation. We will get it that in case 
of fault tolerant quantum computation, entanglement is a prerequisite property which 
makes it possible to get the fault tolerant computation. Then, we conclude that for a 
reliable quantum computation from the first point of making the quantum algorithms to 
the end idea of correctly realizing the quantum algorithms, entanglement should exist to 
assure the quantum nature of the processing. 

In the current technology of classical digital computers, because of that the existing 
silicon-based elements are remarkably reliable, fault tolerancy is not essential to the oper- 
ation of modern digital computers. However, in the early days of the classical computers, 
when gates were realized using vacuum tubes, the situation was completely different. In- 
deed the gate failure was quite high and a great deal of thought went into how to design 
fault tolerant systems. Von Neumann, (jiih . suggested improving the reliability of a cir- 
cuit with noisy gates by executing each gate many times, and using majority vote. Based 
on his conclusion, if the gate failures are statistically independent, and the probability of 
failure per gate is not so high, then any computation can be performed with a reasonable 
reliability. 

In communication through noisy channel, Shannon, (29), proved that as far as the 
error probability is provided to be less than a certain threshold, so called the capacity of 
the noisy channel, a message can be sent with error probability e using only poly(log(l/e)) 
noisy messages, using suitable codes. The threshold or the capacity is an intrinsic property 
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of the channel and does not depend on the message size or on e. 

A similar scenario goes to the classical computation, namely computational threshold 
theorem. A circuit with N error free gates can be simulated with error probability less 
than e using 0(poly(log(iV/e).iV)) gates, each gate fails with probability p, as long as 
p < p t h, where p t h is independent of N and e. 

If the circuit is made with only a single element, with error probability p, then the 
processing fails with probability p. However, we may use recursive codes. Encoding, 
meaning that each gate is copied to three gates, and majority vote is taking on the total 
gates. This is the level-1 circuit and fails with probability of (2) = 15. Therefore the 
failure probability is ~ 15p 2 and if 15p 2 < p then the process is successful. This gives 
p < 1/15. 

For higher level, say level-2, repeat the encoding once more. As a result, it turns out 
that output is incorrect if two or more blocks, each block containing three encoded gates, 
fail. Then, pf aU ~ 15(15p 2 ) 2 = . If we make notation as the number of faulty paths 
to be represented by c = — , then 



level k : cpf aU = (cp) 
then 



pth ■ 

level I : cp f a a = (cp) 2 (1.38) 
level 2 : cp fa u = (cp) 4 

2 k 



PfaU — I P \ ( L39 ) 



Pth \Pth 

The code, which we used above is the recursion code. The computation is performed 
on encoded data, which is never decoded. By definition a procedure is fault tolerant if a 
single component causes at most one error in each encoded block of bits in the output. 

Clearly, to make a circuit fault tolerant requires increase in the size of the circuit for 
each recursion. If d stands for the size of the level-1 circuit, then the circuit size increases 
as ~ d k where k is the recursion level. Therefore, there is an exponential increase in the 
recursion for fault tolerant. However, a double exponential decrease in error probability 
shows that the situation is not so bad. More explicitly, the size of a fault tolerant circuit 
can be computed as follows. To simulate an iV gate circuit, with error probability less 
than e, each gate should have error less than e/N. Thus, the fault tolerant procedure for 
each gate must have failure probability satisfying the following relation 

\ 2 fc 

P 



• < "TF 1 (1-40) 

Pth J Np th ' 

the circuit size Nd k can be computed with solving the above inequalities 

2*log-£-~log-^, (1.41) 
Pth Np th 



then 

log NPth 

log 



V 
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Then we have 

log Npth 



NdX-Ni-^f^, (1.43) 

6 p 

which gives 

AT 

Nd k ~poly(log — ).N. (1.44) 

This proves that the circuit size grows polynomially for making the fault tolerant 
circuit. Of course there have been several assumptions such as free of noise wires. The 
assumptions may make the threshold smaller than the estimated one. 



1.4 Fault Tolerant Quantum Computing 

The issue of fault tolerant quantum computing is more critical because that a quantum 
computer is more fragile against noises. Coherent states are inevitably required all during 
a quantum processing, (Ef), jSj), (S3), JES), (O, d^B)- 

1.4.1 Quantum Error Correction Code 

Classical error correction codes might seem at a glance to be useful for introducing the 
quantum counterpart of error correction codes. However, there are several quantum 
essential properties preventing this idea. 

. Quantum states cannot be cloned. Therefore, encoding, meaning that making several 
copies of a particular state is not directly possible. 

. Errors are continuous. In addition to the bit errors there are also phase errors involved 
in a quantum system. 

. Quantum states collapse when are measured. Therefore, measurement can not give 
a proper way to repair the quantum state from the noises as the state is already 
destroyed. 

In 1995, it turned out that quantum error correction also is possible, (Eih . The objection 
of the quantum measurement destroying the state is resolved by an idea similar to what 
has been used for introducing the quantum algorithms. There is no need to measure the 
quantum state as far as we only need to know about the effect of the environment. Then, 
a measurement on the effect of the environment is sufficient and in this process the part of 
quantum state would be preserved. This is shown by the operator measurement. Actually, 
operator measurement is widely being used for defining measures of entanglement in 
addition to its usage for quantum algorithms. 

Assume a system, which is supposed to be measured. In a straight manner, one 
may measure the system directly, while in a rather elaborated scheme the system can be 
measured if we make the system close to the other system, any ancilla, and make them 
interact to each other and finally measure the state of the ancilla system. This process 
is shown in Figure H~H1 
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Figure 1.8: Quantum circuit for operator measurement. 

Given a unitary transformation U, with eigenvalues ±1 and eigenvectors {u^), the 
operator measurement with ancilla state initially |0), as shown in Figure fTTHj. can be 
described as follows 

(|0»(coK) + ciM) " (1.45) 

|0) + |1^ / I V , I \\ H 

(c |m+) + ci|w_» -> 



V2 

|0)(cq|m + ) + Ci|tt-)) |l)(cp|tt + ) - Ci|M_)) _ff 
y/2 y/2 

2 (c \u+) + Cl |«_» + 2 / (c |m + ) - Ci|«_» = 
Cb|0)|«+) + ci|l)|«_). 

Therefore, if the measurement result of the ancilla measurement is zero the is at \u + ) 
and otherwise for the measurement outcome being one the state of is \u-). By the 
means of operator measurement we have got information on the desired state without 
spoiling the state itself. 

The objections of continuous errors and no-cloning, both can be circumvented by 
the concept of entanglement. Entangled states make the errors orthogonal and then 
distinguishable. Also, entangled states play the same role as the redundant copies in 
classical error corrections. Here, we would like to have an overview on quantum error 
correction codes and fault tolerant quantum computing in order to see the essential role 
of entanglement for this concept. 

A quantum code C[[iV, K}} is a k qubit superposition of an n-qubit Hilbert space. For 
example C[[3, 1]] is a 3— qubit quantum code, which is explained as follows. 

3-qubit Bit-flip Code 

By definition, bit-flip error is the quantum operation defined as follows, (J3fih 

e(p) = (l-p)p + P XpX, (1.46) 

where p is a probability for a bit-flip. 

In a 3-qubit bit-flip error correction code, the state is encoded to 3 qubits. Logical 
is | Ox,) = 1 000) and logical 1 is = |111). Through out the circuit with bit-flip error, 
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the state is changed as follows 



a|0 L > + 6|l L ) 



final state 

000) + 6)111 

001) + 6|110 

010) + 6)101 

100) + 6)011 

011) + 6)100 

101) + 6)010 

110) + 6)001 

111) + 6)000 



probability, 

(1-p) 3 , 

p(l-p) 2 , 

p(l-p) 2 , 

p(l-p) 2 , 

p 2 (l-p), 

p 2 (l-p), 

p 2 (l-p), 

p 3 . 



(1.47) 



Using the operator measurement as introduced above with the unitary transformations 
Ui = IZZ and U 2 = ZZI, the bit-flip error that is occurred on each element is given and 
is called error syndrome. For each error syndrome, then proper recovery operation can 
be applied in order to recover the state from the effect of the noise. 



a|0 L )+6|l L ) -4 < 



State Ui 

a|000) + 6)111) 

a|001) + 6)110) 1 

a|010) + 6)101) 1 

a|100) + 6)011) 

a|011) + 6)100) 

a|101) + 6)010) 1 

a|110) + 6)001) 1 

ajlll) + 6)000) 



Ui Recovery Operation, 

I , 

x ± , 

1 x 2 , 
1 x 3 , 
1 x 3 , 
1 x 2 , 
x u 

/ . 



(1.48) 



As far as there is no information on a or 6, then the measurement results leave the 
state unchanged. 

Fidelity is a measure of distance between quantum states. The fidelity of states p and 

o~ is defined to be 

F(p, a) = trV pVVp 1 ^. (1.49) 



If one of the state is a pure state then the general definition is changed into 



F(\i;),p) = v ^P\pJp). 



(1.50) 



Clearly, the final outcome p is not perfect, because there is some probability of occur- 
ring two or more bit-flip error, that this code cannot correct them. Then, qualitatively, we 
measure the distance of the perfect corrected state and the state corrected with 3-qubit 
bit-flip code. 



n\i>),p) 



= y/l - 3p 2 + 2p 3 

= l-0(p 2 ). 



;i.si) 
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Figure 1.9: Quantum circuit for 3-qubit bit-flip code. 

Then the fidelity, even though is not perfect, still is better that 1 — 0(p), which is the 
fidelity for not using the code. The quantum circuit, which can realize the 3-qubit bit-flip 
error code is shown in Figure lL9l It is also shown that the above bit-flip error correction 
code corrects small rotation angle error. 

3-qubit Phase-flip Code 

Let the error for phase-flip channel be described as follows 

e pf (p)=pp + (l-p)ZpZ. (1.52) 

Correcting a phase flip error can be made in a same manner as the one for bit flip 
error. Recall that HXH = Z and HZH = X, so He p f(HpH)H = Ebf(p)- The encoding 
of the states is as follows 



+ ++), 



(1.53) 



\0l) 

\w = \-—), 

2. Then the syndrome measurement is performed by XXI and 



where |±) = (|0)±|1) 
IXX. 

The reason for the importance of the phase flip error is that the arbitrary qubit 
errors are proved to be combinations of bit-flip (X), phase-flip (Z) and bit-phase flip 
(XZ) errors. The quantum code, which corrects the arbitrary single qubit error is Shor's 
9-qubit code. 

9-qubit Shor's Code 

This code corrects any single qubit error. The encoding is as follows 

(|000) + |111))® 3 



|0l> 



;i.54) 



(|000) - |111))® 3 
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The syndrome measurement operators for bit-flip and phase-flip errors are 

Z\Zi, Z2Z3, Z A Z^, Z^Z & , Z 7 Z%, Z%Zq (1.55) 
X1X2 X3X4X5X6, X^X^X^X-jX^Xg 

It is simple, by direct examination, to check that any single bit flip or phase flip or bit- 
phase flip error gives a unique syndrome measurement result, then can be corrected. Let 
us give an explicit example. Suppose an error on fifth qubit. The initial state changes 
into the following state after the error Y 5 

a(|000) + |111))(|000) + |111))(|000) + |111>) + (1.56) 
6(|000) - |111))(|000) - |111))(|000) - 1 1 11>) -> 
a(|000) + |111))(|010) - |101))(|000) + |111>) + 
6(|000) - |111))(|010) + |101))(|000) - |111>) 

Then the syndrome measurement with the operators defined above, gives result as follows 



Operator 


Syndrome 


Z\Zi 










z^z^ 


1 


z^z^ 


1 


Z 7 Z% 





z%z§ 





X1X2X3X4X5X6 


1 


X^X^XqXjX^Xq 


1 



Or if a phase flip error occurs on the second qubit, Z 2 then 



Operator 


Syndrome 


Z\Zi 










ZaZ<o 





Z'oZq 





Z 7 Z% 





z$z$ 





X\ X2 X3 X4 X§ Xq 


1 


X^X^XqXtXsXq 






It is easy to check that the above syndrome is the same as if Z x or Z 3 happens. Then, 
this code is not an ideal code. 

We realize an easier way for finding the error syndrome. Use the commutator relation 
as [Zi,Zj] = and [Xj,Xj] = 0. Then, any Z or X error has a zero value through 
measurement of the syndrome operators which include the operators that commute with 
error operators. On the other hand, [JQ, Zj] is not zero if % ^ j. Then, the corresponding 
syndrome measurement would give one. This can be checked through the above examples, 
too. 

After finding the error syndrome, the proper operator is applied on the state in order 
to remove the effect of the noise and recover the state. 
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1.4.2 Quantum Error Correction Criteria 

Suppose a quantum channel as e = Ylk^kpEl. Then, let C be a quantum code defined 
by the orthonormal states {IVv)}- There exists a quantum error correcting e on C, iff 

. Orthogonality is satisfied: (ipi\E]E k \ipi) = 0, Vj + k, VI. This means that error should 
take code words to different points. 

. Non-deformation: (^i\E k Ek\t/)i) = d^il. This means that error should not be able to 
distinguish between the code words. 

We escape the proof here. Instead we explain the stabilizer and quantum fault tolerant, 
which we can see clearly the role of entanglement. We will see that these conditions will 
be met again by stabilizers. 



1.4.3 Stabilizer and Quantum Fault Tolerant 

Up to now, we have introduced classical fault tolerant and quantum error correction 
codes, which define the basis for being able to build a quantum fault tolerant containing 
system of quantum gates. However, we have not yet seen how to perform operations 
on encoded data and how to control the properties of quantum error correction inside a 
quantum circuit. These will be met by the method of stabilizer. 

Stabilizer is a mathematical method, however it is the Heisenberg representation of 
quantum computation, instead of the Schrodinger representation, which we are used to 
it and it works with wave functions. We first define the mathematics of stabilizer and 
then will use it to define quantum codes using stabilizers, (jlifih . (jsfih . 

Suppose a Pauli group. The Pauli group G n on n qubits is the group of all n-fold 
tensor products of {X, Y, Z} and {±1, ±i}. For example 

Gi = {±1, ±il, ±X, ±iX, ±Y, ±iY, ±Z, ±iZ}. (1.57) 

By definition, a stabilizer for a vector space V s with basis as (l^i)} is a se t 

S = {geG\ g\^) = \^)^)eV s }. (1.58) 

By convention —I S. Note that S is Abelian, meaning that all the elements of this 
group commute to each other because that they have the same eigenstate. 
For example, for the vector space as follows 

V,-^$P. d.M) 

the stabilizer is 

S = (XX, ZZ) = {XX,II,ZZ,YY}. (1.60) 
Another example is that the vector space for the following stabilizer 

S = {X,Z}, (1.61) 

is the empty space, clearly because that the defined stabilizer elements do not commute 
to each other. 
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For the vector space 

K = {|000>,|111>}, (1.62) 

the stabilizer is 

S={IZZ,ZZI). (1.63) 
As we had for the syndrome measurement. 

Another more interesting example, might be the one corresponding to the stabilizer 
for the vector space 

V a = {(000 + 111)® 3 , (000 - 111) 03 }. (1.64) 
The stabilizer for this case is 

S = {ZiZ 2 , Z 2 Z%, Z^Z 5 , Z 5 Zq, Z 7 Z & , Z 8 Zg, (1.65) 
X1X2 A3X4X5A6, X4X5X6X7X8X9} 

This is the familiar Shor code for error correction. Then it turns out that these stabilizers 
are actually codes, that have been introduced, part of them, as error correction codes. 

Stabilizer Code 

A [[n, k]] stabilizer code C(S) is the vector space stabilized by S G G n , where S is the 
set of stabilizers S = (g±, g n -k)- Here, n is the number of physical qubits and k is 
the number of logical qubits. The minimal generator for the stabilizer is (n — k). For 
example Shor's code [[9, 1]] has 8 generators. 

Given a certain set of errors {E a } G G, if there exists an element in the stabilizer, 
3g G S, such that Eg = —gE, for all E = E^Eb, then the set of errors {E a }, can be 
corrected by the stabilizer code, C(s). 

The proof comes in the following for all G C(s), but we emphasize that the 
explanation is also the same as the previously shown quantum error correction criteria. 

(V>|£|V> = WEgm (1.66) 
= -WgEty) 

= -(mm 
= 

This is equivalent to the quantum error correction criterion as explained by, (ipi\EjEj = 
0,Vl,i 7^ j. Hence, the above explanation says that the error can be corrected, but what 
can be the syndrome of errors? 

For S = (gig2--.gk), and error E, the error syndrome is 

'0 if [E, 9l ] = 0\l = l^k, 
1 otherwise. 

For example, consider the bit error, which has been explained before. The stabilizer 
S = (IZZ, ZZI). For error E = I XI, the error syndrome is {1, 1}, showing that an error 
has been occurred on the second qubit. However, the same stabilizer elements for an error 
E = IIZ, which is a phase-flip error, gives the error syndrome as {0, 0}. Then, phase flip 
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error is not an error or in other words is not a detectable error with this stabilizer. This 
statement is coincident with our previous outcome for the bit flip error correction code. 
Lets assume the stabilizer defined as follows 

S = {XZZXI, IXZZX, XIX ZZ, ZXIXZ}. (1.68) 

Clearly all the elements commute to each other as this should be an Abelian group 
to be a stabilizer. This is the 5-qubit quantum error correction. Suppose a bit flip error 
on the third qubit, E — X 3 . The syndrome would be {1, 1, 0, 0}. The error syndrome for 
a phase flip error e.g. E = Z 4 would be {0, 1, 0, 1}. This code corrects any {X i: Z h Yj} 
errors. 

5-qubit error correction code is called perfect code. The number of possible error for 
five qubits are ( 1 ) x 3 + 1 = 16, addition to one is for no error case. The number of 
syndrome strings is 2 4 = 16. Meaning that the number of syndrome strings is equal to 
the number of possible errors. However, this code has been known to have some bad 
property. It does not give the fault tolerant for some important quantum gates namely, 
Hadamard gate, as will be explained later. 

The other example which can be explained nicely with stabilizer is the Steane 7-qubit 
code. The stabilizer elements S is defined as follows 

gi = IIIXXXX (1.69) 

g 2 = IXXIIXX 

g 3 = XIXIXIX 

g A = IIIZZZZ 

g 5 = IZZIIZZ 

g % = ZIZIZIZ 

This stabilizer is a representation for Steane code, [[7, 1]], and corrects any single error. 
Encoded Operation; Normalizer 

The complete fault tolerant quantum computation requires the corresponding fault tol- 
erant operation, too. In this sense, unitary operation on the stabilizer state is called 
Normalizer. Suppose, a unitary operation U applied on the state \ip). Then 

W-W>- (1-70) 
If is a stabilizer state, g, then in the Heisenberg picture it is equivalently written 

as 

9 - UgUK (1.71) 

The proof is simple as follows 

U\*P) = Ugty) (1.72) 
= [Ugrf]U\1>) 

where [UgW] is the new stabilizer of the transformed state in the Heisenberg picture. 
We want to use this picture in order to understand how to perform operations on the 
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encoded state. Instead of following the state around, we follow the stabilizer around to 
have a fault tolerant quantum computation with fault tolerant operations on encoded 
states. 

For example, suppose the state \ip) = a|000) + 1 1 1 1 ) , which is the encoded state. 
For a fault tolerant quantum computation, operations should be applied on this encoded 
state. For example, for an operation X, then the encoded operation would be the logical 
X = X, which changes the state into = aj 111) + 6|000). For important quantum 
gates, the operations are as follows 

X : a|000) + |111) -> ajlll) + 6|000) (1.73) 

Z : a|000) + |lll) ->a|000) -6|111) 

9 9 
R Y (9) : |000) -> cos-|000) + sin -|111) 

a : looo)^ |000 > + ' 111 > 

v2 

But the question is; Are these fault tolerant operations possible? The requirement for all 
of these to be possible is that, the transformed state after the operation still should be 
an element in the stabilizer. This is explained as follows 

Vg e S,Ugrf e S. (1.74) 

Meaning that for any unitary transformation U the element of the stabilizer should 
stay in the same stabilizer. Satisfying this requirement limits the set of operations to the 
normalizer, N(S), defined as follows 

N(S) = {g e G\ ghg^ eS,Vhe S}. (1.75) 

The state may changes into other elements of the stabilizer but would not leave it. 

Now, recall that gh = ±hg, so ghg^ = ^hgg^ = ±h, but —h S, because by definition 
— / ^ S, thus gh = hg. Therefore, for any normalizer g, [g, h] = 0, for all the stabilizer 
elements h. As a result, the important fact is that the normalizer of the Pauli group is 
generated by the following elements 

N(G) = (H,S,CNOT), (1.76) 

where S = y/Z. This group is called Clifford group and includes all the operations that 
are involved in the stabilizer. However, the group does not include some important gate 
such as Toffoli gate or any arbitrary operation. Therefore, the group is not universal 
for classical computation, because of not including Toffoli gate, and then for quantum 
computation, as it is the subset of the classical computation. 

For example, consider the 5-qubit code, with the stabilizer elements as explained 
before. It is possible to check that for the normalizer elements N(S), X = XXX XX 
and Z = ZZZZZ. However, it is not possible to have the same for Hadamard, namely 
H ^ HHHHH . This is the bad property for the 5-qubit code, as mentioned before. 

Now, recall the Steane 7-qubit code. It is possible to check that in addition to the 
logical X and Z, Hadamard gate also can be written as H — HHHHHHH . So, all 
the Clifford group can be normalized with Steane 7-qubit code. This is the important 
property for this code because the operations are important for fault tolerant quantum 
computation. 
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Complete Fault Tolerant Quantum Computation; Entanglement 



In order to have a complete fault tolerant quantum computation, there are two main 
requirements which should be satisfied. Firstly, all the operations should be performed 
on the encoded data. Meaning that there should be a way to construct the fault tolerant 
operation for any quantum gate. Secondly, the error propagation should be controlled in 
a proper manner. 

In order to satisfy the required conditions and to get the correct fault tolerant quan- 
tum computation, we see that entanglement is inevitably required. Suppose the first 
requirement, namely the universal normalizers. Universal set of gates is possible if we 
could make the fault tolerant Toffoli gate. As this is a 3-qubit gate, so it is difficult to 
work with it. Lets make it rather bit simple, by introducing the T gate as follows 



Therefore, it is possible to construct the Toffoli gate from T gate, if some non-unitary 
steps, such as measurement, are involved. 

The important fact, which has been shown by Shor, is that T gate can be implemented 
using Clifford gates, measurement in the standard computational basis, Z, and prepared 
entanglement. Then, entanglement makes the fault tolerant computation possible. The 
quantum circuit for making a fault tolerant T gate is shown in Figure IT. 101 It is clear from 
the Figure that, the process to make a fault tolerant T gate is so close to teleportation, 
but slightly different. It is computation by teleportation. In this sense, we see very nice 
unification of the concepts which have been explained up to here, and the unification 
point is entanglement. 

The other point, which has been mentioned to make a system fault tolerant is that 
to avoid the error propagation. The problem is serious as to how to avoid propagation 
of errors all during the computation and more seriously in the operational measurement 
or syndrome measurement. This is possible to solve the problem by making the bit- wise 
gates or introducing some ancilla in cat state. However, the more extended explanations 
may be required which would not be presented here but can be found in related papers. 

1.5 Conclusion 

We, on purpose, made this chapter in two parts. First part includes classical complexity 
versus quantum complexity. Some quantum algorithms are revisited. Quantum proper- 
ties such as superposition and entanglement turn out to be prerequisite. The attempts 
for introducing quantum algorithms have been focused on attaining an exponential en- 
hancement over classical scheme. In other words, We want a quantum computer to solve 
polynomially those problems which do not have any up to date known polynomial so- 
lution with classical computers. However, polynomial complexity of classical problems 
still might be more complicated. In this regards, measurement based model has been ex- 
plained followed by parallelisibility. Exactly speaking the scheme is much closely related 
to entanglement. Therefore, entanglement turns out to be important in this regard, too. 




(1.77) 
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Figure 1.10: A quantum circuit for making a fault tolerant T gate. See the text for the 
connection between T gate and Toffoli gate. The first qubits are each encoded to seven 
qubits. The quantum operations CNOT and Hadamard are in the Clifford group and then 
are fault tolerant. The applied T gate on the third encoded qubit is the logical T gate, 
not fault tolerant, but has been checked to be almost fine with some reasonable error. 
The non-unitary parts of the circuit, namely measurements, are inevitably required. The 
circuit is so close to the quantum circuit for teleportation, but one of the two conditional 
Z operation has been changed by X and S operations, both being in the Clifford group 
and then being fault tolerant. The outcome of the circuit is the fault tolerant T gate. 
The important point, which we would like to emphasize here is that the initial second 
and third qubits are prepared as to be entangled states and clearly the success of the 
process is based on the involved entangled states. 

In the second part, we studied fault tolerant computation. Quantum system is much 
fragile and if it is going to be used for any computation, stability of the system should 
be considered with high priority. This is a very nice result that quantum fault tolerant 
computation is possible even with a very strict threshold. The important point, as to 
our interest, is that in order to introduce a fault tolerant quantum computation, entan- 
glement should be inevitably involved. Therefore, we conclude from this chapter, as a 
tour on quantum computing, that entanglement is required for making a reliable quan- 
tum computer, which can implement real quantum operations in order to truly give the 
desired advantages, far beyond any classical processing. 
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Chapter 2 



Quantum Entanglement; 
Characterization and Measurement 



From the time around 1900, when quantum physics has been introduced, superposition 
of states has been known as the main quantum property. Namely, for a definite physical 
property, there are probabilities for the state to be at any possible physical superposition 
of the states and the total state for the property is defined by a superposition of all the 
possible states. This fundamental postulate of quantum physics is originally out of the 
classical concepts, in a sense that in classical physics the state of a physical property is 
assumed to be well defined even prior to the measurement. However the interaction or 
the measurement is the process which makes a quantum system to collapse to one of the 
state even though before the measurement, the state is ought to be in a superposition of 
all the states each might be given with a particular probability during the measurement. 

This is then a kind of paradox for quantum physics. How would it be possible to prove 
the concept of superposition, if we accept that a quantum state just collapses to one of 
its eigenstates with some particular probability and we are not able to extract all the 
possibilities through measurement of the quantum state. One may claim that a quantum 
state originally had been in the state which is detected through the measurement and 
there would not be any practical proof. Then, there would not be possible to verify the 
concept of superposition. 

In a sense, entanglement is a keynote in this regard. The paradoxical behaviour of 
some particular quantum states, namely entangled state, has been primarily pointed out 
by Einstein, Podolskey and Rosen (|37h. Later, Bell (0) intelligently gave an approach 
to detect how entanglement can be verified through measuring a quantum system with 
a physical apparatus. Then, entanglement actually could make it possible to check the 
quantum fundamental concept, superposition. 

We try to get close to the concept of entanglement as it is required for quantum infor- 
mation processing and quantum computation. There are some very important questions 
in this regard. Namely, what is entanglement? What is correlation? Are there some 
check points to detect or even to quantify entanglement? 

Entanglement is known as a resource. This is important to know the reason on why we 
give such an outstanding role for entanglement. This can be explained briefly by giving 
an example. One very important topic in quantum information processing is quantum 
teleportation (|27l). There are so many things that one can do with teleportation. Apart 
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from its application for sending quantum information, teleportation has been known 
to be generalized to some other concepts. Namely, quantum cryptography (jsol ): sending 
information in a safe manner, quantum computation and quantum error correction, which 
is discussed in chapter 1. These facts show that with having only teleportation we can 
do lots of things and this is in a sense appreciating the concept of entanglement. Since, 
teleportation is a pure quantum processing and up to now there is no classical counterpart 
for it. Teleportation is based on entanglement. Then, having entanglement, one can make 
teleportation with which she can do lots of things and in this regard entanglement is a 
kind of physical resource. 

Hence, entanglement is introduced as an important property among quantum states. 
In fact, we need to go further now and to find it out a way to work with this quantum 
property. How to detect and how to measure it quantitatively, are important topics, 
which are open questions still for general cases. However, for special classes of states, 
there are some very excellent papers and remarkable results, which we review here before 
introducing our new measure of entanglement. We will show a new entanglement measure, 
that is a proper case specially for the physical system that has been considered for our 
experiments on realizing quantum entanglement. In other words, we keep reviewing 
the known measures of entanglement and make them step toward introducing the new 
measure. As, we claim on the importance of entanglement, our attempts would be on 
realization of a true entanglement. Then, the realized entanglement in the laboratory 
would be measured with our homemade entanglement measure. The reason is that it is 
introduced in a way that it is specialized for our experiment then it is considerably easier 
to be measured. 

We first start with general explanations in order to get more familiar with entangle- 
ment. The starting point would be entanglement in bipartite system. Later we will make 
the system more complicated but we will see that there would not be any complete or 
perfect explanation for a general case. 

2.1 Entanglement in Bipartite System 

Suppose that we have two identifiable systems, A and B, with finite dimension Hilbert 
spaces, Ha = 4> d an d Hb = 4> d - In particular, Ha = 4> 2 an d Hb = 4> 2 , for bipartite of 
two qubits. 

Here, we would like to study the entanglement for a general bipartite system. Firstly, 
we consider the problem of characterization or detection of entanglement. Later, we will 
go around the problem on how to introduce a valid measure of entanglement and will 
give some examples. 

2.1.1 Characterization of Entanglement; Pure State 

The given bipartite state may be a pure or a mixed state. For each case, separately we 
study the concept of entanglement. 

For a bipartite pure state, by definition, a state is separable if and only if it is a 
product state. A pure state, \ip), is a product state if it can be written as 
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where \iP)a £ H,a and \i[))b £ 7~Cb- If the state is not a product state, meaning that if it 
is not a separable state, then the state is an entangled state. For example 

I0 + ) = ^(|OO) + |11)), (2.2) 

is an entangled state. For this particular example, it is trivial to check the impossibility of 
writing the state as a product of the general states (a|0) + (3\l)) and (a'\0) + /3'j 1)) , 
because it requires the following equations to be satisfied, simultaneously, 

^(100) + |H)) = (a\0) A + (3\1)a) ® (a'\0) B + P'\1) B ) (2.3) 

aa! = —= ad' = 

pa' = PP' = -^=. 

v2 

Then, the state \<p + ) is called an entangled state. 

However, in general case, this is not as easy as the above example. In practice, we 
need to find the right basis vectors to be able to write the state as a product state, and 
only if this is shown to be impossible to find such a product elements for a given state 
then the state would be proved to be entangled. Schmidt decomposition is useful in this 
regard. In other words, the number of terms in a Schmidt decomposition of a particular 
state, called Schmidt rank or Schmidt number, would be useful. 



Schmidt decomposition 

As far as one is dealing with pure states, Schmidt decomposition is a powerful way to get 
the information on the entanglement and also to calculate the amount of entanglement. 
Schmidt decomposition is a decomposition into the biorthogonal basis, and gives the 
smallest possible number of terms for a product basis. Generally, any vector \ip) in 
Ha ®Hb can be expressed as follows 

W) = £Ci,,K&;> (2-4) 

= J2\ k \a' k ,b' k ), 
k 

which is the Schmidt decomposition. The coefficients Cij is the elements in a cIa x ds 
matrix C. If C is a Hermitian matrix, then always it is possible, by using the eigenvectors 
and eigenvalues, to find the diagonal matrix D such that 

C = UDU\ (2.5) 

and to make the Schmidt decomposition. On the other hand, if C are not a Hermitian 
matrix still it is possible to calculate the Schmidt decomposition by the singular value 
method. Generally, it is possible to write 

C = UDV, (2.6) 
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U is a d,A x d,A unitary matrix and F is a 4 x (Ib unitary matrix. Therefore, the state 
can be written as 

d A d B ram(dA,d,B) 

= u ikdkkV k j\ai) <g) \bj) (2.7) 

1=1 j=l k=l 

min(d A ,d B ) / d A \ / d B \ 

k=l \i=l J \j=l J 

min(d A ,d B ) 

fe=i 

where A& = c4fc and we have used the new basis vectors 

/ d A \ / d B \ 

The coefficients are called the Schmidt coefficients. 

The Schmidt number is at most the minimum of d>A and ds, where dA and ds are the 
dimensions of subsystems A and B. But, when some of the Schmidt coefficients are zero, 
then the Schmidt number would be less than min(d^, d B ). 

In a physical point of view, the Schmidt number gives the number of degrees of 
freedom. For a product state, there is only a single degree of freedom and then the 
Schmidt number is 1. Therefore, for a bipartite pure state, the necessary and sufficient 
condition for the separability of the state is that the Schmidt number being equal to 1. 

The important point is that, after writing the state based on the schmidt decom- 
position, then it is so easy to detect the entanglement and measure it. If the Schmidt 
decomposition is given then just by looking at the Schmidt decomposition it is possible 
to get the information on the entanglement of the states. If the Schmidt decomposition 
contains only one term then clearly the state is a product or separable state and then not 
entangled. On the other hand, if a state can not be written in a Schmidt decomposition 
of only one term, then the state is entangled. 



2.1.2 Characterization of Entanglement; Mixed State 

In a similar way, suppose there are two systems, say A and B. There is an interaction 
between the two systems, which may raise entanglement. The state of the system is 
described by the density matrix p. The task here for us is to detect the existence of 
entanglement. 

For a bipartite mixed state, by definition, the state is called separable if it can be 
written as a product state. The product state is given as follows 

£P^PP\ (2-9) 

i 

with pi > and YliPi = 1- For the case of the pure state we could use this defini- 
tion extensively as there is a well established approach to decompose the state, namely 
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Schmidt decomposition. In the case of mixed state, however, there is not such a generally 
applicable approach. 

For bipartite mixed states, there is a necessary condition for separability of states. 



Peres-Horodecki criterion, or PPT (40f) , (|4lh . which stands for positive partial transpose. 



PPT says, if the state is separable, then the partial transpose of the density operator 
with respect to one subsystem is positive. PPT is a necessary condition for separability. 

For a chosen orthonormal product basis {\viVj)} = {\vi) <g> \vj)} for the state p, the 
partial transpose p TB is defined by its matrix elements 

P^,nu = {v m V^\p TB \v n V u ) = Pmvw (2.10) 

For small dimensions, 2x2 and 2x3, the PPT-criterion is also sufficient condition 
for the state to be separable. However, we should say, that for a general bip artite mixed 
state, the situation is not still clear. There exist entangled PPT states, (42). 

The concept of entanglement for mixed state is richer compared to the previously stud- 
ied pure state. There are some mixed state entanglement, which cannot be interchanged 
to the Bell states. The states are called Bound entanglement. 



2.1.3 Measurement of Entanglement; General Axioms 

There are some axioms intuitively that are required for entanglement measures to be met. 
During the process to make any entanglement measure, usually it happens that some of 
the conditions are discarded. Therefore, for each one of the axioms there are several 
modified versions. We will see that while the modified versions are mostly satisfied for 
the entanglement measures on mixed states, for the case of the pure states the essential 
properties in their strong forms are mostly adequate. The conditions on pure state 
entanglement measures are those specialized from the ones on the mixed states. 

One of the most important of these conditions is that, entanglement, apart from 
the method for its quantification, should not increase on average by any LOCC, local 
operation and classical communication. This is called the criterion of the monotonicity 
under LOCC. Monotonicity is an essential requirement as it is even argued that in case 
of satisfying this condition, the other properties such as convexity and so on would also 
be satisfied. 



The axioms for entanglement measures are summarized as follows (|43l). (44) 



(EOa) E(p) = 0, if and only if p is separable. This is useful, though it is very strong 
condition. For example, Bound entangled states have zero distillable entanglement, 
but they are not separable. Therefore, slightly modified version is as follows. 

(EOb) E(p) = 0, if p is separable. 



(Ela) Normalization. For a maximally entangled state in dxd dimensions, e.g. 
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the amount of entanglement is quantified to be E(\(j)^)((j)^\) = logd. Here, 

d 

with the two sets \i) a and \i)b are orthogonal bases. 
(Elb) For a Bell state, the entanglement is £(|0 + )(0+|) = 1. 



(E2a) Monotonicity. For any LOCC operation A, E(A(p)) < E(p). 

(E2b) The equality in the monotonicity condition, A, E(A(p)) = E(p) is for the case 
that A is a strictly local operation. This means that A is a unitary operation or it 
only adds some extra dimensions. 

(E2c) When A is strictly local unitary operation, then E(A(p)) = E(p). 



(E3a) (Continuity) Suppose two sequences of bipartite states, {p n } and {cr n }, living 
on a sequence of Hilbert spaces, {H n }. For all these sequences which satisfy the 
condition \\p n — a n \\i — > 0, then 

E( Pn )-E(a n ) 



1 + log(dim7i n ) 

where, by definition || ... 1^ stands for the trace norm defined as follows 

\ 1 =tr(y/A!A). (2.13) 



(E3b) Continuity condition as explained above can be weakened if the state, e.g. p n is 
a representation of pure states. 



(E4a) (Weak additivity) For all states p and any number n > 1, then 

E(p® n ) 



= E(p). (2.14) 

n 
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(E4b) (Asymptotic weak additivity). The condition of weak additivity is just modified 
to get the asymptotic version. Given e > and a state p, there exists an integer 
N > 0, such that for all integers n < N we have 

^!!)_ £<£(p) <^!!) + , (2 . 15) 

n n 



(E5a) (Subadditivity) For all states p and a, 

E{p®a) < E(p) + E(a). (2.16) 

(E5b) For all states p and m,n> 1, then 

E(p® {m+n) ) < E(p® m ) + E(p® n ). (2.17) 

(E5c) (Existence of a regularisation) Finally, the weakest form for the additivity condi- 
tions, which is satisfied usually by any measure of entanglement is the one explained 
here. It states that for all bipartite state p, there exists a limit as follows 

E°°(p) = lim (2.18) 

n— »oo fl 

for which E°° is called the regularisation of E. 



(E6a) (Convexity) Mixing states does not increase entanglement. Namely, for all bipar- 
tite states p and a 

E(Xp + (1 - \)<t) < \E{p) + (1 - X)E(a), (2.19) 

for all < A < 1. 

(E6b) For any bipartite state and any pure state realization p = Pi ^ 
and J2iPi = l> tnen 

E{p)<Y,PiE{\A)(A\)- (2.20) 
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2.1.4 Generalized Measurement 



In general, a measurement is represented by a set of linear operators {Mi}, called mea- 
surement operators satisfying the completeness relation as follows 



Y J MlM m = l. 



(2.21) 



The probability for the outcome to be m, p(m), is 

p(m) = ^\MlM m \^). (2.22) 
After the measurement, with result m, the state of the system changes into 



(^\MlM m \^) 



(2.23) 



Generalization to the mixed state is also possible. If the density matrix representing 
a particular mixed state, corresponds to a pure state \ip), meaning that p = then 
generalization is as follows 

p(m)=tr(M mP Ml l ), (2.24) 

p(m) is the probability for the outcome m. Then the state changes into the following 
state after the measurement with the outcome m 



M m pMl 

P = 



tr(M mP Ml ' 



(2.25) 



When the density matrix is a representation for a mixed state, instead of a pure state, 
the generalization is still valid if we recall that 



P = ^Pi\A)(A\ = ^Pipi- 



(2.26) 



Therefore, the probability to get the outcome m from the measurement is calculated 
as follows 



p{m) = ^p t tr{M mPl Ml n 



(2.27) 



E 



Pipi 



Ml 



tr I M m 

= tr(M mP Ml). 
The new state after the measurement with outcome m would be 



(2.28) 



M„ 



^2 Pipi 



Mi 



= M m pMl. 
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The projective measurement is a special class of the generalized measurement, where the 
observable, V, is a Hermitian operator and can be decomposed according to the spectral 
form as follows 

V = J2 V rn\v m )(Vml (2.29) 

m 

where v m and \v m ) are the eigenvalues and eigenvectors for V, respectively. The probabil- 
ity for an outcome m is p(m) = \(v m \ip)\ 2 and the state after the measurement is e td \v m ), 
with an arbitrary phase 9. 



2.1.5 Measures of Entanglement; Pure State 

In principle, it is possible to construct any measurement, which converts an initial state 
to a desired state, with a certain probability. For simplicity, suppose that two parties 
Alice and Bob are given a Bell state generalization to a higher dimension would 

be also possible. They want to transform it into the state cos0|OO) + sin0|ll). Note, 
that this is almost general example as the final state would be transformed to a general 
state as cos9\uo,vo) + sm9\ui,vi), by local operations which do not change the amount 
of entanglement. 

To do the task, Alice performs a measurement described by the measurement operators 

\) and (*» 9 ° V (2.30) 
sin fly y cos 9 J v ' 

The measurement operators clearly satisfy the completeness relation. After the mea- 
surement, we may have two states, depending on the outcome, as follows 

|V>i)=cos0|OO)+sin0|ll), \i\) 2 ) = sin0|OO) +cos0|ll). (2.31) 

Therefore, by this measurement, the maximum entangled Bell state has been changed 
into a less or equally entangled state. If 9 = or tt/2, then the output state is a product 
state, which is not entangled. However, if 9 = 7r/4 or 37r/4 then the output state is a Bell 
state. 

The above example, shows generally how a maximum entangled state can be trans- 
formed to a less or equivalent desired state by means of measurement. Hence, suppose 
the dual problem. Meaning that the two parties, Alice and Bob, want to change the 
state l^i) = cos0|OO) +sin0|ll) into the Bell state by measurement. Then we need 
to define the measurement operators. The operator M 1 , defined as follows will give the 
desired state. 

M 1 = C(^|0)(0| + ->(1|), (2.32) 

where C is a constant. The measurement operators should satisfy the completeness 
relation. The new state, after measurement with operator Mi, is 

M 1 |^i) = C(|00> + |ll», (2.33) 

with a probability of 

p(l) = (^i|mJMi|^i) = C 2 ((00|00) + <11|11» = 2C 2 . (2.34) 
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The desired situation is maximizing this probability. However, because of the com- 
pleteness relation, M|M2 = 1 — MlMi, the maximum value should be restricted. In the 
valid range for 9, then we conclude that the constant value would be sin 9. Then, the 
measurement operators are achieved as follows 

M l = tan0|O)(O| + (2.35) 
M 2 = y/l -tan 2 0|O)(O|. 

The probabilities for the outcomes are 

p{\) = (^i|MfMi|^i) = 2sin 2 #, (2.36) 
p(2) = {ip 1 \M^M 2 \ip 1 ) = l-2sin 2 6>, 

and the normalized new states are 

Wi) = ^(|OO) + |ll)) = |0 + ), (2.37) 

Wi) = |oo>. 

Therefore, there is a certain probability to change the partially entangled state into a 
maximally entangled state. The probability depends on how strong is the entanglement 
of the initial state. As much as the initial state is entangled, sin 9 = for entangled 
state, it is a larger probability, p(l) = 1 for entangled state, to change the state into a 
maximally entangled state. 

To conclude, we may use the above perspective in order to give the measures of entan- 
glement. The distillable entanglement E d , sometimes called entanglement of distillation, 
is defined as the maximum yield of Bell states that can be obtained, optimized over all 
possible LOCC. On the other hand, the entanglement cost, E c , is defined as the minimum 
number of Bell states needed to create a given state by means of LOCC. It should be 
clear, that two separated parties cannot prepare an entangled state if they can only com- 
municate classically. But if they have some entangled state, they can convert them to the 
desired entangled states. These definitions should be now clear after the above general 
examples on how to make the both side transformations, from Bell states to partially 
entangled states and vice versa. 

The above definitions are given for finite regimes. However, calculations of Ed and E c 
are hard for finite case. Only, intuitively, it is known that 

E d < E c . (2.38) 

The reason is clear. If it would not be correct the above inequality, then it would be 
possible to create entanglement by means of LOCC, simply by converting Bell states to a 
state not satisfying the above inequality and then converting them back to the Bell state. 
In this process entanglement would be increased(!). 

In the asymptotic limit, the corresponding definitions are 

E D (m ee E?M)) = lim Ed ^ n \ (2.39) 

n^oo n 

EM)) = E^(\^))= lim c{m } . 
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It is shown that for the asymptotic version, both the distillable entanglement and entan- 
glement cost are equal to the entropy of entanglement, E D (\ip)) = Ec(\ip)) = E E (\ip)) : 
defined as follows. 

Recall the example of the Bell basis, \<f) + ), on a composed system of two subsystems 
A and B. The information on the total state does not verify the information on each 
individual spaces A and/or B. This is clear since p A = p B = I. 

The fact, as mentioned intuitively with the above example, would be useful to get 
some idea about the measure of entanglement. Entanglement is very much related to the 
mixedness of the reduced density operator. Since, we know how to measure mixedness; 
typically by entropy; entanglement also would be measured accordingly. 

Therefore, entropy seems to be a proper way for quantifying the entanglement. En- 
tanglement of is defined 

Ety) := S(p A ) (2.40) 
= -tr[p A \ogp A \ 

= -^Pk^ogp k , 

k 

where S(p A ) is the von Neumann entropy of the state p A . Von Neumann entropy for the 
state p is defined as follows 

S(p) = -tr(plogp), (2.41) 

we have also used 

AU = X>I*X*I> (2-42) 

entropy is written based on the eigenvalues. 

For example, consider the state = \ab). Then p A = \a)(a\, which gives E(ip) = 
— 1 log 1 =0. On the other hand, for 1 we can simply calculate i?(0 + ) = — 1/2 log 1/2 — 
1/2 log 1/2 = 1. This value is actually the much that we can have for an entangled state. 

Therefore, after finding out the Schmidt coefficients entanglement is also quantified 

as 

£(V) = -5>>gA^, (2.43) 

Finally, we conclude that if the state is a bipartite pure state, there is a well established 
way to detect the entanglement and also to quantify the amount of the entanglement. 

2.1.6 Measures of Entanglement; Mixed State 

In contrast to pure state, for mixed states, there is no unique way to quantify entan- 
glement, even for asymptotic cases. For pure state, Schmidt decomposition turned out 
to be a very powerful way for detection and measurement of entanglement. Schmidt 
decomposition posses an excellent property that it gives the minimum number required 
for decomposition of states. However for mixed state this is not always true. 
For example suppose an equal mixture of Bell states as follows 

P=\(\<P + )(<P + \ + \<p-) (rl)- (2.44) 
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Each the component represented as the Bell state is entangled, however the total state 
would be also equivalently written as 

p = I(| 00 )(00| + |ll)(ll|). (2.45) 

In this form, it is clear that the state is not entangled, therefore, even though there is a 
possibility to write the state based on entangled state, but this possibility does not posses 
the entanglement of the state. Actually, as far as there is some inexpensive, separable, 
way to make a state, the total state would be generally separable state valued. Even 
though one may try expensivly to construct the same state, but his try would be no more 
appreciated. 

The above description can be emphasized in mathematical view point, too. The state 
p is separable, not entangled, if 3p k > 0, ^ fc = 1, \a k ) G Ha, £ Hb, such that 
p = ^2 \k\o>kbk) (flkbk\- If there is not any pk satisfying the previous statement then the 
state is entangled. This means that all the possible constructions to make the state should 
be checked and this task in general is not possible. Therefore, intuitively, we see the reason 
on why there is not a general and easy approach for measuring the entanglement of a 
mixed state. 



Entanglement Cost and Distillable Entanglement 

The entanglement cost E c is the minimum number of Bell states required to produce the 
state by means of LOCC. The distillable entanglement is the maximum number of 
Bell states that can be distilled by an optimal LOCC distillation protocol. Because that 
a perfect conversion between mixed states in the finite regime is generally not possible, 
therefore for mixed state the asymptotic versions, Ed and Ec are considered. It is shown 
that for any measure E which satisfies certain natural conditions for asymptotic measures, 
the entanglement cost and distillable entanglement provide upper and lower bounds 

E D <E < E c . (2.46) 

For Bound entangled states, the situation is clear as Ejj = and Eq is finite. Gen- 
erally speaking, based on the axioms that each entanglement measure should satisfy, the 
upper and lower bound is defined. 

i) If the entanglement measure satisfies (Ela), (E2a), (E3a) and (E4a), for all the states 

p, then 

E D (p) < E(p) < E c {p). (2.47) 

These conditions are very strong. Up to now, only the "squashed entanglement" 
(jisl ) has been known to satisfy them all. 

ii) For an entanglement measure satisfying (Ela), (E2a), (E3a) and (E5c), then for all 

states p, the regularized version E°° is bounded 

E D (p)<E co (p)<E c (p), (2.48) 

these conditions are easier to be satisfied and they are satisfied by "entanglement 
of formation" and "relative entropy of entanglement" . 
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Distance Based Measures 



The class is introduced with properties that satisfy the conditions (EOa), (E2a) and 
(E2c), namely the measures that are included in this class would be zero if and only if 
the state is separable, the expectation value does not increase under LOCC and they are 
left unchanged under local unitary operations. 

Suppose that D(p,a) is a distance function representing the distance between two 
states p and a. Here, a G S(TC), with S(TL) to be the set of separable states. Then the 
entanglement of the state p is defined to be as follows 

E(p)= m£D(p,a). (2.49) 

<t£o (H) 

Relative Entropy of Entanglement: The quantum relative entropy, S(p\\a), can be 
used as the distance function D(p, a). 

S(p\\a) = tr[p{\ogp- log a)]. (2.50) 

The required conditions to be satisfied by D(p, a), now should be satisfied by S(p\\a). 
Relative entropy is not a metric and it is not even symmetric. But still it has some good 
properties. It is nonnegative, and zero for identical density operators. The same unitary 
operations on both states leave it unchanged. Then it can be used as a distinguishablity 
between the states. 

Suppose that we are given a large but finite number of n quantum states, all the same, 
being either p or o. Our task is to find the state by measurement. The probability for 
inferring from optimal measurements on the composite system, that the given state is p, 
while it is indeed a is (j4fih 

P n (a -> p) = 2- nS ^ a) . (2.51) 

It is clear that this expression is not symmetric. Relative entropy can be calculated 
from the eigenvalues and eigenvectors of density operators. Suppose that and |rj) are 
the eigenvalues and eigenvectors of the state p and Sj and |sj) are the eigenvalues and 
eigenvectors of the state a. Then generally we have 

s (p\W) = ^{r^logri - (si\p\si) log^}. (2.52) 

i 

If the eigenvectors for the two density matrices are the same, i.e. when p and a 
commute, then rearrangement of the eigenvectors gives result as follows 

S , (p|| < 7) = ^r i (logr i -logs i ), (2.53) 

i 

that is much easier for calculation point of view. 

Then relative entropy of entanglement is defined to be as follows (j4 

E r (p)= miS(p\\a). (2.54) 

Relative entropy of entanglement satisfies the condition (EOa), because that S(p\\o-) = 
only when p = a. It also satisfies the normalization conditions (Ela) and continuity on 
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pur e states (E3b). It is nonincreasing under LOCC, satisfying the condition (E2s) (J4 

63), (B, 0), Q. 

For the regularization of the relative entropy of entanglement we have 

E R (p) = E? = lim -Al—L. (2.55) 

n^oo n 

The above measure also satisfies automatically (EOa), (Ela), (E2a), (E5a) and (E6a), 
as these conditions extend to any regularization of a measure that satisfies them (|51l). 

In order to generate other distance based measures, the quantum relative entropy can 
be replaced by other forms of distances. In the following part, we have listed some of 
these measures and the properties as well. 

Bures metric based measure: The metric is defined as (jH3), (|H3), (0), (J!3) 

D B (p\\a) = 2-2^/F(p,a), (2.56) 

where F(p, a) is called Uhlmann's transition probability, defined as follows 

F(p,a)^[tr{(^pV^) 1/2 }} 2 . (2.57) 

The corresponding entanglement measure satisfies (EOa), (Elb) and (E2a). 
Hilbert- Schmidt distance based measure: The distance is defined as (|47l ). ([55T ) 

D HS (p, a) ee ||p - <x||^ 5 = tr[(p - a) 2 ]. (2.58) 



however, it has been shown by Ozawa (J56j) that the above distance does not satisfy the 
sufficient conditions for a nonincreasing function under LOCC. But the distance defined 
as follows is shown to have the required properties (JoTf ) 



D T (p, a) ee ||p - a\\ = tr[V(/> - <r) 2 ]. (2.59) 



Entanglement of Formation 

Recall that in the asymptotic limit, the entanglement cost for pure state is given by the 
entropy of entanglement. In a sense, entanglement of formation is extension of this effect 
to the mixed states (|58f). (|59f). Then, for a given ensemble of pure states, e = {pi, 
the entanglement of formation is defined as 

£/(£) = ^PvMlV'i)), (2-60) 

i 

which is actually an average of the entropy of entanglement for all the states in the 
ensemble. 

A mixed state can be realized by a multitude of pure state ensembles. The entan- 
glement of formation for different ensembles might be different to each other. But, by 
definition, we consider the most economic ensemble, meaning that the one which is mini- 
mized among all the possible cases. Then the definition of the entanglement of formation 
for mixed state is as follows 

E f (p)=infJ2PiE E (m, (2.61) 
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where the infinimum is taken over all ensembles e = {pi, that realizes the state p. 

Speaking about the conditions which can be satisfied by entanglement of formation, 
the definition of the entanglement of formation can be useful to show the satisfaction of 
several of the conditions. Namely, (EOa), (Ela), (E2a), (E3a), (E5a) and (E6a), all are 
shown to be satisfied. Additivity has been shown for some _special cases to be satisfied, 
however for a general case still it is not proved (0), (jfHI), (0), (JH), (0)- 

In the definition of the entanglement of formation, we see the inifinimum which makes 
calculation of the amount of the entanglement of formation hard. However, for special 
case of two qubits it is possible to calculate it by using the concept of concurrence (0), 
©. 

By definition, spin flip for a pure state \ip), is defined as 

$) = Y\r), (2-62) 

where Y is the Pauli operator, \ip*) is the complex conjugate of \xjj), when it is expressed 
in a fixed basis such as (|00), 1 01) , |10), |11)}. For a general p spin flip is 

p= (Y®Y)p*(Y®Y). (2.63) 

The Hermitian matrix R is 



R^^VPPVPi (2-64) 
with eigenvalues, Aj, in descending order, concurrence C is defined as 

C(p) = max{0, Ai - A 2 - A 3 - A 4 }. (2.65) 

Then the entanglement of formation is calculated as 

E f (p) = e{C{p)). (2.66) 

Here, the function e is defined as 



e (C)=/»(l±^0^), (2.67) 



2 

with h standing for the binary entropy function as follows 

h{x) = —x logx — (1 — x) log(l — x). (2.68) 

Then, entanglement of formation, for this special class of states can be calculated. 
There have been also attempts on generalizing the scheme, though the complete case is 
still an open problem. 



Negativity 

This entanglement measure is introduced with the aim of introducing a computable mea- 
sure (0). All the entanglement measures including distillation of entanglement and 
entanglement cost, have been introduced up to now, include in their definition some form 
of optimization and therefore are difficult to be measured. However, negativity is based 
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on the trace norm, which has been introduced before and can be calculated using stan- 
dard linear alg ebra packages. There are two quantities. The first quantity is negativity 
defined as (@), 0) 




(2.69) 



which is equal to the absolute value of the sum of negative eigenvalues of p B . The other 
quantity is logarithmic negativity, defined as follows 

£jvr(p)= log ll/^Hi. (2.70) 

However, this is not strictly monotone under LOCC, only it is for a subclass of LOCC. 

Some of the required conditions are satisfied by Af. It is zero for all the separable 
states, but it is also zero for all the PPT entangled states. Therefore, it satisfies (EOb) 
but does not satisfy (EOa). It does not satisfy the normalization criterion (Ela), as it is 
not one but 1/2 for the Bell states. It is monotone under LOCC but the conditions, of 
additivity (E4a) and subadditivity (E5a) are not satisfied, since the superadditivity is 
satisfied as follows 

Af(p ® a) = Af(p) + N(o) + 2M(p)H(o). (2.71) 

On the other hand, logarithmic negativity Ej^ satisfies some criteria, too. Similarly, 
it satisfies (EOb). It does not satisfy, monotonicity (E2a), nor convexity (E6a). But it 
satisfies strong additivity as follows 

E N (p ®a) = E N (p) + E N (a). (2.72) 

Then weak additivity (E4a) and subadditivity (E5a) are satisfied. 



Squashed Entanglement 

As we explained in the above sections, among the entanglement measures that have been 
introduced up to now, there is not a perfect one, in a sense that it satisfies all the required 
conditions. Squashed entanglement (j45T ) is the one which satisfies all the conditions. This 
measure is defined as follows 

E sq (p) = M{-I(A; B\E)\p ABE extension of ptoH E }, (2.73) 

where p = p AB = tr e{p ABE ), as the infinimum is taking over all extensions to a third 
subsystem E. I(A; B\E) is the quantum mutual conditional information defined as 

I(A; B\E) = S(p AE ) + S(p BE ) - S{p ABE ) - S{p E ). (2.74) 

Except for (EOa), which is unknown yet, and instead it satisfies (EOb). Squashed en- 
tanglement satisfies the other conditions in their strictest manner, namely (Ela), (E2a), 
(E3a), (E4a), (E5a) and (E6a). This measure of entanglement is new and still is being 
developed. 
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Witness Entropy of Entanglement 

This entanglement measure is based on the concept of entanglement witness (joll). (jo^h . 
Entanglement witness is a Hermitian operator, which has positive expectation values for 
all separable states p, tr (Wp) > 0, and a negative value for at least one entangled state 
a, tr (Wa) < 0. Therefore, entanglement is detected if a negative value is detected for the 
observable. Conventionally, the entanglement witness is normalized to have tr (W) = 1. 
The entanglement witness is called to be optimal, W op for a state a if 

ti{W op a) < tr{Wa), (2.75) 

for all entanglement witnesses W . 

The witness entropy of entanglement is defined to be as follows 

E w {p) = logp/^maxfO, -tr(Wp)], (2.76) 

where D and d are dimension of the total Hilbert space and the smallest of the dimensions 
of the Hilbert spaces of the subsystems. 

Except the additivity, this measure satisfies other conditions, namely (EOa), (Ela), 
(E2a), (E3a), (E5a) and (EGa). 



2.2 Entanglement in Multipartite System 

Entanglement also can exist among more than two parties. Then it is called multipartite 
entanglement. The situation is much more complicated compared to the bipartite case 
and there are so many open problems which are still under progress. Nevertheless, some 
of the essential concepts and entanglement measures can be generalized to multipartite 
system from bipartite one. 

Generally speaking multipartite entanglement complication, apart from its several 
reasons, is also because that multipartite system is composed of subsystems, called here 
partitions. Entanglement on the whole system does not pertain the entanglement on the 
individual subsystems. Though in some cases, it is pertained on each subsystem too. 
Then, the problem itself has an intrinsic complexity. 

For example consider the Greenberger, Home and Zeilinger, GHZ, state as follows 

|GHZ) ee —(|000) + |111)). (2.77) 
v2 

This state is the maximum entanglement between three subsystems. However, if any 
of the subsystems is traced out then the density matrix sharing by the two remaining 
parties is represented as follows 

p=i(|00>(00| + |ll>(ll|), (2.78) 

which is separable. Thus, the maximal entanglement on three parties as represented 
above does not include any bipartite entanglement. 
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Some of the notions from bipartite system can be extended to multipartite system. 
Generalization of separability condition gives a corresponding condition for the separa- 
bility of multipartite state. Namely, a state p is separable if it can be written as a convex 
combination of product states 

i 

with Pi > and YliPi = 1- 

A separable state may be arisen from ignorance of a subsystem of the total multipartite 
entangled state. Because, as it has been shown by a simple example above, an entangle- 
ment multipartite state may or may not be an entangled state if one of the subsystems is 
traced out. By definition, a multipartite state that loses its entanglement and becomes 
a separable state by tracing out a substate is called, multiseparable, in accordance to the 
separable state with the above definition, |7oh . 

On the other hand, the PPT criterion for separability simply extends to the multipar- 
tite entanglement. Because, if a multipartite state is a separable, then transposing any 
of the substates should satisfy the positivity. 

However, for measurement of the entanglement, it is not always possible to get the 
proper measure by extension of the bipartite measures of entanglement, in general. En- 
tanglement of distillation and entanglement cost have been introduced as two important 
measures for the case of bipartite state. However, the definition of these two measures 
are based on the Bell state, and Bell state is a bipartite state. It is not so clear how 
to extend the Bell state to the higher level of states. Therefore, these two entanglement 
measures are not easily applicable for multipartite state. Hence, there have been some 
attempts on this regards. For almost the same reason, entanglement of formation also 
can not be adequate as it is based on the entanglement cost. 

In order to define a measure of entanglement, in a similar manner to the bipartite 
case, we need to have some conditions which should be satisfied with the entanglement 
measures. It can be easily accepted that all the conditions have been introduced pre- 
viously for the case of bipartite entanglement can be applied for multipartite case too, 
except that the normalization condition which imply that measure should be normalized 
for Bell state is not applicable for a clear reason. 

Distance based measures of entanglement do not have any problem in their definitions 
to be extended to the multipartite case. But even in this form, the intrinsic meaning of 
the entanglement measure is vague compared to the bipartite case. The above example 
on GHZ state shows that bipartite entanglement is zero for the maximally three-partite 
entangled state GHZ. On the other hand, there are multipartite not maximal but en- 
tangled states with non-zero bipartite entanglement. Therefore, measurement on the 
entanglement of the multipartite state does not have a clearly known constitute on its 
subsystems. 

2.3 Continuous Variable 

Almost recently, entanglement has been considered in continuous variable regime too. Our 
study up to now has been on finite dimensional Hilbert spaces, though. Generalization 
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of the previous discussion to the continuous regime has not been yet completed. But it is 
easy to accept that for an infinite dimensional state, there are states with infinite amount 
of entanglement. Recall that the entropy of entanglement for a maximally entangled state 
is logci, then for the case where the state is infinite, it gives infinite entanglement (jTlh . 

0- 

The definition of separability in a continuous variable bipartite system is as follows 



PiA A) ® P?\ (2.80) 

where pi is a probability distribution. But in an infinite dimension, states only can be 
approximately described with a form as the above one. However, PPT criterion is still 
valid for infinite dimension. For some special class such as Gaussian state, it is also a 
sufficient condition in addition to be a necessary condition |7^). (74), |75l). |76l). 



(Ha). 

The subject of the multipartite entanglement is still an active field of research and 
more elaborated results may be represented in this field. 



2.4 Operational Entanglement Measures 



Barbara Terhal in an excellent review paper (78), addresses very important note as fol- 
lows, "Understandably, when considering more complex physical systems, the dividing 
line between what is entangled and what is not entangled, may become somewhat fuzzy. 
The guideline in deciding these matters, I believe, should be the question: Do we have 
an operational form of quantum entanglement? What resource does the particular state 
constitute in quantum communication and computation?" 

The concept of operational entanglement measure is also very important. Here, it 
might be better to explain it in the part for Bipartite entanglement measures. But, as 
far as operational entanglement measure turns out to be much useful for our case as to 
introduce a proper entanglement measure, we would like to explain it separately. 

Operational entanglement measure can be explained as follows: firstly, take your 
favorite application of entanglement. For example, superdense coding ^7$). Then, study 
how well or bad it performs for different sets of states, for example states involved in 
nuclear magnetic resonance implementation of the particular example here, superdense 
coding. Then, define an entanglement measure from this and check that it satisfies at 
least some of the conditions for entanglement, above all LOCC monotonicity (^3). 

We have used almost the same prescription in introducing a new entanglement mea- 
sure. Superdense coding has been selected and the states for NMR experiment have 
been studied. However, we have modified the scheme in a sense that we demonstrate our 
entangled measure as based on the concept of entanglement witness. We will explain it 
more in the next chapter after explaining the NMR system. 

One slightly similar entanglement measure in this regard has been constructed based 



on the capacity of superdense coding (|79l). by Hiroshima (80), (|8l[ ). An asymptotic 
entanglement measure for any bipartite state has been introduced as based on superdense 
coding capacity, that has been optimized with respect to local quantum operations and 
classical communication, LOCC. General properties for an entanglement measure also 
have been proved for this particular case. 
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Chapter 3 

NMR Quantum Information 
Processing and Quantum 
Computation 



Nuclear magnetic resonance, NMR, is a well established field in physics, originally, and 
in chemistry, widely. NMR technology goes back to 1940's and constantly has been used 
largely for variety of purposes. In chemistry, it is a major tool to study the structure and 
properties of solids, liquids and gases, (jiih. (84), (HH). 

There are several good reasons on the wide attentions to NMR as a means to realize a 
quantum computer, (|86f). (|88h. (fijl). In this chapter, we try to give some points on 

NMR as a fascinating physical system for realizing a quantum computer. We first start 
from the criteria which are required for any valid candidate for realization of a quantum 
computer. These criteria are known as DiVincenzo criteria, (|9oh. Then we will see NMR 
as a physical system, which naturally satisfies most of the required criteria. However, 
there are some critical bad behaviour for NMR, and we will explain them. Finally, we will 
prove that some NMR realization of quantum information processing, namely superdense 
coding, is also accessible with a classical computer. Therefore, it is vague to assume 
the corresponding experiment as a true quantum processing. Then we will introduce 
an appropriate way for detecting the entanglement, the scheme will be shown to have 
advantages over the conventional ways for detecting entanglement as it uses easy to access 
operators for NMR and determines the status of entanglement in the system, with only 
a single run experiment. 

It turns out that repairing NMR from this disadvantages might be possible, however 
may not be so simple, (91), (|92l). Therefore, in the next chapter we will try another 
physical system, electron nuclear double resonance, ENDOR, (|93h . ENDOR pertains 
most of the advantages of NMR system, while the existence of an electron spin also 
makes the required experimental conditions for constructing high spin polarized states 
much easier . 



3.1 DiVincenzo Criteria 

The minimal requirements for a quantum computer have been originally introduced by 
DiVincenzo. Even though after it, it has been some additional criteria added by others. 
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Here we list the original criteria, 



3.1.1 A Scalable Physical System with Well Characterized Qubits 

A set of physical systems is required, in order to give a respresentation of quantum bits. 
The physical system should have two distinct levels as to be a representation of a two- 
level quantum bit. The famous examples are spin-1/2 particles and polarized photons. 
Scalability also is a very important property, which should be satisfied by the physical 
systems. While a quantum computer with a number of qubits up to 50 explores 2 50 
computational paths in parallel, well out of reach for any classical computers, for e.g. 
factorization problem of a 400 digit number, a few thousands logical qubit are required. 
Hence, adding up the required number of qubits for error correction, the number is almost 
100 times more. Then, the problem of scalability is always very challenging for all the 
candidates of a quantum computer. 

3.1.2 A Universal Set of Quantum Gates 

Operations in a quantum computer are defined in terms of quantum gates. Therefore, 
it is very important for any realization of a quantum computer to be able to give a 
valid realization for quantum gates. At least, a universal set of quantum gates should be 
possible to be performed by the corresponding physical system. 

The time evolution of a quantum system is determined by its Hamiltonian. Therefore, 
in order to realize quantum logic gates, it must be some ability to control the Hamiltonian 
over time. In other words, the resulting time evolutions correspond to the computational 
steps of an algorithm. 

Theoretically attempts have been on simplifying the required steps and also to rear- 
range the requirements and to make them based on the existing terms in the Hamiltonian 
of the system. Also, experimentally, engineering of the Hamiltonian is of a very top pri- 
ority field of study. 

3.1.3 The Ability to Initialize the State of Qubits to a Simple 
Fiducial State 

This is a fact that if the input of a computation is random, the output is of little use. 
Then it is highly required to have the ability to prepare reliably a pure input state. For 
a classical computer it is so easy to make all the bits initialized, however for a quantum 
computer somehow it can be one of the most challenging requirements as we will see that 
this issue is so problematic for NMR quantum computation. 

3.1.4 A Qubit-specific Measurement Capability 

A computation is useful only if we can read out the result of the processing. Because of 
the postulates of quantum physics, it is impossible to get all the information about the 
state of a qubit by a measurement. However, quantum algorithms make it possible to 
extract the required information through the projective measurements in the basis. The 
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basis, can be specially the computational basis as |0) and |1) or any other set of basis, 
because that we can interchange the basis by applying unitary operations. 

This criterion sometimes seems to be challenging for NMR, because in NMR there is 
an ensemble state and all the measurements give an average on the state. However, it has 
been shown that almost all known quantum algorithms can be adapted to be implemented 
on (NMR) ensemble computer. 

3.1.5 Long Relevant Decoherence Times 

Long decoherence time, much longer than the gate operation time is required. Decoher- 
ence time of a physical system representing a qubit should be much longer than the gate 
operation time in order to have possibility to perform several required operations in a 
computational process. This requirement is important and in so many quantum physical 
candidates for realization of a quantum computer, it is actually very critical subject. 
However, in the case of NMR quantum computing and quantum information processing, 
we need not to be worried about decoherence time. Since nuclear spins pertain very long 
decoherence time usually seconds or order of seconds while the order of computational 
time is microsecond. 

3.2 Liquid State NMR Quantum Computing 

In this section, we study, very briefly, whether liquid state NMR satisfies the required 
criteria, as summarized in the above section. We will not go to the details and will 
emphasize those requirements which are important mostly for this contribution. Our 
interest is mainly initialization. This criterion causes some disadvantages for NMR. We 
will see that as far as spin polarization is very low and the initialization to a particular 
pure state is not possible for NMR, the credibility of the particular quantum processing, 
by means of NMR is not well approved. 

3.2.1 Scalable and Well Characterized Physical System of Qubits 

The basic system used in NMR quantum computing is an iV-spin system. The system 
is a molecule which involves iV magnetically distinct qubits. Usually the spins are half 
spins, though the higher spins are also workable. 

In a liquid NMR experiment, the sample including the specific molecules are dissolved 
in a liquid solvent. Sample is in a special tube, called test-tube, which is specially assigned 
for NMR. The tube is placed in a strong static magnetic field, Bq. We usually call the 
direction of the strong magnetic field the z direction. 

Because of the strong static magnetic field, and the spin behaviour of the particles, 
each spin precesses along the magnetic field. The Hamiltonian for the system of N = 2 
spins can be written as follows 

H = hu A I zA + hu B I zB + hu AB I zA I zB . (3.1) 

In the above equation, I z is the spin angular momentum in the z direction, = 
— £ = A, B, where 7g is the gyromagnetic ratio for spin i. The last term gives to 
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first order, the spin-spin coupling, because that spins are also coupled to each other via 
a through-bond electric interaction. 

In order to have some intuition about the values in the above Hamiltonian, let us 
give some numerical values. Suppose a molecule includes two carbon- 13, spin half, called 
spins A and B. NMR strong magnetic field is typically 9.3 Tesla. Then, the resonance 
frequencies would be y 1 ~ ^ ~ 100 MHZ. Suppose the magnetic field to be 11.7 Tesla. 
Then a proton in this field would precess with a resonance frequency of j- ~ 500 MHz. 
Typically, carbon-carbon coupling is >i ^ s - ~ 100 Hz. The proton-proton coupling is several 
Hertzs. 

Then in NMR, nuclei represent the two-level physical systems required for quantum 
computation. The characterization of the spins would be possible if we consider the 
resonance frequencies for each different spin. However, the weak point for NMR, in this 
part, is its scalability. The number of nuclei in each molecule is of course restricted 
and cannot be increased to infinity. There have been some attempts to resolve this 
disadvantage of NMR. The cellular automata has been proposed by Lloyd, (Joih . which 
might be useful in this sense. However, still the problem of the scalability, in its real 
meaning as the number of qubits should be possible to be increased extendedly while the 
computation should be a fault tolerant, for liquid state NMR still is an open problem. 
Here, the main attention might not be focused on this part as we will see that even 
for small number of qubits still there are some points which should be considered more 
seriously. 



3.2.2 Universal Set of Quantum Gates 

In a theoretical point of view, gates or operations on qubits are fixed and qubits are 
introduced to the gates. However, in NMR quantum computation, qubits are being 
represented by spins that are fixed in molecular structure and quantum gates are defined 
by pulses that are selectively applied on each spin. 

Generally speaking, spins are manipulated by applying much smaller radio-frequency 
field, rf, Bi in the x — y plane to excite the spins on their resonant frequencies uoi. More 
explicitly speaking, the applied pulse is an electromagnetic field with the strength Bi, 
which rotates in the transverse plane at u; r f, at or near the spin precession frequency u . 



The Hamiltonian of rf is, (|96|), (|97|), (|9 



H T f = —hwi[cos(uj T {t + (j))l x + sin^rft + (p)l y ], (3.2) 

where iO\ = jBi and <fr is the phase of rf field. 

The motion of a nuclear spin that is subjected to both the static magnetic field Bo and 
a rotating magnetic field Bi is rather complex, if it is described in the usual laboratory 
coordinate system. However, comparably easier, we find that in the rotating frame and 
at uj T f = u Q , the spin evolves under an effective field 

B = Bi cos(<p)x + B x sin(0)y. (3.3) 

By varying phase, 0, and the magnitude of B 1; the rotation angle and axis can be 
controlled. 

9 = -yB lP uj, (3.4) 
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where pu is called pulse width or pulse length. Then the state of the system is calculated 
to be as follows 

\ m = e^/ 2 |£(0)>. (3.5) 

Therefore, the effect of e iStax ^ 2 is just a rotation of the state about the x axis by the 
angle St, where St represents the integrated power of the applied resonant rf field. As 
far as , only the relative phase between pulses applied to the same spin is a matter of 
importance then by changing the phases of the rf field by 90°, we can similarly implement 
Y pulses. Typically, for the system introduced above proton pulse length is 5 microseconds 
for 9 = 90°. 

Rotation about the z axis is also possible after implementing any arbitrary rotation 
about x and y axes, since 

R z (0) = R x (-ir/2)R y (9)R x (n/2). (3.6) 

However, there are more elaborated approaches to implement an arbitrary Z rotation. 

Operations can be applied particularly on the selected spin, by applying rf field on 
resonance frequency with the particular spin. Other techniques such as pulse shaping also 
would be used for this purpose. 

Two-qubit gates are implemented by using the pairwise interactions between spins 
in the same molecule. There are two dominant systems of interactions; direct dipolar 
coupling, and indirect through-bond electronic interactions. Dipolar coupling is described 
by an interaction Hamiltonian of the form 

U% = ^^[Ia-Ib - 3(I A .n)(I B .n)], (3.7) 

where n is the unit vector in the direction joining the two nuclei, and I is the magnetic 
moment vector. However, dipolar interactions rapidly average away in a low viscosity 
liquid. 

Through-bond interaction, J-coupling, is an indirect interaction. The magnetic field 
seen by one nucleus is perturbed by the state of electronic clouds, which interacts with 
another nucleus through a Fermi contact interaction. The Hamiltonian is 

n J AB = hji A i B = hji z i z + ^[i+i- + i-i+], (3.8) 

where J is the scalar constant ; but generally a tensor. The flip-flop term is negligible 
for weak coupling or heteronuclear species. Therefore, up to a very good approximation 
it would be as follows 

H J AB = hJl zA l zB , (3.9) 
and for a duration of t, the evolution would be 

Vj(t) = exp[-i27rJI*I*f]. (3.10) 

Arbitrary single bit operations in addition to the CNOT gate make a universal set 
of gates. The ability to construct the arbitrary single bit operation has been already 
discussed. CNOT gate also is implemented in several ways, e.g. as shown in the followings 

CNOT AB = R yA (-90)R zB (-90)R zA (270 = -90)R zAB (180)R yA (90). (3.11) 
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This is shown as follows, (J3jJ 
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(3.12) 



that is the CNOT operation up to an irrelevant overall phase. 

Therefore, any arbitrary operation can be performed on the selected spin as far as 
we have known how to make a universal set of gates. However, the trivial operation of 
"identity" seems to be non-trivial in NMR, since evolutions of spins are always on. This 
can be done using a common NMR tool known as refocucing. 

The point for NMR quantum computing and quantum information processing is that 
NMR is an already rich field of research going back to 50 years ago. Then occasionally, 
the advanced techniques of NMR have been proved to be much useful for the quantum 
computational purposes specially constructing quantum operations. 



3.2.3 Initialization of the State 

The issue of initialization is much complicated in NMR quantum computation and quan- 
tum information processing. As much as we discuss the technology of NMR, we will see 
that inherently NMR suffers from the lack of a true initial state even though it has been 
some proposals for artificial realization of an initial state. The fact which we will discuss 
more extensively is that the problem of initialization for NMR makes it possible to prove 
the classicality of the realization of some particular NMR experiment, which has been 
claimed to be a quantum based processing. We discuss the issue of initialization for NMR 
by recalling some fundamental aspects. 

Liquid state NMR quantum computation and quantum information processing have 
been done conventionally and mostly at a thermal equilibrium condition. The thermal 
equilibrium state of the NMR system is described by density matrix p as follows 

where (3 = 

At room temperature, k^T is much larger than the energy differences of the levels, 
typically 10 5 larger. Then, the state is approximated, with a very good accuracy for 
approximation, as follows 

(3-14) 
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where N is the number of spins. For e.g. 
then, Q, 



AknT 



N = 2 and suppose for a moment that ua. ~ 4o>b, 
\ 

(3.15) 
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This is a mixture of all possible pure states in the computational basis. 

Recall the fact that the NMR observables are traceless. Then the identity term in 
the above equation is not observable in NMR. The thing which is observable is only the 
difference in the populations of different states. 

Microscopic state for NMR might be in a large variety of states, however as far as 
only the ensemble averages are measured, there is no way to distinguish between them. 

Notably, and we will discuss it more in details later, entanglement is a microscopic 
property for NMR system. Then it seems appropriate to discuss the entanglement of 
NMR only in case of the direct access to individual states but this is not possible in 
NMR. This may look as a serious paradox then how to establish the entanglement of the 
NMR states. One way to this purpose, which we have used for our proposal on detec- 
tion of entanglement in NMR, is based on the operational entanglement measures. Any 
quantum processing, which is based on entanglement would not be reasonably performed 
successfully with a system of separable states. This fact will be shown to be much useful 
for NMR. 



3.2.4 How to Solve the Problem? 

Pseudo-pure state generally has been introduced by different people in several ways. 
Preparation of any pseudo-pure state inherently requires some non-unitary operations, 
such as summation of different experiments. 

Pseudo-pure state in any form of its realization has not been validated perfectly for 
our purposes for several reasons. Firstly, the number of resources, molecules or experi- 
ments, required for making a pseudo-pure state scales exponentially with the number of 
qubits. Therefore, an exponential number of e.g. molecules are required to make signals 
readable. In the current NMR experiment, with an Avogadro number of molecules and 
spin polarization ~ 1CT 5 , up to 20-qubit is possible to be worked with. By increasing 
the number of qubits in a same experimental condition, would give a very weak signal, 
smaller than the standard deviation in reading the results. Therefore, signals cannot be 
detected. 

Pseudo-pure state involves generally non-uniform operations for its realization hence 
its evolution is not clearly approved. Finally and more importantly to sake our interest, 
it is naive to think that pseudo-pure state gives entanglement even with low polarized 
spins. Pseudo-pure state might seemingly look working well, but the large amount of 
mixed states, any way, still remains in the NMR system. Recall, the discussion on the 
measure of entanglement for mixed states. There might be some way to construct a 
separable mixed state from entangled states. However, as far as the final mixed state is 
separable, the state is a separable. 

It would be all perfect if it could have been possible to make pure states directly from 
NMR mixed states. Direct purification of the NMR mixed states at the ground state 
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requires very strict experimental conditions of the temperature down to milli Kelvin. 
Apart from the complexity of how to reach easily such a low temperature, the sample 
would be no more a liquid sample but would change into a solid one. Solid state sample, 
then needs more elaborated operations on it. Dipolar couplings would be important in 
the Hamiltonian and need to have special care on manipulations of spins due to the more 
complex Hamiltonian, whereas the strict experimental conditions also would not be stable 
easily. 

Therefore, instead of getting involved on preparation of the experimental condition for 
making a pure state, people have considered indirect ways. There are several approaches 
to this purpose, among them are optical pumping, para- hydrogen and dynamic nuclear 
polarization or DNP, 0), (0), (@). 

DNP is the transfer of polarization from the electron spin in a free radical to the 
nuclear spins. DNP has resulted to a very high spin polarization up to 0.7. Though the 
molecule was not useful directly for an NMR quantum computation. The fact used for 
DNP is that the magnetic moment of electron spin is ~ 10 3 larger than nuclear spin, 
then nuclear spin polarization can be enhanced by electron spins. There has not been 
reported any successfully demonstration of the scheme for NMR quantum computation 
and quantum information processing. 

Finally, we summarize that the initialization is a significant problem for quantum 
computation and quantum information processing by means of NMR. Pseudo-pure state 
has intrinsically problems, specially for the problem of entanglement, as will be discussed 
more in details in following sections. Enhancement of polarization is possible whether 
directly or indirectly. While the direct approach needs very much elaborated experimental 
conditions, none of the indirect several techniques has been perfectly successfully done 
for special purpose of quantum computation, yet. 

3.2.5 Readout 

Measurement in NMR is done with an rf coil placed close to the sample. In order to 
readout a spin, a readout pulse, i? x (90) is applied on the thermal equilibrium state. As 
a result the nuclei are tiped from the z axis into the x — y plane. Then, nuclei precess 
along the z axes and generate a magnetic induction signal known as the free induction 
decay, FID. FID is captured by the pickup coil. Then FID is Fourier transformed to 
obtain the spectrum. In the obtained spectrum, different spins, qubits, in a molecule are 
spectrally distinguishable via their Larmor frequencies, u>i, and the amplitude and phase 
of the different spectral lines give information about the corresponding spin states. The 
extracted information depends on the convention, typically a positive absorptive line is 
due to a |0) and negative absorptive line is due to the spin in |1). 

Suppose the state of the system at the start of the measurement is described by 
p. Then, mathematically the experimentally detected signal can be calculated from the 
transverse magnetization for a single spin (98), 

M = duhTr[p(ia x + a y )\, (3.16) 

here d is the volume density of the detected spin. Then, for a .fT-turn solenoid coil with 
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quality factor Q and area A, and magnetic flux 0, we have 



V(t) = QK d f t (3.17) 

= QK^-u MA 
at 

~ QK'yBfi e~ t/T2 {d'yhTr\p(ia x + ay)]}A. (3.18) 

The corresponding voltage is complex. Both the magnitude and phase are measured. 
Then the envelope of the FID decays as e _ */ T2 , and yields to a Lorentzian line after 
Fourier transform as follows 

I ™ (3 19) 

l + (co-cu ) 2 l + (cu-cu ) 2 ' 1 • ] 

and has a linewidth at half height of the following form 

The magnetic signal of a single nuclear spin is too weak to be directly detected. 
Therefore, NMR experiment should be done using a large ensemble of identical molecules, 
typically on the order of 10 18 , disolved in a liquid solvent. The system is an ensemble 
system, rather than a single iV-spin molecule. 

This raises to some difficulty Suppose an NMR implementation of a particular quan- 
tum algorithm. Then, the information should be given by the distribution of the random 
numbers detected through the measurement. However, signal coming from the NMR en- 
semble state is an averaged over the ensemble state. Average of random numbers clearly 
is useless. This problem has been solved by slightly modifying quantum algorithms in a 
way that the averaged result also gives the required outcome. 

The other problem is that a density matrix representing the state of an ensemble 
system is completely described if all the elements are detected and known. However, 
NMR measurement of the operator —il l x — I* selects specific entries in the density matrix, 
called single quantum coherence (SQC) elements. The SQC elements connect basis states, 
which differ by only one quantum of energy. Therefore, in order to get the total density 
matrix elements, several 90° pulses are systematically applied in order to rotate each 
spin around the x and y axes to make observable all of the terms in the deviation of the 
density matrix from the identity. The process to get the total density matrix through 
the measurement is called state tomography. In a general scheme, without any prior 
knowledge on the state, the state tomography involves many experiments, then it is 
impractical for experiments involving large number of qubits. 

In the following sections, we will show that a proper scheme to measure the en- 
tanglement of states is the one which gives the required information on the status of 
entanglement with a less number of experiments as compared to the state tomography. 
Even though, it is always possible, for small degrees of freedom, to measure the status 
of entanglement after state tomography, but the scheme which employs less number of 
experiments is of course more appreciated. 
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3.2.6 Decoherence 



The decoherence process of uncoupled nuclear spins is well described by combination 
of two phenomena: longitudinal and transverse relaxations (|98h. (jiih. (|85h. These two 
processes are closely related to generalized amplitude damping and phase damping, re- 
spectively. 

Relaxation of nuclear spins is caused by fluctuations in the magnetic field experienced 
by the spins. Whether the magnetic field fluctuations contribute to energy exchange with 
the bath or phase randomization depends on the time scale of the fluctuations. Roughly 
speaking, two processes are happened as follows. Fluctuations at resonance frequency 
of the nuclear spins lead to efficient energy exchange with the spins. However, slow 
fluctuations or fluctuations at zero frequency, give rise to phase randomization. 

The issue of decoherence generally is very critical for any realization of a quantum 
computer. However, NMR is very much appreciated on this matter. Therefore, we are not 
going to g et involved in this topic. The extended discussion might be found in literature, 
(josh . (0), (jiih . We only recall as stated above that two main relaxation processes have 
been explained as T\, spin-lattice relaxation time and T 2 , spin-spin relaxation time. They 
are measurable quantities e.g. by standard inversion recovery and Carr-Purcell-Meiboom- 
Gill pulse sequences, fj96h . 

Generally speaking T2 < T\. For liquids with low viscosity T2 ~ 7\ and for viscous 
liquids T2 <C T\, due to long range of spin-spin interactions. In common organic liquids 
T 2 ~ T\ ~ 1, while in noble gases 7\ is several hours. In a liquid NMR experiment, 
e.g. with carbon-13 labeled chloroform, proton and carbon represent T 2 < T\, of several 
seconds. 



3.3 What is the Problem With NMR? 

Through the above sections, we understand the reason that why NMR has been intro- 
duced as one of the most viable methods for realization of a quantum computing, liquid 
state NMR has b een u se d ex te nsive ly for i m plem enting even relatively complicated quan- 



tum algorithms, JmJ, fiflj), JlflJ, JlflJ, £ 



However, the status of NMR has not been completely approved compared to the other 
candidates, such as photons and ion traps, in that it treats an ensemble system composed 
of a large number of molecules. Under the curre nt exp er imen ta l con ditions, the states 
are mixed rather than quantum pure states, (jlOSl ). (jlOfil ). (1H2), fjlOSt ) . 

Accordingly, nearly all previous attempts to demonstrate NMR quantum informa- 
tion processing (QIP) and quantum computations (QC) relied on pseudo-initialization of 
highly mixed states. However, with or without employing the pseudo-pure states, math- 
ematical arguments refute the credibility of NMR QIP and QC except at very low spin 
temperatures by verifying separability of the states. 

The issue of entanglement of the states, has been the main focus of so many researches. 
Nevertheless, still the exact and unified solution has not been accepted widely. The 
separability of NMR mixed states, density matrix, is mathematically proved. But, there 
are also some discussions from the other side. For example, it has been claimed that 
a certain mixed state would be possible to be explained by several density matrices. 
More particularly, suppose a system of mixed states. Conventionally, we know that 
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all the density matrix representations for this system are equivalent. Therefore, if a 
density matrix among them all has been proved to be separable, then the system of 
mixed state also is separable. This is because that the different representations cannot 
be distinguished from each other by any physical means. However, there has been some 
discussions also to show that this can not be always true. The different density matrices 
can be distinguished by showing the fluctuation, through measuring some particular 
observables. Then separability of a particular density matrix representation for a certain 
mixed state does not verify completely the separability of the mixed state, 

One more reliable way to solve the discrepencies, we believe that would be to show 
that any particualar implementation of a quantum algorithm by means of NMR is also 
accessible for a classical system. Therefore, in this manner, there would be no quantum 
beneficial to use the NMR system because the corresponding processing would be classical. 
The exact proof for all the cases is not still done. However, particularly, we can show 
that an implementation of a quantum algorithm, namely superdense coding, is possible 
with classical system. 

We show that with NMR mixed states, even if the states are definitely separable, still 
there is a considerably high probability to detect signals with appropriate signs which 
apparently imply transfer of two-bit message even though we encode only one qubit. 
The reason for detection of the signals is proved to be because of the large number of 
molecules being involved in the NMR ensemble state. Therefore, a system composed of a 
huge number of molecules should be necessarily prepared in order to realize the transfer 
of two bits of information. Then, taking account the required number of resources, 
molecules, it is totally non-sense to claim on the quantum advantage over the classical 
counterpart. However, for an exact demonstration of SDC with NMR, in the sense that 
it can realize enhancement over the classical communication, nuclear spin polarization 
should be increased above a certain threshold, and this threshold coincides with the 
mathematical criterion for non-separability of the density matrix. 



3.3.1 Ideal SDC with Nuclear Spins 

Let us consider a pair of nuclear spins 7 = 1/2 and S = 1/2 placed in a static magnetic 
field -B , and suppose for a moment that the system is initially in a pure state l^o) — 1 00) . 
Then, the procedure of SDC, whose quantum circuit is described in Figure 13. lj is as 
follows. Firstly, the entangling operation is applied on the state of the nuclear spins. 
The entangling operation is composed of a Hadamard gate (H) on the I spin and then 
a controlled-NOT gate (U cn ) whose control and target qubits are the I and S spins, 
respectively. The quantum state \ipi) after the entangling operation is represented as 

\ipx) = U ent |Vo) = ^|(l 00 > + I 11 )) = \Pao), (3.21) 

where we have defined U en t = U cn (Hj ® I5). Here, |/3oo) is known as one of the four Bell 

states 

l/U - '^> + ™*> , (3.22) 

where z,x = 0,1 and x — 1—x. Now, suppose that the nuclear spin I is given to Bob and 
the other nuclear spin to Alice. Bob encodes a two-bit classical message zx by applying 
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Figure 3.1: A quantum circuit for superdense coding. Two nuclear spins / and S are 
involved. After the first Hadamard and controlled-NOT, the spin I is given to Bob and 
S to Alice. Bob encodes the spin / by applying the unitary transformation U zx and then 
sends it off to Alice. Alice has now the two nuclear spins in her disposal. She applies the 
second controlled-NOT and Hadamard gates and measures the spin magnetizations. As 
a result, she obtains the encoded message, meaning that two-bit message is transferred 
while only a single spin has been encoded. 

the unitary operation U zx = Z Z X X on the I spin and then sends off the encoded qubit to 
Alice. The effect of the unitary transformation by Bob is to toggle = |/3 o) m t° the 
other Bell state. That is, 

|^ 2 ) = U* x |A)o) = |/U- (3-23) 

Then, Alice applies the disentangling operation, which is composed of a U cn followed by 
H. The state after the disentangling operation is 

|V>3> = U B cll|^2> = \zx), (3.24) 

where we defined Ueeii = (H/®Is)U cn . Finally, she performs measurement of the resultant 
magnetizations, Z/ and Z s and extracts the results as (— l) z and (— l) x . 

If it was possible to execute all the above process in an NMR apparatus, then the 
signal intensities detected through the measurement in this ideal case of the fully polarized 
initial state would have been maximum and information on the choice of U zx could be 
extracted with a success probability of unity. In this case, the maximum entanglement is 
shared between Alice and Bob, which causes transfer of two classical bits of information 
from Bob to Alice via the I spin alone. Communication by SDC is efficient by a factor 
of two as compared to the classical communication. 

3.3.2 Transfer of Two-Bit Message by Encoding a Single Spin 
in NMR 

NMR is a spectroscopy, manipulating the nuclear spins placed in a strong, typically 10 
Tesla, static magnetic field B . The corresponding total spin Hamiltonian in an isotropic 
liquid state is composed of several terms, but the most important parts are considered 
as the Zeeman Hamiltonian and spin-spin interaction parts. So that we take 

TO— 1 TO — 2 

H = - + 4 E E (3.25) 

i=0 i=0 j>i 
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where Ui is the Larmor frequency and Zj is the Pauli operator for the z component of 
the ith spin and is the coupling constant between ith and jth spins in a molecule. 
In the high field approximation, which is valid for most practical cases, in the thermal 
equilibrium at temperature T, the polarization q is given by 

^ = tanh(^g^), (3.26) 

where 7$ is gyromagnetic ratio. Accordingly, the probabilities Pi and qi of finding the ith 
spin in the states |0) and |1), respectively, are given by 

Pi = — g— , % = — • ( 3 - 27 ) 
Then, the density matrix in thermal equilibrium is given by 

Pm-qubit = P 1 ® ® P" 1 = ®£ X P\ (3.28) 

where p l is the density matrix for the ith nuclear spin in a molecule and is written as 

P l : =Pi|0)(0| + gi |l)<l|. (3.29) 



At room temperature, the energy splitting between the ground state and the excited 
state is much smaller than the thermal energy and as a result, large number of spins 
(qubits) are in statistical mixture with low spin polarization. Accordingly, the information 
processing should be evaluated statistically. 

For nuclear spins J = 1/2 and S = 1/2 the initial state is a mixed state described by 
partially polarized density matrix with the form of 1)3.28)1 and (J3.29)) with m = 2 as 

Po = (pj|0)(0| + 9j |l)(l|)®(p s |0)(0| + ^|l)(l|) (3.30) 
= Pips\00)(00\+ Pl q s \01)(01\ + g/ ps|10>(10| + g /gs |ll)(ll|. 

After performing the entangling operation the state of the two qubits ()3.3())1 changes to 
a general Bell diagonal state of the form 

Pl = UentPoUL (3-31) 

= PiPs I Poo) (Poo \+Piqs\ P01 ) (P01 1 
+ qiPs\Pio)(Pio\ + qiqs\Pn)(Pn\- 

This is the shared state between Alice and Bob. Suppose that Bob owns the state of 
the nuclear spin I. Then, Bob applies the unitary operation XJ ZX . The state p2, after the 
encoding operation is still a general Bell diagonal state as follows 

p 2 = tWiti (3.32) 

= PlPs\Pz,x)(Pz,x\ +Piqs\Pz,x)(Pz,x\ 

+ qiPs\Pz,x){Pz,x\ + qiqs\Pz,x)(Pz,x\- 
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He hands over the nuclear spin to Alice, who applies the disentangling operation. Then, 
the state P2 changes into 

P3 = U B eiiP2U^ cll (3.33) 

= p I p s \z,x)(z,x\ +p I q s \z,x){z,x\ 

+ q I ps\z,x){z,x\ + q I q s \z,x){z,x\ 

= (pi\z)(z\ + qi\z)(z\) ® {ps\x){x\ + q s \x)(x\). 

Measurement of the spin magnetizations Z/ and Z$ is done on the total ensemble state 
composed of n molecules. Then, measurement on the product state ®f =0 p^ is composed 
of separate but simultaneous measurements on the spin magnetizations and is represented 
as Yu7=i ® ^5 Note, that stands for the density matrix p 3 of the iih molecule. 
Measurement of the spin magnetizations gives results as binomial probability distributions 
over (— n, —n + 1, —1, 0, 1, n — 1, n), with the mean values to be as follows 

M = (-l) z np I + (-iynq I = (-l) z ne I , (3.34) 
p s = (-lYnp s + {-lfnq s = {-lfne s . 

Let us make an assumption that z, x — 0, which does not influence the generality of the 
discussion. The corresponding variances are characterized by 

<7j = Anpiqj = n(l — ej), (3.35) 
a 2 s = Anp s q s = n(l - e|). 

Recall that oj and as are essentially the widths of the range over which the outcomes 
are distributed around the mean values ()3.34j) . The relative widths of distributions are 
characterized by 



pi nei ei^jn 

&s _ y/n(l - e|) 1 



(3.36) 



Us ne s e sy /n 

Then it is clear that the relative widths decrease as \fn with increasing the number of 
molecules n. Thus, the greater the number of molecules, the more likely it is that an 
observation gives a result which is relatively close to the mean values 1)3.34)1 . 

Now, it is also required to calculate the corresponding error probability of detection 
of a negative value for z = 0. For simplicity, let us consider only z\ , which is the 
measurement on the spin magnetization of the nuclear spin / in the ith molecule. It 
should be clear at the moment that results are applicable to the signals detected through 
measurement on the other spin magnetization in the molecule, Z^. The error probability 
is defined as 

n 

P e = P(J2 Z f < 0|z = 0). (3.37) 

i=X 

We would like to calculate the error probability to show that for a range of n of the 
current NMR experiments, this quantity is negligible. From the DeMoivre and Laplace 
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theorem, 

{53), 

if n is large enough then generally we have 



P <j a < £? z ? } ^ J < /3 j ~ _L ^ e ~4 da;. (3.38) 



For the error probability P e , therefore 



P e = P<{-oo< E ^ - <~> (3-39) 



e^^"dx 

oo 



1 -^r^J 

— e 2<T -r — 



2n /i/ 

1 (ne?) 1 



If the number of molecules are large enough, the result of the measurement on the spin 
magnetizations Zj and Zj on the ensemble state gives results very close to the mean 
values 1)3.34)1 with negligible error probabilities. In a conventional NMR experiment, the 
number of molecules are n ~ 10 18 , with low spin polarization e ~ 10 -5 . Then we calculate 
P e <C 1CT 100 , to be negligible. Therefore, in an NMR apparatus, if signals can be detected 
despite of the low spin polarization the choice zx can be evaluated with probability of 
almost one through the signs of the signals (|3.34|h 

Detection of signals depends on the NMR apparatus. The NMR signal intensity is 
defined by amplitude 

V s = |v^/7p^% n e;, (3.40) 



and the noise is determined through Nyquist formula, ([lilt ) 



V N = ^/4k B TRAv, (3.41) 

where Q is the quality factor of the resonance coil, V is the volume of the coil, R is 
the resistance, ujj is the Larmor frequency of the nuclear spin J, and Av is the amplifier 
bandwidth. The noise which is defined by (|3.41j) is entirely classical noise generated by 
the equilibrium fluctuations of the electric current inside an electrical conductor. Nev- 
ertheless, for an NMR experiment with two qubit liquid ensemble at room temperature 
with V = 1cm 3 and Q = 10 3 , the number of required molecules is bounded by iV > 10 16 , 

(E3). 

Therefore, in the case which we are facing, n is large to get strong enough signal 
intensities to be detectable with an NMR apparatus. Deviations are very small ()3.35|1 
and the signals are very close to the mean values ()3.34j) with the error probabilities to be 
negligible. Thus the signs of the detected signals give information on the choice of zx, 
and hence the encoded two-bit message, regardless of the separability due to the very low 
polarized initial states. Then we emphasize the results of the measurement of the spin 
magnetizations Zj and Zs to be as follows 

(Z/> = (-l) z e I} (Z s ) = (-lfe s . (3.42) 



63 



We note that this does not necessarily mean that the experiment described here is a 
demonstration of SDC with NMR mixed states, even though the signs of the detected 
signals appear to give the two-bit message while the I spin alone has been applied by 
encoding pulses. SDC relies on entanglement of states in order to realize transfer of two- 
bit information by only encoding a single qubit. In the case of the NMR experiment, 
as explained here, the large number of molecules are inevitably required for detection of 
encoded two-bit message. Then taking account the required number of molecules it is in 
principle erroneous to call the corresponding NMR experiment a demonstration of SDC. 

3.4 Appropriate Entanglement Witness for NMR 

Given the density matrix (|3.32|) . we obtain a probability for a successful SDC to be pips, 
with which a two-bit message can be transferred inside any individual single molecule. 
Therefore, for this experiment to represent the quantum information advantage and to 
outperform the classical one-bit communication through one-qubit channel, the success 
probability, pips, must exceed 50%. This is because two-bit information can be obtained 
by one-bit classical communication and a one-bit of random guess, which comes true with 
a probability of 1/2. Then SDC is beyond the classical achievements only if the inequality 

Pips > 1/2 (3.43) 

is satisfied. It is worth noting that this condition exactly coincides with the condition for 
the non-separability of p2 derived from e.g. the negativity criterion. 

Therefore, only when F, to be defined as follows, imposes a negative value, the NMR 
SDC is successful and there exists entanglement. 

F = l/2-p lPs . (3.44) 

Detection of entanglement through finding a negative value for an observable is remi- 
niscent of entanglement witness. Entanglement witness is a Hermitian operator W = 
which has positive mean values for all separable states p, Tr(Wp) > 0, but a negative 
mean value for at least one entangled state <7 ent , Tr(Wa ent ) < 0. In other words, entan- 
glement is detected if a negative mean value is obtained through the measurement of the 
entanglement witness. 

Using the spin polarizations, we rewrite F as 

F = i-i(l + ej)(l+es) (3.45) 

= ^-^(l + l<Zi>|)(l + |<Zs>|), 

where we used ()3.42|) . The absolute values are required for the evaluation of the function 
F for different choices of zx. 

Measurement on the state p% with the observables Zj and Zs is equivalent to mea- 
surement on the state p2 (or pi in the special case of x, z — 0) in the Bell basis, because 

(Zj) = Trp 3 (Z 7 ®I 5 ) (3.46) 
= Trp 2 (X I ®X s ) = (W 1 ) 
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(Z s ) = Trp 3 (I/®Z 5 ) 

= Tip 2 (Z I ®Z s ) = (W 2 ), 



(3.47) 



where the two observables W\ and W 2 are defined as follows 

Wx = U^Z, ® Is)U Be ii = X 7 ® X s , (3.48) 

^2 = 4 ell (I/ ® Z 5 )U Bc ii = Z 7 ® Z s . (3.49) 
From ()3.45|) . ()3.46|) and ()3.47|) . then F is further rewritten as 

F = /«Wi), (W 2 )) = \ - 1(1 + |(Wi>|)(l + |<W a >|). (3.50) 

Measurement of the two observables W\ and W2 on the ensemble system described by 
p 2 is related to the measurement on the state p% of the spin magnetization, which affords 
information on the spin polarizations. Therefore, separate and simultaneous measurement 
of the observables W\ and W 2 is possible in a single experiment, and tells the existence 
of entanglement. That is, if (Wi) and (W 2 ) satisfy 

F = f((W 1 ),(W 2 )) <0, (3.51) 

then the state is entangled. In this sense, W\ and W 2 are regarded as a new class of 
entanglement witnesses. These observables are easily measurable in a single run NMR 
experiment and give information on the status of entanglement through quantitative 
evaluation of the function F ([3.50)1 . 

We also note, that entanglement can be detected in principle through measuring 
the conventional entanglement witness. For p 2 , entanglement witness is derived by the 
conventional approach as 

W = -£(h®I s + (-ifX, <g> X s + (-1)*(-1) S Y 7 ® Y s + (-1)% ® Z s ). (3.52) 

If we assume the ability of implementation of any form of the unitary transformations 
with NMR pulse sequences, then the conventional entanglement witness as defined by 
()3.52|) also can be measured in a single run and by measuring the spin magnetizations (see 
appendix A). We emphasize here, that the scheme which is introduced in this contribution 
for detection of entanglement, still has advantage in the sense that for different choices of 
z and x it is applicable by just evaluation the function with the absolute values. In other 
words, there is no need to change the experimental operations. However, the conventional 
entanglement witness even though is proved here to be measurable in a single run NMR 
experiment by only measuring the spin magnetizations, still requires different pre-applied 
unitary transformations, which depend on the choices of z and x. 

3.4.1 Single Run Detection of the Conventional Entanglement 
Witness 

Theorem: Suppose that we have the observable W = WW U = W ' (aZi®Is+bIi®Zs+ 
cl (g> I)U with arbitrary coefficients a, b, c G R and a unitary transformation U G U(4), 
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and have the conventional entanglement witness W = ® Is + (— ® X s + 
{-Ifi-lfYi ®Y S + {-VfZi ® Z s ). Then there exist sets of a,b,c, and U such that 
Txp2(pi,Ps)W = Trp2(pi,ps)W for all possible pi,ps > 0.5 when we specify (z, x) to one 
of (0,0), (0,1), (1,0), (1,1). ' 



Proof: Let us write U = J^t'^o o u ki\k) I 11 this P ro °f? we make use of the following 
equations: for x = z = 0, we have 
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(3.53) 



W Q is written as 
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b + cj 



(3.54) 



In the following sentences, we restrict the problem to the case of x = z = 0. 

The eigenvalues of W are different to the ones for W. Then the two observables 
are essentially different to each other. However they can still give the same value after 
calculating the traces. 

The problem to be solved is as follows: 

Problem for x = z = 0: Suppose that we have the observable W = U^W U = 
U*(aWi + bW 2 + c)U = U^iaZi ®I S + bli <g> Z s + cl <8> I)U with coefficients a, b, c G R 
and a unitary transformation U £ U(4), and have the conventional entanglement witness 
W = \(h ®I S -X I ®X S + Y I ®Y S -Z I ® Zs). Then find a set of a, b, c, and U such 
that Tr p 2 (pi,ps)W = Tr p 2 (pi,ps)W for all possible pi,ps > 0.5. 

Solution: Firstly, consider a maximally mixed state p 2 {\, \) = I ® I /A for the input. 
We have 

Tr P2 (l,l)W = l, (3.55) 



and 



'2 7 2' 



Trp 2 (i, 1)W = Mp 2 (i l -)U^W = c. 



(3.56) 



Therefore c = 1/4 holds. We have 
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(3.57) 



where a = a + b and {3 = a — b. 
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Secondly, we change the problem to an equivalent problem of finding eigenvalues. 
Consider the unitary transformation 
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(3.58) 



With this transformation, we change the basis: 

Trp 1 (p I ,p s )W = Trp 2 (pi,p s )W 

with 

Pi(pi,Ps) = U ] HCN p 2 {pi,ps)U H cN = 



1 

2 



PiPs 
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(3.59) 



(3.60) 



We also have 

Trp 1 (p I ,p s )V^W V = Trp 2 ( Pl ,p s )W 

with V = UUhcn- We write that V = J2ki v ki\k)(l\- 

Now the problem is to find a set of a, f3 and V that satisfy 

Tr Pl ( Pl , ps^W'-V^WoV) = 0. 



(3.61) 



(3.62) 



(3.63) 



Because Pi(pi,Ps) is a diagonal matrix, the above equality hold for Wpi,ps if and only 
if the diagonal elements of (W — V^W Q V) are 0. This condition is expressed in the next 
equation: 

W - V^WoV = A, 



where 



This leads to that 
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(3.65) 
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Therefore the problem has been led to an equivalent problem: 



Equivalent problem for x = z = 0: Find an Hermitian matrix B = 



f-3 b 01 b 02 b 03 \ 

1 b 12 b X3 

1 b 2 3 

\b* 03 K 3 b* 23 1 / 



h* 
°01 

h* h* 
u 02 u 12 



whose eigenvalues are written as ±4a, ±4/3 where a, (3 e R. 

It is obvious that this problem has solutions because Tr£> = Tr diag(4a, 4/3, —4/3, —4a) 
and any unitary transformation in U(4) preserves the trace. 



To find an explicit values of a, b, and the elements of U, we solve the equation: 

B = \/ t diag(4a, 4(3, -4(3, -4a) V 
( kol 2 



^fcO%l \Vkl 



VklV* k0 v k2 v* k0 v k3 v* k0 \ 
2 v k2 v* kl v k3 v* kl 

2 



Vk0V k2 v kl v k2 \v k2 \ v k3 v k2 
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\v k0 v* k3 v kl v* k3 v k2 v* k3 \v k3 \ 2 / 
where A = 4a, A 1 = 4(3, A 2 = —4(3, A 3 = —4a. This leads to the equations: 
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Suppose that a > § and (3 — 0. Then we have the system of algebraic equations: 



4a No 2 - 4a v 30 2 = 


-3, 


(3.72) 


4a v i 2 — 4a f 3i 2 = 


1, 


(3.73) 


4a N2 2 _ 4a \v 32 2 = 


1, 


(3.74) 


4a w 03 2 - 4a \v 33 2 = 


1. 


(3.75) 



As far as we only want to show the possibility of the measurement of the conventional 
entanglement witness by only measuring the spin polarizations in a single run NMR 
experiment, it is enough to give an example and to avoid the unnecessary complexities 
due to a general proof. 

An example V ex (a) of V, as a function of a, is written as 
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Therefore an example of the set of a, b, and U as functions of a is that 
a e Ja) = b e Ja) 



Vex(a)Ulj CN . 
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Suppose that the state to be measured is a maximally entangled state. Measurement 
with the conventional witness observable W gives (W) = —1/2. We impose the observable 
W to give the same value. We observe that for the maximally entangled state, Z/ and 
Zg should take ±1. Recall that generally c = 1/4. Then, as an example, we can choose 
a ex = fr ex = 3/8 and the corresponding unitary transformation is obtained as follows 



V P J~) 



/0 



V 



1 1 r in2/3 

1 1 r -in2/3 





1 „-i?r2/3\ 



1 r iw2/3 



1 




(3.78) 



For different choices of zx, the elements of the matrix A should be changed in a way 
that the unsimilar diagonal element —3/4 stands for the matrix component A x+2z ^ x+2z . 
Then similar calculation shows that the components of U cx for different choices of zx 
interchange according to each case. 

Therefore, we showed that it is also possible to decompose the conventional entangle- 
ment witness to the separate but simultaneous measurement on the spin magnetizations. 
However, we emphasize that still the entanglement witness, which we introduced in this 
work is more handy as it covers all the cases for different z and x, by just evaluation of 
the function F with the absolute values without additional requirements on the differ- 
ent experimental operations. Whereas, as far as the conventional entanglement witness 
is concerned, different sets of a, b and more over different unitary transformations are 
required for different choices of z and x. In other words, by only changing the signs of 
the absolute values a and b but with a fixed unitary transformation, generally it is im- 
possible to get the equality for the expectation values of W and W. This can be proved 
by contradiction. 



Proposition: We already have the density matrix p x = di&g(pips,piqs, qiPsi Qils) an d 
the conventional entanglement witness family W\ = diag(— |, |, |, |), W 2 ' = diag(|, — |, |, |), 
Wz = diag(|, |, — |, g), Wi = diag(|, |, |, — |) for p\. Suppose that there is a specified- 
and-fixed unitary transformation U and coefficients a, b G R that satisfy 

TrpxWi' = Tr Pl (aU^Z ® IU + bU ] I ® ZU + -I ® /). (3.79) 

Then, it is impossible to satisfy the equality 

TrpiWV = Tr Pl (a'U*Z ® IU + b'UU ®ZU + h®I) (3.80) 

for V/c G {2, 3, 4} by using the coefficients a' = ±a and b' = ±b. 



Proof: We prove the proposition by using a proof by contradiction. 

Assume that both of Eqs. ()3.79j) and ()3.8()|) hold. Let us denote the diagonal elements 
of WZ®IU by x 00 , Xu,x 2 2, x 33 and the diagonal elements of U^I®ZU by y 00 , yn, y 22 , 2/33- 
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We have the system of equations led from Eq. (|3.79j) : 



4 JJ " JlJ 2 
Setting k = 4 in Eq. (|3.8U|) leads to the system of equations 



^ + ax 00 + by 0Q = --, (3.81) 

1,1 . 

- + ax u + fa/ii = -, (3.82) 

- + ax 22 + by 22 = -, (3.83) 

1,1 . 

+ ax 33 + by 33 = -. (3.84) 



^ + a'xoo + b'yoo = ^, (3.85) 
I + a ' Xu + b ' yn = I ; (3. 86 ) 

- + ax 22 + b'y 2 2 = ^) (3-87) 

- + a'x 33 + b'y 33 = --. (3.88) 
These two systems of equations lead to that 

(a + aOxoo + ^ + Oyoo^-^, (3.89) 

(a + a / )aru + (6 + 6')3/u = i (3.90) 

(a + a')x 22 + (b + b')y 22 = -, (3.91) 

(a + a')x 33 + (b + b')y 33 = ~. (3.92) 

The candidates of (a',b') for G {2,3,4} are (a, —6), (—a, 6), and (—a, —6). Hence each 
of the candidates must satisfy the above system of equations if the assumption is correct. 
Nevertheless, when (a', b') = (—a, —b), Eq. ()3.89j) does not hold. The assumption has led 
to a contradiction. 

We conclude as follows. The conventional entanglement witness is measurable in a 
single NMR experiment and by only measuring the spin magnetizations, provided that 
any required unitary transformation can be applied prior to the measurement. However, 
the unitary transformation has to be changed in accordance to the choices of z and x and 
in this sense the new entanglement witness, which we introduced in this contribution is 
proved to have its favorite advantage. 



3.5 Conclusion 

Liquid state NMR with very low nuclear spin polarizations prohibits the existence of 
entanglement. Then, with this physical system, it is absolutely impossible to demonstrate 
a faithful non-local quantum information processing, which requires entangled states. 
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Although two-bit information is correctly detected in NMR SDC experiment, it is not 
based on the existence of entanglement but relys on the large number of molecules. 

For a completely reliable demonstration of NMR SDC within a single molecule, spin 
polarization should inevitably be enhanced over a certain threshold and this threshold 
coincides with the condition for non-separability of the states. According to the results, 
we introduced a new class of entanglement witnesses, with several advantages particularly 
for NMR. The introduced entanglement witness is measurable in a single run experiment 
and generally is applicable for all the states without any requirement on any extra ex- 
perimental operation. Detection of entanglement through the conventional entanglement 
witness is also proved to be possible in a single NMR experiment, however for different 
states under investigation, it requires different pre-application of somehow complicated 
unitary transformations. 
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Chapter 4 

ENDOR-Based Quantum 
Computing 

4.1 Introduction 

In this thesis, the main purpose for our study has been mainly understanding the concept 
of entanglement. As a theoretical point of view entanglement has been studied in chapters 
1 and 2. Furthermore, we focused on a true realization of entanglement. Current NMR, as 
it is usually used for QC is studied in chapter 3 and has been shown to be far beyond the 
proper apparatus that can give, at least easily, quantum entanglement. Enhancement of 
the spin polarizations in NMR seem s to be in reach considering the amount of researches 
in this field, for example see and references therein. Though, we are looking for 

a rather bit different approach for realizing quantum entanglement. In order to get 
entanglement, spin polarizations should be enhanced. Electron spin, because of larger 
gyromagnetic ratio, three orders of magnitude larger than nuclear spin polarization in 
equilibrium, can be very useful in this regard. At the same time we would like to keep 
nuclear spins because of their inevitable advantages for QC that have been proved in 
NMR QC. The spin technology which satisfies our interests and involves electron and 
nuclear spins is called Electron Nuclear DOuble Resonance, ENDOR. 

Recall that NMR is a spectroscopy with high resolution, but rather low sensitivity. On 
the other hand, electron paramagnetic resonance, EPR is a highly sensitive spectroscopy 
but it does not represent much high resolution. ENDOR is a double resonance technique 
that combines the high resolution and nuclear selectivity of an NMR experiment with 
the high sensitivity of an EPR experiment. 

In ENDOR, first developed by Feher in 1956, (|93h. a change in the intensity of an EPR 
signal is observed as an RF frequency on nuclear resonance applied to the system. In 
order to get ENDOR, we first perform EPR. Then, the magnetic field is set and fixed at 
the resonance position of the EPR absorption signal. For the cw mode of the experiment, 
high power microwave is applied, giving rise to partial saturation of the EPR signal. 
For low power experiments, usually the EPR signal grows with the square root of the 
microwave power. But, if the power of microwave increases, the signal diminishes. This 
is called saturation and is important to be fulfilled to get ENDOR. Actually, ENDOR is 
desaturation of EPR by applying RF irradiation. As a result of desaturation, EPR signal 
intensity is recovered. The corresponding changes in the observed EPR signal is called 
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endor, EH), £3), US)- 

Indeed, the double resonance experiment, such as ENDOR, requires more efforts than 
simpler EPR or NMR experiments, but ENDOR comes with excellent reward on your 
further efforts that is giving you more insight into spin dynamics of electron and nuclear 
systems in the sample. ENDOR in terms of spin technology is based on a "quantum 
transformation of detection" from nuclear precession to electronic one. Particularly, pulse- 
based electron-nuclear multiple resonance is a rapidly developing field in terms of spin 
technology. 

We study ENDOR aiming on realizing an entangled state with ENDOR systems. To 
some extent, we will examine the implementation of quantum operations with ENDOR. In 
other words, full quantum computing with ENDOR has been aiming to be established. We 
get results of selective manipulations of spins, implementations of quantum operations, 
and measurements. However, quantum entanglement is just to be realized as the required 
experimental conditions are fulfilled and there is no fundamental restriction for it, as will 
be explained in the following section. 

In the following sections, some of the experiments have been done for the sake of sam- 
ple measurements for quantum computing and quantum information processing, decoher- 
ence time measurements and so on. The experimental evidences on quantum information 
processing and quantum computation with ENDOR will be given. We are mainly inter- 
ested on pulsed ENDOR for QC. However, continuous wave (cw) ENDOR and pulsed 
EPR will also be used for sample measurements and decoherence time measurements. 
For sample preparation and measurements in order to get the required information of 
the particular sample for quantum computation, we first start with cw ENDOR. Strictly 
speaking, in terms of magnetic resonance spectroscopy, pulse-based ENDOR has never 
been almighty, but cw and pulsed ENDOR technologies are complementary. 

For any physical system as a candidate for the realization of a quantum computer, 
there are some fundamental criteria, known as DiVincenzo criteria that should be met. 
Molecule-based ENDOR systems are also expected to meet these criteria in order to be 
a realistic physical system for QIP and QC; see Table l4~T1 for a list of DiVincenzo criteria 
and properties of the ENDOR system, in this regard. 

In the molecule-based ENDOR QIP and QC, molecular electron spins in addition to 
nuclear spins have been introduced as quantum bits (qubits). In a thermal equilibrium, 
the relative populations in the ground states with molecular electron spins are more 
than 10 3 times larger than the corresponding excited states in the presence of a static 
magnetic field or the ones with zero-field splitting if system involves more than two 
electrons. Therefore, with ENDOR systems, achievement of the required experimental 
conditions for preparing the initial state for QIP and QC, compared with NMR systems, 
seems to be substantially easier with the current technology. 

Any physical system for QIP and QC should be chemically stable during computa- 
tional processing. We have prepared robust organic open-shell entities against long and 
high-power irradiations of radio frequency and microwave even at ambient temperature, 
as some among those exemplified in this contribution. In addition, the corresponding 
decoherence time of the qubits is expected to be long enough compared with the compu- 
tational or operational time. As a result of the existence of the nuclear spins in molecular 
open-shell entities, there is a wide possibility to work with samples with long decoher- 
ence times. Proper samples with long decoherence time and their synthetic procedure 
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should be considered, in advance. In this context, realistic molecule-based multiple-qubit 
systems require component synthetic chemistry combined with knowledge of molecular 
magnetism. Long enough decoherence times for several samples have been measured 
during the course of this work. 

One important point which we will discuss in a separate section in the following is 
that, for a matter of initial state we still use a pseudo-pure state. For spin with low 
spin polarization, the pseudo-pure state has been pointed out to have several problems. 
It does not solve the problem of the separability of the states. Generally, it is not 
scalable because of the exponential increase in the required number of operational steps 
for realizing a pure state. However, if polarization of the spins increases to some extent, 
then the states would be entangled. Entanglement can exist for high polarized states 
although the state might be prepared as pseudo-pure states. The scalability problem for 
making pseudo-pure states is a different issue than entanglement and we might consider 
solving the scalability problem some other time specially when dealing with experiments 
with a rather larger number of qubits. For a large number of qubits, several proposals 
would be examined e .g. a lg orith mic cooling, which very recently has been shown to be 
useful for NMR QC, £3), £3)- 

ENDOR systems are essentially useful as we will show in some following section. We 
will introduce the required steps for ENDOR QC as sample preparation and measurement 
of magnetic properties, decoherence time controlling and measurements, readout. We will 
also give experimental results supporting the discussion for each section, as well as finally 
some quantum computation experiments with ENDOR will be given. 
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Figure 4.1: Energy levels and corresponding EPR/ESR and ENDOR resonance transi- 
tions in the presence of a static magnetic field. The definition for the level is as follows; 
|00) and |10) denote \M S = +1/2, Mj = +1/2) and \M S = -1/2, M/ = +1/2), respec- 
tively. Similarly, |01) and |11) denote \M S = +1/2, Mi = -1/2) and \M S = -1/2, Mj = 
—1/2), respectively. Also, see the text for the use of the notation of the sublevels. 

In QC-ENDOR, manipulation and processing on the qubits as well as the readout 
processing can be realized by introducing radio frequency pulses on the nuclear spins or 
by the microwave frequency pulses on the electron spins (see Figure 14. 1|) . 

It is known that realization of a particular class of quantum gates known as universal 
gates can be enough for implementation of any other quantum gates. One-qubit gates 
in addition to a non-trivial two-qubit gate, e.g. Controlled-NOT gate (CNOT), give a 
universal set of quantum gates. For implementations of quantum gates, it is possible to 
perform this task by introducing the relevant pulses. From the experimental side, for the 
QC-ENDOR system particular quantum gates have been demonstrated by means of the 
two qubits, as it is described in the following section. Multi-qubit gates for molecular 
systems involving a larger number of qubits should be considered for the relevant physical 
system in terms of the particular form of the spin Hamiltonian for the corresponding 
sample. There are several typical types of quantum operations based on multi-qubit- 
gates, depending upon spin Hamiltonians for realistic QC-ENDOR experiments. 

For instance, we have implemented quantum gates between a single electron and a 
single nuclear spin. For rather more general quantum operations, proper sample, electron 
spin bus molecule, has been selected. For this form of samples, even if qubits are not 
directly coupled to each other, for instance nuclei, still it is possible to perform universal 
computation. Electron spin is a degree of freedom for this class of samples, for which 
it can be coupled to any one nuclear spin, in a controlled manner. Therefore, two- 
qubit gates can be defined and easily generalized to other pairs of qubits. The original 
idea of coupling via a bus is generally applicable foe any other degree of freedom in a 
sample for quantum computing. This can be typically a photon that is coupled to several 
quantum dots embedded in a semiconductor structure. We will use this idea for rather 
more upcoming advanced experiment for a complete quantum computing with ENDOr. 
However, for preliminary experiments, also for a matter of implementation of quantum 
entanglement, quantum gates can be easily demonstrated between electron and nuclear 
spin and essentially the situation is very simple as the number of qubits are not very 
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Table 4.1: DiVincenzo criteria for ENDOR Quantum Computing. 





DiVincenzo Criteria 


ENDOR 


Qubit 


Identifiable, well character- 
ized qubits are required. 


Molecule-based electron and 
nuclear spins in molecular 
open-shell entities based on 
chemical syntheses and iden- 
tifications. 


Initialization 


Possibility to be initialized to 
a simple and fiduciary state. 


System of electron and nu- 
clear spins would be cooled 
down in a very strong mag- 
netic field, that the condi- 
tions for entanglement would 
be achieved. For prelim- 
inary experiments, pseudo- 
pure states will be used. For 
large number of qubits, the 
other approaches, e.g. algo- 
rithmic cooling, will be used. 


Decoherence time 


Long relevant decoherence 
times, much longer than the 
gate operation time is neces- 
sary. 


Long decoherence times of 
the nuclear spins and elec- 
tron spin have been available 
for demonstration of quan- 
tum operations between the 
two spins. Proper molecu- 
lar entities with long deco- 
herence times for multi-qubit 
operations are also possi- 
ble, for which stable isotope- 
labeled open-shell molecules 
have been designed and syn- 
thesized. 


Quantum operation 


Universal set of quantum 
gates is required. 


Quantum gates between a 
single electron and a sin- 
gle nuclear spin have been 
demonstrated. Electron spin 
bus molecule for ENDOR QC 
is prepared for which more 
general quantum operation 
can be applicable. 


Measurement 


The ability of measurements 
on qubits to obtain the re- 
sult of the computation is re- 
quired. 


The current measurement 
scheme is an ensemble one, 
for which modification of 
quantum algorithms are re- 
quired. This is not a problem 
as it has been already solved 
out for NMR QC. However, 
we may use also STM-based 
electron magnetic resonance 
detection. 
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large. 

In the ENDOR based QC, the readout or the measurement part has been done through 
electron spin echo detection. For detection of entanglement, as far as the global infor- 
mation on the state of electron and nuclear spins is required, modified detection scheme, 
Time Proportional Phase Incrementation, TPPI has been used. In the following sections 
these points will be explained in more details. In virtue of TPPI, we will be dealing with 
spin manipulation of electron-nuclear systems. 

It has been proved that for the case of pure states, entanglement is the necessary 
requirement for quantum exponential speed-up over the classical counterpart. For mixed 
states the statement has not been proved yet, but highly believed to be correct. Therefore, 
one very important issue which should be examined for any physical system for QIP and 
QC is the entanglement status. 

Realization of the entanglement between an electron spin and a nucle ar spin has 
been reported in an ENDOR experiment by using pseudo-pure states, f)l20h . However, 
it should be stated that the entanglement they achieved is " pseudo-entanglement". It 
means that the states actually are not entangled but the quantum operations for making 
entanglement have been applied and if the required experimental conditions i.e. low 
enough temperature below critical temperature and high magnetic field, are satisfied, 
then quantum entanglement would be established. 

In order to determine the required experimental conditions for establishment of entan- 
glement, we should consider the Hamiltonian of the system and determine density matrix 
which explains the system of interest. However, the full Hamiltonian of the system might 
be very complicated. Simplified Hamiltonian would not give all the correspondence be- 
tween entanglement and experimental conditions, but still might be interesting to some 
extent. We calculate the conditions of temperature and magnetic field in order to es- 
tablish entanglement. The measure that we use for calculation of entanglement is the 
one which we introduced previously, in chapter 3. Recall that through our discussion on 
SDC, states for two spins / and S are entangled iff pips > 1/2. Here pi and ps are the 
probabilities for nuclear and electron spins, respectively, to be at the zero state and are 
related to the spin polarizations as follows 

Pi = — g— > Ps = — — ■ I 4 - 1 ) 

Now, substituting the corresponding spin polarizations, we get the required conditions 
for entanglement as follows 

Ent. = e/e 5 + e 2 + e s > 1. (4.2) 

Here e/ and es are the spin polarizations for nuclear and electron spins. Recall that 
the above inequality for entanglement is achieved for the case that a Hadamard gate 
and a CNOT gate are applied to the state (not pseudo-pure state) of the system. How- 
ever, as far as a pseudo-pure state e.g. effective pure state can be written as a convex 
combination transformed thermal states, and any measure of entanglement is a convex 
function therefore entanglement is reduced under convex combination, then an individual 



transformed thermal state may possess more entanglement than the combination, ()12lh . 
For any other particular class of quantum operation, entanglement might be calculated 
e.g. scheme that has been explained in chapter 2. 
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Figure 4.2: Entanglement diagram as a function of temperature. Vertical axis is the 
amount of entanglement defined bv 14.21 The area in which the states are entangled is in 
shadow, a) W-band b) Q-band. c)W-band under the condition of high polarization for 
nuclear spins d) Q-band under the condition of high polarization for nuclear spins. 

Referring to the above measure of entanglement, the experimental condition for exis- 
tence of entanglement for a magnetic field of 95GHz, W-band experiment, is a temper- 
ature of 0.83 K. However, if a transfer of spin polarization from electron spin to nuclear 
spin is performed the required temperature is 5.17 K, which is well in reach with current 
technology with a W-band ENDOR spectrometer operating at liquid Helium tempera- 
ture. Figure l4~2l shows the required temperature for existence of entanglement, given by 
inequality 14.21 after substituting temperature in spin polarizations, for magnetic fields of 
95 GHz, W-band, and 35 GHz, Q-band experiments, for two cases of the experiments 
without and with transfer of spin polarization. But, how to transfer electron spin po- 
larization to nuclear spin. This is a matter that is discussed in the following section. It 
will be clear that it is an advantage of pulsed ENDOR spectroscopy that the scheme for 
transfer of electron spin polarization to nuclear spin is performed comparably easier and 
by applying relevant pulses. 



4.2 Transfer of Spin Polarization from Electron to 
Nuclear Spin 



In the case of working with pulsed ENDOR, transfer of spin polarization is easily done by 
applying appropriate pulses. Following, we explain a very simple approach with which 
by applying some particular pulses the final nuclear spin polarization becomes almost 
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as large as electron spin polarization. The approach that is explained here express the 
possibility of achieving almost equal spin polarization on electron and nuclear spins. The 
approach is shown in Figure 14 .HI and works as follows. 

Suppose an isolated system composed of an electron spin S — | and a nuclear spin 
I — |. The corresponding four energy levels are shown in Figure fOT a). In order to 
make a nearly equal spin polar ization on two electron and nuclear spins, we propose the 
following pulse sequence, (jl22T ). 

In the Figure FOj the horizontal length of each energy level refers to the corresponding 
population of that level. Nuclear spin polarization is neglected in comparison with the 
electron spin polarization. 

— = exp = 1 - — — 4.3 

= 1 - e. 

Figure l4~3"l shows the polarization of states using the above naming for e = ff . Initially, 
electron spin polarization, e$, is equal to e, while nuclear spin polarization is zero, ej = 0. 
In this Figure e is exaggerated for the sake of easy understanding. Firstly, apply an 
on- resonance microwave (^13) 7r pulse. The state after the pulse would be as shown in 
Figure l4~3T b). Now, the polarizations are interchanged as follows 

es = (4.4) 
e/ = e. 

Applying an rf tt pulse (0*12) would interchange the population between two energy levels, 
one and two, but may not affect spin polarizations. Now, we make an assumption, that 
can be satisfied by carefully choosing proper samples. Suppose that spin-lattice relaxation 
time of electron, T le , is less than the relaxation time for nuclear spin, T ln . Therefore, after 
a waiting time of T le the system may go back to equilibrium on electron level and may 
yield the energy levels as shown in Figure l4~3f d). For instance, finally we have electron 
and nuclear spin polarizations equal to ^f- and ^-(1 — e). Hence, we have 

es = e (4.5) 
ei = e. 

The above simplified approximation for calculating the population is based on a rate- 
equation method. The above spin polarizations are the values that are required in much 
higher energy levels as shown in Figure 14.31 There might be some more elaborated pulse 
sequences to get high spin polarizations on both electron and nuclear spins, without a 
serious assumption on the sample properties. The above explained approach, however 
is an experimentally easy accessible approach for the current technology of ENDOR 
spectroscopy, and as far as the pulse sequence for QC with ENDOR is not complicated 
would work properly. 



4.3 Providing Some Theoretical Concepts 

We have been mainly engaged in two experimental tasks for the realization of a quantum 
computer by molecule-based ENDOR. One has been attempts for the preparation of the 
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Figure 4.3: A scheme on how to make nearly equal spin polarizations on nuclear spin 
and electron spin, equal to the initial electron spin polarization. 
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experimental requirements for demonstrating quantum entanglement between a molecular 
electron spin and a nuclear spin by the use of a simple organic radical in the single crystal 
or the electron delocalized systems. High spin polarizations on both the electron and 
nuclear spins are essentially required to achieve entanglement between the two spins in 
the molecular frame. 

While the preparation of all the experimental requirements is in progress, the efforts 
still can be maintained on some other aspects of the research. The other aspects are 
materials challenge to design and synthesize stable open-shell molecular entities suitable 
for QIP/QC ENDOR experiments. Sample selection and measurements are very heavy 
experimental tasks. Novel molecular open-shell systems with stable isotope labels suitable 
for our purposes have also been designed and synthesized. The critical temperature can 
be tuned by invoking stable high-spin molecular entities in near the future. 

ENDOR consists of EPR, electron paramagnetic resonance, followed by NMR. EPR is 
a spectroscopic technique, which detects species that have unpaired electrons. Materials, 
which have unpaired electrons are various, including free radicals, many transition metal 
ions, and defects in materials. 

Here we would like to first explain theoretical fundamentals of EPR/ENDOR, (jl23l ). 



(124). It should be noted that almost the same explanation goes for NMR. However, 
differences in the magnitudes and signs of the magnetic interaction involved for NMR 
and EPR/ENDOR, naturally lead to divergences in the experimental techniques that 
have been employed. I would like to organize a comprehensive discussion, including 
NMR too, in this section, but because of the space and the main interest topic for this 
work, I would better keep going with the most important points from EPR/ENDOR, that 
are required in order to understand the following experiments on ENDOR QC including 
material search and so on. 

Electron with S — 1/2 has a magnetic moment /i e , which is proportional to the 
magnitude of the spin, that is 

/i e = -gPS. (4.6) 

Here, g is a dimensionless constant, generally a tensor, called the electron g factor. 
(3 is the electronic Bohr magneton defined as (3 = with — e and m representing, 
respectively, the charge and mass of the electron. hS is the spin angular momentum 
vector. 

If we apply a steady magnetic field H, there is an interaction between the field and 
the magnetic moment fi e as follows 

H = -n e ■ H. (4.7) 
If the applied static magnetic field is in the z direction then 

H = g(3HS z . (4.8) 

As far as S = 1/2, there are two allowed orientations of spin, parallel or antiparallel 
to H, and are represented by \(3) corresponding to m s = —1/2 and \a) corresponding to 
m s = 1/2, Figure IP 

In order to introduce transitions between the two nuclear spin levels, an oscillating 
electromagnetic field is applied to the system. Absorption of energy occurs providing the 
magnetic vector of the oscillating field is perpendicular to the steady magnetic field H 
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Figure 4.4: EPR energy levels. 



Table 4.2: Different EPR spectrometers with corresponding frequencies. 





Frequency/ 
GHz 


Magnetic 

Field/ 

G(g=2) 


L 


1-2 


500 


S 


2-4 


1000 


X 


9-10 


3400 


K 


24 


8500 


Q 


34 


12100 


w 


94 


33500 



and also provided the frequency v of the oscillating field satisfies the resonance conditions 

hu = g(3H. (4.9) 

The resonance frequencies for different magnetic fields are listed in Table l4~2l 
It should be stated that the oscillating field is equally likely to produce transitions 
up or down. However, due to the fact that the population of (3 electron is slightly larger 
than a in a thermal equilibrium, then we can only detect absorption of energy from the 
radiation field. 

There are two possibilities. Conventionally, we could make the magnetic field, B 
constant and vary the frequency, microwave until the signal is detected on the condition 
of resonance. However, because of the restriction on the Microwave technology, this 
option is not practical. Instead, we apply a constant microwave frequency and vary the 
magnetic field close to the resonance condition. 

Suppose the magnetic energy level around 3000G. Then, in a thermal equilibrium, 
from the Boltzmann distribution, we have 



N, 



up 



N, 



exp( 



down 



kT ' 



(4.10) 



There would be a net alignment of the spins, however very tiny. For example, in 
a magnetic field of 9.8 GHz, X-band, and at room temperature, with a sample which 
contains 10, 000 spins, on average 5, 004 spins would be parallel and 4996 spins would be 
antiparallel. The resulting difference population would be 8, giving the spin polarization 
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almost 1CT 3 . Then, levels are almost but not exactly, equally populated. If a microwave 
magnetic field is applied perpendicular to the static magnetic field, some transitions 
are acquired between the magnetic levels. Microwave stimulates transitions in both the 
directions. Since the population in the lower level is slightly larger than the upper level 
then there will be a net observed EPR signal. Under the steady condition the process 
would be soon stopped as the population would be equalized. The system is subject to 
interactions between spins of the surrounding described by characteristic T\ and T 2 , spin 
lattice relaxation process and spin-spin relaxation process, respectively. These make the 
longitudinal and transverse transitions. 

After saturation of EPR, ENDOR is performed by monitoring the EPR signal inten- 
sity, while sweeping an RF signal to drive the NMR transitions in the cw conventional 
mode of experiments. Shortly speaking, there are some main steps to acquire an ENDOR 
spectrum. 

1. Acquire an spectrum and determine the magnetic filed values at which the EPR signals 
occur. 



2. Set the magnetic field to the resonance position of the absorption EPR signals. 

3. Increase the microwave power and partially saturate the EPR signal. 

4. Monitor change in the EPR signal, while sweeping the RF. For instance, it should 

be considered as the signals are too weak and need RF amplifier to make them 
detectable. 



4.3.1 The Magnetic Hamiltonian, EPR/ENDOR Spectra 

In this part, we consider a simple system of one proton and one electron spin. This can 
be illustrated by Hydrogen atom. We first derive expressions for the energy levels of this 
atom in a magnetic field. Then, possible spin transitions and EPR/ENDOR spectra will 
be calculated. 

Let us give an overview of what will come in following. First of all an appropriate wave 
function should be defined that can describe the four energy levels from one electron and 
one proton. These energy levels are degenerate in the case that there is no magnetic field. 
Then applying magnetic field H makes the system to involve the interaction between H 
and S in addition to H and I. We would also study coupling between spins, hyperfine 
coupling. This particular simple example possesses most of the complications that we 
would require to understand and make it clear how to read EPR/ENDOR spectrum. We 
should be careful however that the simplification due to this particular example might 
be dangerous. Specially because that the hydrogen atom is spherical, then the spin 
Hamiltonian for the hydrogen is not anisotropic, though we need to study the anisotropical 
behaviour of g and A tensors. Therefore, l ater w e will e xten d this discussion in a way 



that can include the anisotropic terms, too, (|125l ). (jl26l ). ([123T ). (J,124). 
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The Zeeman Energy 

The nucleus and the electron both interact with the steady magnetic field and yield the 
following energy 

H = gf3HS z -g N p N HI z (4.11) 

g and are the g factors for the electron and the nucleus. The second term, the nuclear 
energy is much smaller than the electronic, first term, part. 

The Isotropic Hyperfine Coupling 

The magnetic moments of electron and nuclei are coupled via the so-called contact inter- 
action, which represents the energy of nuclear moment in the magnetic field produced at 
the nucleus by electric currents that are associated with the spinning electron, it has the 
form 

Hx = al ■ S = a(I x S x + I y S y + I Z S Z ). (4.12) 

There is a magnetic coupling between the magnetic moments of the electron and 
nucleus which is entirely analogous to the classical dipolar coupling between two bar 
magnets. But the dipolar Hamiltonian averaged out to zero whenever the electron cloud 
is spherical. 

Totally the complete Hamiltonian is 

H = g(3HS z - g N (3 N HI z + aS • I. (4.13) 

It is the task of the experimentalists to find the proper Hamiltonian with parame- 
ters. To make it more general we should consider that Hamiltonian includes terms as H-S, 
H • I and S • I. The first two come in the 7io and the last term in TrL\. In order to solve the 
corresponding hamiltonian, we would require perturbation theory. The eigenfunctions of 
Ti.o, 0i, •••,04, correspond to the "zero-order" unperturbed energy values, ei,...,64. Here 
7ii is the perturbation term that yield 

. (m\Hi\n) . . 

i>n = <t>n - 2j _ 0m (4.14) 

m^=n 

E n = e n + {n\Hi\n) - 

the last equation for energy includes the first-order correction in its second term and the 
second-order correction in its last term, it is clear that in the first order of perturbation 
the energy of 0„ is changed without changing the form of n , however in the second-order 
of perturbation the changes of energy also mixes n with other function n . This will be 
clear in the experiment with DPNO that we will show in the following sections. 

Zero-Order Energy 

Recall that we denoted the two possible spin function of an electron S = 1/2, with \a e ) 
and \P e ) corresponding to m s = +1/2 and m s = —1/2, respectively. These come from 
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the eigenequations as follows 

S z \a e ) = +l/2|a e >, (4.15) 
S z \(3 e ) = -1/2|&>- 
Similar fact is going on for nuclear spin / = 1/2, as follows 

I z \a N ) = +1/2K), (4.16) 
Wn) = -1/2\{3 N ). 

Therefore, at this level, the appropriate basis functions are 

0i = \a e a N ), 02 = \a e p N ), (4.17) 

03 = |/5 e ajv), 04 = \PepN)- 

Initially the four states are degenerate, but after applying the steady magnetic field 
represented by Hamiltonian 7i then the states will separate as follows, e.g. 

H \a e (3 N ) = gf3HS z \a e )\f3 N ) - g N (3 N H\a e )I z \[3 N ) (4.18) 

= ^9PH\a e ) \p N ) + -g N (3 N H\a e ) \(3 N ) 

= (^gPH+^g N H)\a e p N ). 

Then the zero-order energy values are 

ei = ^g[3H -^g N (3 N H \a e a N ), (4.19) 

£2 = ^g(3H +^g N (3 N H \a e /3 N ), 

£3 = ~^gpH-^g N p N H \/3 e a N ), 

e 4 = ~g/3H + -g N /3 N H \(3 e (3 N )- 

The First-Order Energies 

At this level Tii works as a perturbation on the zero-order and it gives results as follows 

Ti ± = aS ■ I = a(S z I z + SJ X + S y Iy), (4.20) 

for the first-order, the first term is only involved and the latter terms are for the second- 
order and will be discussed later, hence, we have 

(a e a^\ciS z I z \a e aiy) = —a, (4-21) 

{(3 e a N \aS z I z \a e (3 N ) = --a, 

(f3 e a N \aS z I z \(3 e a N ) = ~a, 

{(3 e (3 N \aSJ z \(3 e (3 N ) = -a. 
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The Second-Order Energies 

This is due to the term a(S x I x + S y I y ). We recall that shift operators are defined as 
follows 



S + — S x + iSy, 

S = S X tSy. 

Therefore, considering also the relation for shift operators, we have 

sj x + s y i y = ^(s + r + s-i + ). 



The nonvanishing matrix elements are 



{a e [3 N \S + I \/3 e a N ) 
(f3 e a N \S~I + \a e f3 N ) 



1, 
1. 



(4.22) 



(4.23) 



(4.24) 



The complete matrix of the operator Hi is 



al • S = -a 

4 



/l 0\ 
0-120 

2-10 

\0 1/ 



Then perturbation theory changes the states into new states as follows 

a 



i>2 = \a e pN) + 

Corresponding to the energies 

E2 
E 3 



2(gf3H + g N p N H) 
a 

2(g(3H + g N (3 N H) 



\(5 e a N ), 



\a e (3_ 



N 



11 1 

(- 9 W + - gN p N H) - - + i{gl3H l (M3NH y 

11 1 ^ 

-(- 2 9PH + -g N p N H) - - - H!lUI " !lxjxll) 



(4.25) 



(4.26) 



(4.27) 



We just remind that for the simple case as the hydrogen atom, we would solve the problem 
even without using perturbation theory. 



The First-Order Spectrum 

Suppose that an oscillating magnetic field of strength 2Hi cos cut acts on the atom. For 
transitions applying a proper e.g. microwave frequency may cause allowed transitions 
like a e aN — > Pe^N, or forbidden transition, a e ^N — ► fieOiN, with very small probability. 
On the other hand, RF pulse causes allowed transition a e a^ — > a e f3^. 
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If the oscillating magnetic field is applied along the z direction, it will not cause any 
energy absorption, because it just modulate the energy levels of the spin system. On the 
other hand, if it is applied across the steady field H in the x direction we have 

V(t) = 2(gPH 1 S x -g N p N H 1 I x )cosu;t (4.28) 
= 2V cos u>t. 

The transition probability from state n to the state m is equal to 

P nm = ^2 | (n\ V\m) \ 2 5(uj mn - u). (4.29) 

If Hi corresponds to the electron spin then resonance transitions of electron spins are 
aquired as follows 

P nm = ^g 2 P 2 H*\(n\S x \m)\ 2 8(u mn -u). (4.30) 
If we use the shift operators, then we have 

(a e a N \S x \(3 e a N ) = ^(a e \S + \f3 e ) = ^, (4.31) 

with the transition probability 

P = £pg 2 l3 2 H*6(u mn -u). (4.32) 

We replace 7 = g/3. Also, to be closer to the real situation, we write the above 
equation as based on the general line shape function g(u>). 

P= l -ix 1 2 Hlg{uo). (4.33) 

And the frequency for the transition would be 

hu a = (± g pH-±g N p N H+±a)-(-±gPH-±g N p N H-±a) (4.34) 
= g(3H + \a. 



hv h = (^g(3H+^g N p N H + ^a)-(-^gPH+^g N p N H+^a) (4.35) 
= gPH-±a. 

Then the first order EPR spectrum consists of two lines separated by a, the hyperfine 
constant. Now, the magnetic field is set to the top for the EPR. Then microwave power is 
increased to saturate the EPR signal. ENDOR comes by monitoring EPR signals, while 
sweeping the RF irradiation. For instance it should be considered as the signals are too 
weak and need RF amplifier to make them detectable. Suppose that the microwave is 
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fixed on the lower state with energy —l/2g/3H and the RF is applied. The resonance 
frequency would be given in a same manner as above but this time for ENDOR as follows 

h/i a = {-^g(3H + -g N (3 N H + -a) - (-^g/3H - -g N /3 N H - -a) (4.36) 
= g N f3 N H + - 



hy. h = {-\g(3H +\g N l3 N H - -a) - {-\g(3H -\g N (3 N H + ~a) (4.37) 
= 9nPnH - - 

The Second-Order Forbidden Transitions 

The transition in which both electron and nuclear spin states are changed is strictly 
forbidden in first-order, but in the second-order, I4.2(j[ it is weakly allowed, provided the 
oscillating magnetic field is polarized parallel to the static field, rather than perpendicular. 
In the following calculation we used A = n , au , a a m . We used the states which we 

& 2(gf3H+g N l3 N H) 

already calculated for the first order perturbation, then we have 

(il>2\S z \i/;3) = (a e (3 N + \(3 e a N \S z \p e a N - \a e (3 N ) (4.38) 
= (a e /3 N \S z \(3 e a N ) + \(P e a N \S z \p e a N ) 
- \(a e p N \S z \a e p N ) - \ 2 ( 2 (f3 e a N \S z \a e p N ) 

= 0--A--A-0 
2 2 

= -A. 

A is small, then A 2 is neglected. 

The corresponding transition probability is 

P = 2\ 2 -K 1 2 Hlg(uj). (4.39) 

We would like to close this part here and move to more precise explanation particularly 
for liquid and solid state sample as will be directly used for sample study for ENDOR QC. 
But first we would like to note that there are much more corrections which we neglected 
here. Also, generally g and a are anisotropic and represented by tensors. This will be 
clear by studying the orientational speciment (single crystal), the principal values of the 
g tensor and/or a tensors can only be determined from solid state study. 



4.3.2 Liquid State Sample Study 

Here, we concentrate our attention on liquid state sample study. We shall still use the 
general discussion that is given above and only add some more important points. We will 
see that incase the number of nuclei is large, ENDOR gives much more clear spectra as 
compared to EPR. 

Electron nuclear hyperfine interaction, in general, is represented as follows S • A • .1, 
where A is the hyperfine tensor. This tensor is decomposed into the sum of two terms 
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aS ■ I + S • A ■ I, where a is the isotropic constant part, and A is the magnetic dipolar 
tensor. In solid state it gives the anisotropic term that shall be discussed in the following 
section. However, the average value of A is zero and the observed hyperfine splitting gives 
the isotropic coupling constant directly. 

The spin Hamiltonian for a system composed of an electron spin and one nuclear spin 
can be written as follows, in a very simple form. We have assu med th at q and A are 



isotropic scalars as would be the case for a liquid phase sample, (I114T ). ()128l ). Also, we 
have the approximation of the first order perturbation theory or high field limit. In other 
words, the strong field is taken in the z direction and mj and mj are valid quantum 
numbers. This approximation ignores second order hyperfine terms, h 2 A 2 /g(3H. 

U = g(3HS z + SAI - g N (3 N HI z . (4.40) 

The first term is electron Zeeman effect and actually it is the most dominant term in 
the Hamiltonin. This interaction determines the magnetic field at which the EPR signal 
is centered. The information on the electronic structure of the radical can be extracted 
through this term. However, some important information such as the molecular structure, 
spin distribution, or dynamics of the sample would not be possible to be measured by this 
term. Nuclear Zeemann interaction is represented with the last term. It gives information 
on the structure, or identity, of the nuclear spin. This term is comparably small. Usually, 
it is not possible to be measured in a typical EPR experiment because of the selection 
rule and more adequate ENDOR experiment would make it detectable. 

The second term in the above Hamiltonian is due to the interaction between electron 
and nuclear spins, hyperfine interaction. Very essential information on the structure, 
identity and number of atoms involved in the radical, their distances from the unpaired 
electron and electron spin densities are available through the hyperfine couplings if some 
quantum chemical considerations are made. The other term, not appearing in the above 
equation, is the electric quadropolar interaction for the nuclear spin. P is the quadropolar 
matrix, generally a tensor and vanishing in solution. The nucleus with a spin larger than 
1/2 gives a non-zero quadropolar term. 

We define the electron and nuclear spin resonances as follows 

v e = g(3H/h (4.41) 
v n = gw^wH/ h 

v e and v n are the electron spin resonance and nuclear spin resonance, respectively. Then, 
energy eigenstates are calculated as follows, see Figs l4.5l and 14.61 

Em s , mi /h = v e m s - v n mi + am s mi. (4.42) 

Selection rules give the allowed transitions for EPR and ENDOR as follows 

Am s = ±l,Am ; = 0, (4.43) 

for EPR, and 

Am s = 0, Am/ = ±1, (4.44) 

for NMR. 



90 




Figure 4.5: Energy levels for a system composed of one electron spin m s = 1/2 and a 
nuclear spin mj = 1/2. Here, v n > \A\/2, and A < and g n > 0. 
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Figure 4.6: Energy levels for a system composed of one electron spin m s = 1/2 and a 
nuclear spin mj = 1/2. Here, \A\ > 2u n , and A < and g n > 0. 
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Figure 4.7: A simple illustration of ENDOR Spectra. 



The typical EPR spectrum for this system, composed of an electron and a nuclear 
spin 1/2 is two lines being centered at v e and separated by see Figure |4~T1 For NMR, 
there are two possibilities. If v n > \A\/2 then NMR signals are two picks, centered at v n 
and separated by \ A\. If u n < \A\/2 then the two lines are centered at \A\/2 and separated 
by 2v n . This is a direct explanation of energy eigenstates, as given by the above equation. 
This is shown in Figure 14.71 and the energy levels are shown in Figs. I4.5H4.6I Through 
the calculation, parameters for any molecule being as the sample can be measured. We 
show in the following section a particular experiment, which is done in liquid phase. 

In order to show the difference between EPR and ENDOR for the case that the number 
of nuclei is large let us assume, for simplicity, the situation in which a free radical involves 
two nuclei. The Hamiltonian of the system could be 

H = g(3HS z - g N (3 N H(I lz + I 2z ) + aiS.Ii + a 2 S.I 2 . (4.45) 

For EPR the terms describing the nuclear Zeeman interactions can be omited be- 
cause they do not affect EPR. However, they do affect ENDOR. Up to the first order 
perturbation, for EPR, we have 

H = g(3HS z + S z (aJ lz + a 2 I 2z ). (4.46) 

The energies are proportional to 

E = m s [g(3H + aimi + a2rri2\. (4.47) 

Usually, the frequency is fixed and field is swiped in a small range on either side of 
the field that corresponds to the condition on resonance for electron H*. We may convert 
the coupling constants a x and a 2 into Gauss. (1 Gauss= 2.80 MHz, for g of free electron.) 
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Then, we have 



H = H* - aimi - a 2 m 2 (4.48) 
= H* — aifrii. 

i 

In EPR spectrum, for two nuclear spin sample, there are four lines. However, it often 
happens that several nuclei are equivalent to each other, usually occupying symmetrically 
equivalent positions in the molecule. Then, particularly in the above example we have 
only three lines, and then relative intensities are 1:2:1, different to the previously 
equivalent intensity. Generally speaking, for one electron and n equivalent e.g. protons, 
the EPR would show n + 1 lines whose relative intensities are proportional to the coeffi- 
cients of binomial expansion of (1 + x) n . For more than one electron also there would be 
some corresponding overlaping. For instance, naphtalene negative ion contains two sets 
of four identical protons. One set would give EPR with five lines with relative intensities 
as 1 : 4 : 6 : 4 : 1. The other set would separate lines. Hence, complete spectrum contains 
25 lines and its analysis yields the two hyperfine couplings. By the intensity pattern the 
number of equivalent nuclei would be known. 

However, ENDOR, as a double resonance technique, offers superior resolution to the 
hyperfine coupling constants compared to EPR. The nuclear g values, gN is also measured 
in ENDOR, fascillating the identification of nucleous. The only information lost because 
of the fewer number of signals is the number of equivalent nuclei. But this can be achieved 
in advance through EPR. 

The dramatic resolution enhancement achieved by ENDOR results to a large extent 
from the fact that two resonance conditions have to be fulfilled simultaneously, as to be 
clear from the Hamiltonian of the system. Then, it can highly simplify the spectrum. 

The Unpaired Spin Density 

The unpaired electron density or spin density p{v N ) at the nucleus or at some other point 
in space is a probability density, measured in electrons/ (angestroms) 3 . In a real situation, 
a molecule contains many electrons, all of whose spins are couples together, and it is not 
in general possible to say there is just one "unpaired electron" in a certain orbital, while 
the other electrons are perfectly paired to each other. 

The correct contact Hamiltonian for molecule is represented by the operator as follows 

H c = y9P9nPn 6 ^ ~ r N)S k .I, (4.49) 

k 

where the sum runs over all the electrons, and the value of the splitting constant is found 
by averaging over the full many-electron wave function. The result can be written as 
follows 

47T 

a = —gf3g N f3 N p{r N ), (4.50) 
where p(r N ), the unpaired electron density at the nucleus is defined as 

p(r N ) = J 0* 2S ^(r k ~ r N )^dr, (4.51) 
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and the spin component of the state ip is taken to be S z = 1/2. p{tn) is simply the 
difference between the average number of electrons at the nucleus which have spin a 
and the number with spin /3; the operator 2S z k gives a factor ±1 in the above equation 
depending on the spin, while the 5 function ensured that the electron is at the nucleus. 
If the nucleus is in a position where an excess of electrons have j3 spin the value of p(r/v) 
will be negative, and one may speak of a negative spin density. 

On the other hand, the unpaired electron, in a molecular radical is often distributed 
over many atoms and gives proportionally smaller splitting. Thus, the probability that 
the electron is on the particular nuclear spin is quite small. Then, the spin density in an 
orbital, pu, is a number, representing the functional population of unpaired electron, on 
an atom. There is a direct proportionality between the spin density and the hyperfine 
splitting. For instance, McConnel's relation is a good, not perfect, relation to provide a 
means of interpreting the spectra of virtually all aromatic free radicals. We will see some 
relations on this regard in the section for experiment with liquid state sample. 



4.3.3 Solid State Sample Study 

Many unstable radicals can be trapped in solids, and give highly interesting results. The 
radicals may be formed by ultraviolet irradiation, electron bombardment, or exposure 
to nuclear radiations. In this work we consider single crystals. The radicals are almost 
always regularly oriented relative to the crystal axes, and this provides the opportunity to 
rotate the crystal in the magnetic field and measure the anisotropic magnetic interactions. 
Although the presence of anisotropic tensors makes the analysis and interpretation of 
spectra more difficult than in solution, the results give valuable information about the 
structure and orientation of the unstable molecules. 

Here we would like to have an overview on how to read spectra in a si ngle c rysta l 



and to find out the anisotropic g and hyperfine tensors. We have used (|123[ ) and ()124l ). 
however other textbooks might be useful too. 

In a particular crystal sample study, generally there would be no prior information 
on the principle axes of hyperfine and g tensors and these would be determined through 
the experiment. The elements of g and A might be evaluated relative to an orthogonal 
system that can be identified by obvious structural features of the crystal. Then, through 
diagonalization of the g and A tensors, principle axes would be determined. Therefore, 
some considerations on the crystal systems and corresponding angles would be helpful, 
and might be found in some related text books. 



Evaluation of the g Tensor 

In an experiment as you may see in the following section, g tensor evaluation for in- 
terested sample, DPNO, has been done through EPR measurement. To the first order, 
one can obtain g-only experimental curves by taking measurements at the center of the 
symmetrical hyperfine multiplets. 

Second order effects can, however, skew the center of the hyperfine patterns, when 
the hyperfine splitting is large. In case the first order is just contributed the analyses 
would be as follows. Measurement of the sample would be done in different planes. For 
instance, suppose measurements are taken in the xy plane where H forms an angle 9 with 
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the x axis, the effective g 2 B for various 9 values is given as follows 

9e = (g xx cos9 + g yx sin9) 2 + (g xy cos9 + g yy sm9) 2 (4.52) 
+ (g xz cos9 + g yz sin9) 2 

= g 2 x cos 2 9 + g 2 sin 2 9 + 2(g xx g yx + g xy g yy + g xy g yz ) sin 9 cos 9, 

where 

gl = glx + gl y + glzi ( 4 - 53 ) 

g 2 y = g 2 yy + gy X + gy Z , (4-54) 

and where g x and g,, are the effective g values along the x and y axes. The above equation 
can be rewritten as the following form which is more convenient for plotting. 



g 2 = K^x, y) + K 2 {x, y) cos 29 + K 3 (x, y) sin 29, (4.55) 



where 



K 1 (x,y) = ^(g 2 x + g 2 y ), (4.56) 



K 2 (x,y) = ^(g 2 x -g 2 ), (4.57) 

K 3 (x, y) = g xx g yx + g xy gyy + 9xz9 y z- (4.58) 

Expressions for the other coordinate planes can be obtained by obvious exchange of 
the coordinates. 

Measurements along the three axes give the three effective values g x , g y , and g z , 
and measurements in the three planes away from the axes at any 9 can give only three 
additional constants. Because the tensor is symmetric, g xy = g yx , and so on, there are only 
six unknown tensor elements to be evaluated, and there are six undependent equations 
linking these unknowns to the measurable constants. Thus all elements of the g tensor 
are determinable. 

The effective gg values to be used in the 14.551 may be obtained from the relation 

ge = T^r, (4.59) 

where Hg are the resonance field values at the fixed frequency uq for the different orien- 
tations 9. 

In case that the matrix elements of g expressed in the arbitrary system x, y, z, as 
described in the preceding paragraphs, the resulting tensor can be diagonalized by a ro- 
tational transformation of the tensor to the principal axes. Let k x , k y , and U z be the 
direction consines between the principal axis i and the refrence axes x, y, and z, respec- 
tively, and let gi be the corresponding principal value. Since the off-diagonal elements 
vanish in the principal system, the transformation may be expressed as 

g X x gxy gx Z 

g y x g yy g yz I \ h y I = gi \ k y I • (4.60) 
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This is the problem of finding the eigenvalues and eigenvectors. Solution for this 
problem gives three values of g^ which are the principal values of g. We use the designa- 
tions g u , g v , and g w , where g u is the smallest, g v is the intermediate and g w is the largest 
principal element of g and where u, v, and w signify the principal axes. Substitution of 
each of the principal values in turn for gi yields a set of three equations involving the 
three direction cosines for the particular principal axis. These equations with the auxil- 
iary equation lf x + lf y + lf z = 1 may be solved for the direction cosines for the particular 
principal value i = u,v,w that is substituted. 

As a method of procedure, it is usually simpler to obtain the elements {g 2 )ij of the 
squared g matrix from 14. 5B1 than those of the unsquared g. One can then diagonalize the 
squared matrix and take the square root of the diagonal elements to obtain the principal 
g values. Since the g values are always positive, this procedure causes no ambiguity in 
sign. The relationship of the squared to the unsquared elements is easily found from the 
matrix product rule: 



(g 2 )ij = (i\g 2 \j) = (^M^K^b) = 

k=x,y,z k 



(4.61) 



The resulting expressions are simplified because the g matrix is symmetric, that is, 
9ij = 9ji- F° r example, 



(g )xx gxxgxx + gxygyx ~t~ gxzgzx g X x 9xy 9, 



9x 



Similarly, (g 2 

[9% 



gl {g 2 



9l, and 



gxxgxy ~\~ gxygyy ~\~ gxzgzy (g 



lyx 



K 3 (x,y) = K 3 (y,x). 



(4.62) 



(4.63) 



With the corresponding relationship for the other coordinates, g square is expressed 
in terms of the observed constants by 



(g 2 )xx (g 2 )xy (g 2 )xz 

{g 2 )y* {g 2 )yy (g 2 )yz 

{g 2 )zx {g 2 ) zy {g 2 )zz 



gl 



K 3 (x,y) K 3 (x,z) 
K 3 (y,x) g 2 K 3 (y,z) 
K 3 (z,x) K 3 (z,y) g 2 z 



(4.64) 



where the g 2 , g 2 , and g 2 z values along the respective axes are related to the K constants 
by 

g 2 x = K^y) + K 2 (x,y) = K L (z,x) - K 2 (z,x), (4.65) 



gl = K i{v, x ) + K 2{y,x) = K x {x,y) - K 2 (x,y), 



(4.66) 



g\ = K^z, y) + K 2 (z, y) = K^y, z) - K 2 (y, z). 



(4.67) 



Then solving the eigenvalue problem gives the three values of gf that are the squared 
principal elements, g 2 , g 2 and g^, of the g tensor. 

The direction cosines, or eigenvectors, for the principal elements are, of course, the 
same as those for the squared principal elements and hence may be obtained from solution 
of the set of equations 



[K x {x, y) + K 2 (x, y) - g 2 ]l ix + K 3 (x, y)l iy + K 3 (x, z)l 



0. 



(4.68) 



97 



K 3 (y,x)l ix + [Ki(y,x) + K 2 (y,x) - gf\l iy + K 3 (y,z)l iz = 0, 



(4.69) 



K 3 (z,x)l ix + K 3 {z,y)l iy + [Kfax) + K 2 (z,y) - g\\l iz = 0, (4.70) 

where i = u,v,w and where the values g 2 ,, g 2 , and g 2 obtained from eigenvalue problem 
as above, are in turn substituted for g\. The solution of these equations yields only the 
ratios of the direction cosines; the addition ^+^+^ 2 = 1 is required to give their specific 
value. For instance, the analysis of the DPNO, we used the both ways as described above. 

Evaluation of Hyperfine Coupling Tensor 

In a similar manner to that explained for the evaluation of the g tensor as explained 
above, here also we would scape the theoretical background that might be found in text 
books and will show how to evaluate hyperfine tensor from the detected spectrum. We 
neglect the effect of g tensor anisotropy. This assumption is valid for most free radicals. 
For these cases, g tensor anisotropy is sufficiently small that its effect on the nuclear 
coupling tensor is negligible. Then, almost the same process as for g tensor would be 
performed for A. 

We start calculation for the effective coupling in the xy reference plane. 

A 2 g = A\ cos 2 + Ay sin 2 9 + 2(A xx A yx + A xy A yy + A xz A yz ) sin 9 cos 9, (4.71) 

where A x and A y are the effective coupling along the x and y axes, which are expressed 
in terms of the squared tensor elements by 

A\ = A\ x + A 2 xy + A 2 XZ , (4.72) 

■^v = ^yy + + tfiz- (4.73) 
Note that these equations have the same form as those for the g tensor and hence 
may be expressed as 

A] (x, y) = d (x, y) + C 2 (x, y) cos 29 + C 3 (x, y) sin 29, (4.74) 

where 9 is the angle between the x axis and the applied field, and where 

Ci(x,y) = l(Al + Al), (4.75) 

C2(x,y)= 1 -(Al-Al), (4.76) 

C 3 (x, y) = A xx A yx + A xy A yy + A xz A yz . (4.77) 

Similar expressions for the other coordinate planes may be obtained by obvious 
changes of the coordinates. 

By measuring the hyperfine splittings in the three coordinate planes and by squaring 
the resulting values, one can evaluate the six independent C constants. The elements of 
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the squared hyperfine coupling may be obtained from these measurable constants by use 
of the relations 

(A 2 ) xx = A 2 X = C x (x,y) + C 2 (x,y) = C x {z,x) - C 2 (z,x), (4.78) 





= K 


= Ci(y,x) + C 2 (y,x) = 


Ci(x,y) - 


C 2 (x,y), 


(4 
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which are similar to those for squared g values. One substitutes the numerical values of 
squared hyperfine coupling and solves for the three eigenvalues A 2 = A^, A 2 ,, A 2 ^, which 
are the squared principal elements. The direction cosines for each of the principal axes 
i = u,v,w would be the eigenvectors. The direction cosines of the principal axes of the 
nuclear tensor are often the same as those of the g tensor, but not generally. 



4.4 Practice of EPR/ENDOR 



It might be easier to understand EPR/ENDOR, if we have a look at the components 
for an EPR/ENDOR spectrometer. The heart of the EPR/ENDOR is the resonant 
cavity, which contains the sample and directs the microwave to and from the sample. 
The absorption of microwave by the sample is monitored, amplified and recorded by the 
detection with modulation systems. To provide a resonant condition requires a stable, 
linearly variable, and homogeneous magnetic field of arbitrary magnitude. 

The main components of an ENDOR spectrometer is shown in Figure 14.81 The first 
significantly visible part is a magnet, which turns the electric energy levels. Then, the 
important part is a microwave resonant cavity, an ENDOR RF coil, which are metal box 
in order to amplify weak signals from the sample. Sample is placed in the cavity. In our 
cw experiment, a TM mode cavity has been used. The important part is bridge that 
includes the electromagnetic radiation source, microwave, and the detector. In addition, 
similar to other s pectr os copi es, con sole contains signal processing and control electronics 



and a computer, (jll4 ). (|l2ffl . (|l27fl 



Microwave bridge contains the microwave source and the detector. The microwave 
radiation is originated by a vacuum tube called a Klystron. In our experiment a GUN os- 
cillator has been used. The reason why we fix the frequency of microwave and change the 
magnetic field is that microwave resources can be tuned only over a very small frequency 
range comparing with magnitudes of interactions appearing in the spin Hamiltonian. 

The microwave is app lie d to the sample inside the cavity through a small hole in the 



cavity called an iris, (JUJ), ()126l ). Then iris should be properly adjusted in a way that the 
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Figure 4.8: Block diagram of (pulsed QC-)ENDOR setup, X-band (10 GHz). 
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microwave power changes at the detector only because of the absorption of microwave 
energy by the sample. In order to adjust the level of microwave power, an attenuator 
is placed in between the output of the microwave source and the sample. Attenuator 
controls the microwave power by blocking the flow of microwave radiation. There is also 
a circulator inside the bridge, which determines the trajectories for which each input 
and output microwave should follow. In order to detect the signal, there is a Schottkey 
barrier diode, which converts the microwave power to an electric current. For an accurate 
qualitative signal intensity measurements as well as optimal sensitivity, the diode should 
be fixed to the linear working range, meaning the range where the current is proportional 
to the square root of the microwave power. This happens at higher levels, greater than 1 
milliWatt. 

A microwave resonant cavity is made specially with dimensions, that is on resonance 
with the fundamental frequency of the microwave. For example, for a cavity dimensions 
of 1.7 x 0.9 x 0.45 in 3 , the fundamental frequency occurs at a microwave frequency of 
9.5 GHz. Out of resonance, no microwave would be derived. Hence, any absorption of 
the fundamental resonance of the cavity by sample will disturb the resonance frequency, 
which can in turn be detected by the proper circuitry. 

Sample size also may turn out eventually to be of particular importance. EPR/ENDOR 
signal intensity depends on the total number of spins in the active part of the cavity. Then, 
it will be better to use the largest sample tube in order to get higher signal intensity. 

For particular practice of ENDOR, only there are some other experimental equipments 
that are added to EPR. DICE ENDOR is for generating RF field. Then, a proper RF 
power amplifier is also required. In Figure I4.H| the laser part is for our future work of 
high spin polarization and selective addressing. 

The work that is explained here is summarized as following steps. 

1. We measure cw EPR and ENDOR of liquid and solid state samples to see whether 

there could be any signal from those samples, also to get sample information that 
is required for QC. The information should be composed of elaborate analyses of 
molecular magnetic properties of the open-shell entity under study. 

2. Then throughout the pulsed EPR experiment, we measure T\ and T 2 to see whether 

the sample can be a proper choice for quantum computing and quantum information 
processing for a matter of decoherence time. Recall that one of the conditions for 
a proper sample is the long decoherence time. 

3. Finally, the experiments on implementations of quantum operations, making pseudo- 

pure state and measurement issue by means of pulsed ENDOR will be explained 
and more importantly, the experiment on entanglement preparation between spins 
involved in ENDOR experiment will be discussed. 

We first start with experiments on samples in liquid phase. Samples in liquid phase, 
which enjoy rich nuclear qubits are for future usage for solution-ENDOR QC experiments. 
More elaborated examinations of the magnetic property such as linewidth of ENDOR sig- 
nals as a function of temperature, etc. are required. Then we will give some experimental 
results on solid state samples, particularly appropriate for low temperature ENDOR QC, 
for experiment with quantum entanglement. 
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Figure 4.9: Molecular structure of 2-(N-aminoxyl-N-t-Bu)-benzimidazole, BABI. 




Figure 4.10: EPR spectra of 2-(N-aminoxyl-N-t-Bu)-benzimidazole, BABI. 

4.4.1 Experiment on 2-(N-aminoxyl-N-t-Bu)-benzimidazole, BABI: 

The particular sample, which we have studied is 2-(N-aminoxyl-N-t-Bu)-benzimidazole, 
BABI, in liquid phase. The molecular structure is shown in Figure 14.91 This molecule 
is an open shell molecule with hydrogen bonding in its crystalline state. From solution 
EPR and 1 H-ENDOR spectra observed in mineral oil at ambient temperature, all the 
14 N- and 1 H-hyperfine coupling parameters of their compounds were determined. An 
EPR spectrum is shown in Figure 14*. 101 and the ENDOR spectra is shown in Figure l4*.lll 

In order to discuss electronic and molecular structures of their compounds, theoretical 
calculations based on the density functional theory were carried out. The observed hyper- 
fine coupling constants have been compared with those by the DFT, density functional 
theory, calculations, Tabl d4.3l 

The simulated EPR spectrum which is based on the calculated parameters is shown 
in Figure 14.121 Very well coincidence exists between the two EPR spectra, one from 
experiment and the other from simulation. 
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Figure 4.11: ENDOR spectra of 2-(N-aminoxyl-N-t-Bu)-benzimidazole, BABI. 



Table 4.3: Experimentally determined hyperfine coupling constants compared with those 
by DFT calculations, a/MHz. Previous work refers to []. 



Site 


This Work 


Previous Study 




Experiment 


UB3LYP / 
cc-pVDZ 


Experiment 


UB3LYP / 
cc-pVDZ 


N(l) 


1.164 


1.134 


2.162 


0.881 


N(2) 


6.040 


7.865 


5.980 


6.387 


N(3) 


28.752 


25.990 


28.525 


19.791 


H(l) 


-1.524 


-1.735 


2.471 


-1.775 


H(2) 


-2.446 


-2.192 


1.432 


-1.977 


H(3) 


0.231 


0.537 


1.039 


-0.174 


H(4) 


-3.170 


-2.843 


0.590 


-2.774 


H(5) 


0.676 


0.701 


0.309 


0.076 




Figure 4.12: Simulation of EPR spectra of 2-(N-aminoxyl-N-t-Bu)-benzimidazole, BABI. 
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Table 4.4: Spin densities from experiment and calculation. 



Site 


Experiment 


UBLYP/cc-pVDZ 


N(l) 


0.0138 


0.0240 


N(2) 


0.0717 


0.1293 


N(3) 


0.3414 


0.3494 


C(l) 




-0.0698 


C(2) 




-0.0248 


C(3) 


0.0396 


0.0368 


C(4) 


-0.0037 


-0.0120 


C(5) 


0.0514 


0.0476 


C(6) 


-0.0110 


-0.0176 


C(7) 




0.0358 


H(l) 


-0.0011 


-0.0020 


H(2) 


-0.0017 


-0.0014 


H(3) 


0.0002 


0.0005 


H(4) 


-0.0022 


-0.0018 


H(5) 


0.0005 


0.0004 







0.5026 




Figure 4.13: Spin densities on each site. 

Spin density can be calculated through McConnell's equation as follows 

a (H) = 1426.2 x p(H) (4.84) 

a (N) = 84.2 x p(N) 

a (ttC-H) = -61.8 x p(ttC-H) 

The spin densities calculated for each site is shown in Tab l d4 .41 The spin density on each 
site is shown in Figure 14.131 by circles with areas representing the magnitude. 

4.4.2 Solid State Sample Study for ENDOR QC; Malonyl Rad- 
ical 

Malonyl radicals incorporated in the single crystal of malonic acid were generated by 
X-ray irradiation at ambient temperature as shown in Figure d. 141 Spin Hamiltonian pa- 



104 



HOOC. XOOH frradiS HOOC. -/COOH 

c >~ c 

@ * tRT - (H) 

malonic acid extremely stable malonyl radical 



Figure 4.14: Generation of malonyl radical by irradiating the sample. 



Table 4.5: Magnetic tensors parameters of malonyl radical. 



s 


G X x Gyy Gzz 


a A xx /MRz °A w /MHz a A zz /MRz 


1/2 


2.0026 2.0035 2.00331 


-29 -61 -91 



rameters of malonyl radical under study, as summarized in Table 14 ,5[ have been reported 
by McConnell and coworkers. 

Figs. 14. 16114.201 show the experimental results on malonyl radical at X-band and Q- 
band. The corresponding energy levels and pulsed EPR spectra and Pulsed ENDOR 
spectra are shown to be different particularly in case of pulsed ENDOR for X-band and 
Q-band. The magnetic field for Q-band is more than three times larger than X-band. 
Then the nuclear magnetic field becomes three times lager and its correspondance relation 
with the hyperfine coupling changes and yields to a different form of ENDOR spectra. 

4.4.3 Solid State Sample Study for ENDOR QC; Diphenylni- 
troxide 

Open-shell molecular entities as molecular devices such as quantum computing and quan- 
tum information processing have emerged in interdiciplinary areas, giving us new aspects 
of molecular spinics or molecular spin science. During the coarse of implementation 
of pulsed ENDOR based quantum computing and quantum information processing, we 
have studied extremely stable diphenylnitroxide, DPNO, profiles which are necessary for 
molecular designing for the above molecular devices. 

One nice feature of DPNO, which makes it potentially a proper sample for quantum 
information processing and quantum computation is the decoherence time. Through the 
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Figure 4.15: Principal axes of malonyl acid. 
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Figure 4.16: Energy levels for malonyl, X-band. 
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Figure 4.17: Energy levels for malonyl, Q-band. 
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Figure 4.18: Pulsed EPR spectrum of malonyl radical in the single crystal, X-band. The 
arrow indicates the static magnetic field for the ENDOR measurements. 
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Figure 4.19: Pulsed EPR spectrum of malonyl radical in the single crystal, Q-band. 




Figure 4.20: Pulsed ENDOR spectrum from malonyl radical with the magnetic field set 
as in Figure l4~THl X-band 
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Figure 4.21: Pulsed ENDOR spectrum from malonyl radical with the magnetic field set 
as in Figure l4~TT)| Q-band. 
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measurement of decoherence time, it turns out that DPNO is extremely stable sample 
and even more importantly it is possible rather simply to control the decoherence time by 
controlling the dilution of the sample. Therefore we have been trying to get full analysis 
of DPNO. 

DPNO has been studied previously by several groups. To our knowledge, the first 
pioneering work has been Deguchi, 1961. This work is an ESR work both on pure DPNO 
single crystals and on DPNO diluted in benzophenone single crystals, even before the 
x-ray structural analysis appeared in 1968. The exchange line broadening was analyzed 
in this work. Then, later high resolution hyperfine EPR spectra in solution were an- 
alyzed by Yamauchi, 1967. The g tensor for DPNO diluted in benzophenone lattices 
and the hyperfine tensor of the nitrogen nucleus of the nitroxide site were determined 
by conventional EPR spectroscopy as early as 1970's by Tien-S. Lin. The principle val- 
ues for g tensor of DPNO and hyperfine tensor A for nitrogen have been reported to 
be, g xx = 2.0092, g yy = 2.0056, g zz = 2.0022, A xx = 1.9 G, A yy = 3.6 G, A zz = 23.8 G. 
Brustolon's group revisited DPNO in the benzophenone lattice by cw ENDOR spec- 
troscopy in 1988. According to their paper, the g tensor has the principle values of 
g xx = 2.0079, g yy = 2.0040, g zz = 2.0014. Only preliminary results of the nitrogen hy- 
perfine coupling and quadropole tensors in the abc system were given, though the tensor 
analysis was not given. Also, 1 H ENDOR/general-TRIPLE spectroscopy was carried out 
in isopropanol at 210 K and radical pairs were detected from the concentrated mixed 
crystals. In between the two latter introduced works, Yamauchi et al. did proton EN- 
DOR in ethylbenzene at 203 K, 1987. The substituted ortho-methyl effect was influential. 
Nevertheless, complete analyses of the magnetic tensors of DPNO from the experimental 
side have never been documented yet because of inhomogeneously broadened EPR lines 
in the solid state due to the existence of protons. 

There are several reasons on our interest on studying DPNO. We have selected 
Diphenylnitroxide as an electron-spin bus QC/QIP molecule as a central hardware for 
a molecule-based quantum spin computer. This is due to possibility of the sample for 
controlling the decoherence time. DPNO is selected also because that it is a stable 
molecular-spin physical system with well characterized magnetic tensors under strong 
microwave and radio frequency irradiations. The irradiation in this regard is sometimes 
highly required for quantum information processing and quantum computation. DPNO 
is very easy-to-dilute molecular spin in diamagnetic molecular lattices which makes it a 
quite perfect sample to control decoherence time. We have prepared mixed single crystals 
of a desired host-guest with variety of concentration ratios. DPNO is also selected as it is 
feasible to have measurements in a wide range of temperature. This point specially turns 
out to be important if we consider the critical experimental conditions for realization of 
a quantum entanglement as it is required to go down to a few Kelvin. Furthermore, it 
is also possible to extend the ENDOR based quantum computation with this sample to 
larger number of qubits as this sample is a proper sample for extension to multi-nuclear 
spin-bus systems by stable isotopic labeling and A tensor engineering. All the points 
indicated here, if collected to each other make it a valuable task to study and extract full 
analysis of DPNO. 

The accurate magnetic tensors of DPNO in the benzophenone lattice have been ac- 
quired, for the first time, in this work, by single crystal cw ENDOR spectroscopy, X- 
band, at ambient temperature. The analyses include an accurate determination of the 
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Figure 4.22: Different samples that have been used for acquiring data on magnetic 
tensors for DPNO. 
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Figure 4.23: Benzophenone molecule is isostructurally substituted by DPNO . 

quadropolar tensor of the nitrogen nucleus. In order to come over the broadening problem 
due to existence of protons, we have synthesized DPNO, DPNO-g?io and stable isotope 
labeled ones, including DP 15 NO-dio, according to the method modifying the procedure 
described in the literature, Figure l4~2"2l 

For sample preparation, DPNO has been magnetically diluted in benzophenone single 
crystals, benzophenone molecule is isostructurally substituted by DPNO at a desired 
concentration, see Figure 14.231 External form of benzophenone single crystals, as a host 
molecule, is shown in Figure l4".241 The crystal symmetry of benzophenone is orthorhombic 

o o o 

with a = 10.30 A, b = 12.15 A and c = 8.00 A, with four molecules in each unit cell. 

All the angular dependencies of EPR and ENDOR spectra were carried out in the 
crystallographic abc coordinate system at ambient temperature with a Bruker ESP350 
spectrometer (a TM mode cavity and helical coil for RF irradiation). In each crysta- 
lographic plane, two of the molecules are equivalent. Therefore, two sets of molecules 
exist while for measurements on axis, all the four molecules of a unit cell are equivalent 
and then there is not any site splitting due to different sites, see Figure 14.251 We have 
used this property in order to define the axis during our experiment. For DPNO, as 
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Figure 4.25: Single-crystal EPR spectra of DPNO, DPNO-rf 10 , and DP 15 15NO-rfi in a 
benzophenone single crystal. The horizontal axes are magnetic field in milli Tesla 
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Figure 4.26: Angular dependence of single- crystal EPR spectra of DPNO, DPNO-ciio, 
and DP 15 NO-<iio i n the crystallographic ab plane. 

will be shown later site splitting is not so clear because of the proton broadening of the 
spectra, while this fact for each plane is more clear in case of DP 15 NO-c?io. Note that 
eventhough we have made measurements at ambient temperature, the ENDOR spectra 
can be observed at other desired temperature in the range of 350 K to liquid helium 
temperature. 

For a complete analyses for finding the g tensor, each transition has been assigned 
based on the general solution for the resonance fields and ENDOR frequencies as given 
by the second-order perturbation theory. The g tensor was acquired from DP 15 NO-<iio. 
Figure l4\2fil -Figure 14.281 show the angular dependencies of single-crystal EPR spectra of 
DPNO, DPNO-dio, and DP 15 NO-g?io in the crystallographic ab, be and ca planes. 

Complete ENDOR measurement in each crystallographic plane, for all the samples 
has been carried out. Figure 14.291 shows the EPR and ENDOR angular dependency 
measurements on DPNO-rfio- Furthermore, l H ENDOR/general-TRIPLE spectroscopy 
has been carried out, see Figure [QUI and Figure l4~3*Tl 

After assigning all the transitions for the resonance fields and ENDOR frequencies 
of the observed spectra in the abc crystallographic systems for the g tensor and nitro- 
gen A and Q tensors, then the least square fitting has been done. Finally, SIMPLEX 
optimization of the observed ENDOR frequencies has been performed in the abc crystal- 
lographic systems for the g, nitrogen A and Q tensors, comparing with the theoretical 
ones obtained by the numerical diagonalization based on the hybrid eigenfield approach. 

Attempts to calculate the isotropic g-value of DPNO by quantum chemical compu- 
tations of high level have been appearing, recently. We have employed the scheme with 
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Figure 4.27: Angular dependence of single-crystal EPR spectra of DPNO, DPNO-ciio, 
and DP 15 NO-<iio in the crystallographic be plane. 
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Figure 4.28: Angular dependence of single-crystal EPR spectra of DPNO and DPNO-(iio 
in the crystallographic ca plane. 
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Figure 4.29: EPR spectrum of DPNO-<iio and ENDOR of nitrogen 14 of the same molecule 
in ab plane. Around b axis, the perturbation treatment is broken. 
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Figure 4.30: Typical ^-ENDOR and TRIPLE spectra of DPNO-/ii in a benzophenone- 
h w single crystal. Arrow indicates the second frequency for pumping. 
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Figure 4.31: Typical 1 H-ENDOR and TRIPLE spectra of DPNO-Ziio in a benzophenone- 
h\Q single crystal. Arrow indicates the second frequency for pumping. 

Table 4.6: g tensor obtained from the experiment as compared with the calculated one. 
*:ADF2005.01, basis set: QZ4P, LDA: VWN and XC: Becke Perdew. 





giso 


gxX 


Syy 


gzz 


Experiment 


2.0069 


2.0110 


2.0065 


2.0033 


Calculation 


2.0062 


2.0109 


2.0057 


2.0021 



different basis sets in order to calculate the parameters of this particular sample. Table 
!4.6l shows the experimental and quantum chemical calculational results for g tensor values 
and Table l4~B1 shows the nitrogen A and Q tensors in comparison to the calculated ones. 

For this results, the molecular structure in accordance to the crystallographic axis is 
shown in Figure 14*. 321 

For theoretical calculation of g tensor, various correction terms have been introduced 
to the isotropic g value of the free electron. 



.Rel 



As 



„Dia 



is the g value of the free electron, 



„Para 



(4.85) 



correction including Aj 



r KE 



Ag 



MV 



etc. 



. = 2.0023. Ag Kel ~ -0.0002 is relativistic 
Relativistic correction is tiny here as we are 



dealing with nitrogen, proton and carbon, that none of them has a large atomic num- 
ber. Anisotropic corrections, diamagnetic correction, Ag Dm ~ 0.0001 and paramagnetic 
correction, Ag Pam ~ 0.004 also have been included. Diamagnetic term is small due to 
considerably small nuclear gyromagnetic ratio as compared to the paramagnetic term. 
The latter one is more influential and is calculated as follows for molecular orbital (MO) 
spins 

Ag^ = ^e IWfc)l , (4.86) 

k 6 
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Figure 4.32: Molecular structure for diphenylnitroxide in Benzephenone. 



Table 4.7: A and Q tensors Obtained from the experiment as compared with the cal- 
culated ones for different basis set. Also the spin densities on nitrogen and oxygen have 
been shown. 



Basis 




6- 

31g(d) 


6- 

311g(d,i 


EPR- 
,)II 


EPR- 
III 


DZVP 


dzvp: 


5 TZVP 


Exp. 


Set 




















A/ 


XX 


56.195 


51.386 


55.084 


58.360 


66.592 


68.643 


57.490 


68.16 


MHz 


yy 


3.744 


-5.107 


-1.665 


-2.037 


6.488 


7.018 


-2.366 


8.55 




zz 


2.830 


-6.097 


-2.430 


-2.731 


5.796 


6.362 


-3.121 


6.56 


A iso /Ml 


I 


20.923 


13.394 


16.996 


17.864 


26.292 


27.341 


17.335 


27.13 


^■aniso/ 


XX 


35.272 


37.992 


38.088 


40.496 


40.300 


41.302 


40.155 


40.87 


MHz 


yy 
















-19.33 






17.179 


18.501 


18.661 


19.901 


19.804 


20.32 


19.700 






zz 
















-21.53 






18.093 


19.491 


19.426 


20.595 


20.496 


20.97 


20.455 




Q/ 


XX 


-2.481 


-2.744 


-2.758 


-2.610 


-3.110 


3.344 


-2.820 


-2.48 


MHz 


yy 


1.834 


2.054 


2.043 


2.014 


2.245 


2.434 


2.160 


1.89 




zz 


0.647 


0.690 


0.715 


0.596 


0.866 


0.910 


0.660 


0.59 


P 


N 


0.331 


0.317 


0.333 


0.295 


0.327 


0.322 


0.328 




P 





0.495 


0.497 


0.478 


0.497 


0.486 


0.484 


0.482 
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Figure 4.33: Contribution of the paramagnetic terms to g tensor of DPNO; idea for g 
tensor engineering. Spin densities on each spin site is shown in colors, red for positive 
spins and blue for negative spins. Radiuses are for the magnitude of the densities. 

where £ is the spin-orbit coupling constant. 

Then there is a notable point here that with quantum chemistry, it is possible to 
control the g tensor through location of the excited states, g tensor engineering. For 
ENDOR quantum computation and even more elaborated schemes such as ELDOR, 
Electron-Electron DOuble Resonance quantum computation and specially the case for the 
larger number of qubits, it is highly required to introduce spins with different g tensors 
as to be distinguishable from each other. Then the concept of g tensor engineering turns 
out to be very important. For this particular molecule, the status is shown in Figure 
14.331 It is then clear that a transition of a /3 spin gives a considerable contribution to the 
g tensor. This effect can be used for engineering spins with different tensorial values. 

We emphasize here that the computational chemistry is very important specially for 
g tensor engineering. Although the results reported in the above table compared to the 
experimental values have been selected after calculations with several basis sets. The 
proper basis set then has been selected as the one which yields to values closer to the 
experimental results. But the basis as reported here as the proper basis set, can not be 
generalized as it is special case of the molecule and for other samples the situation might 
be reported to be completely different. For theoretical calculations of A and Q tensors, 
we have considered the corresponding equations. 

We are ready for measurement of DPNO for pulsed ENDOR QC. Experiment has 
been done at X-band and Q-band, at different temperatures. Energy levels for DPNO 
at Q-band is shown in Figure 14.341 Figure 14.351 shows the pulsed EPR measurement of 
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Figure 4.34: The energy levels for DPNO. 



Table 4.8: NMR frequencies of various nuclei. 



Nucleus 


v n at 9.5 GHz (3390 G) 


v n at 35 GHz (12,488 G) 


X H 


14.4 MHz 


53.2 MHz 


2 H 


2.22 


8.16 


14 N 


1.04 


3.83 


15 N 


1.46 


5.39 


19p 


13.6 


50.0 



DPNO at Q-band and 50 K. Pulsed (Davis) ENDOR spectrum of DPNO is shown in 
Figure 1061 

We did also experiment with DP 15 NO— Pulsed EPR is shown in Figure 14.371 and 
pulsed ENDOR (Mims) is shown in Figure 14.361 This molecule provides three qubits 
and to some extent proves partly our claim that this molecule is an electron spin bus 
molecule and is possible to increase the number of qubits for this molecule. However, in 
Figure H.38I the spectrum is not so clear in comparison to the previous ones. The reason 
for this is also easy to understand. See the table 14.81 that shows nuclear frequencies for 
X-band and Q-band, for some spins. Nuclear spin frequencies of 15 iV is much smaller 
than e.g. proton. Therefore, in order to manipulate 15 N we need a very high power rf 
amplifier, more than two times of the one that we had for this experiment. Unfortunately, 
at the time of the experiment, the required rf amplifier has not been available. 

In order for quantum computation, decoherence time profiles for several sample have 
been measured by pulsed-based EPR. Spin-lattice relaxation time, T\ for DPNO is 
392/msec and for DP 1 5NO-c? 10 is 42.5/ msec, at 10K. measurement on spin-spin relaxation 
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Figure 4.35: Pulsed EPR measurement of DPNO. 
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Figure 4.36: Pulsed ENDOR measurement of DPNO. 
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Figure 4.37: Pulsed EPR measurement of DP 15 NO-di . 
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Figure 4.38: Pulsed ENDOR measurement of DP 15 NO-d ; 



time T 2 gives results for DPNO to be 0.777//isec and for DP^NO-dio to be 0.489///sec. 
These values are quite fine for introducing quantum computation and information with 
this sample. 

4.5 Decoherence Time Measurement 

Decohrence time measurement is done with pulsed EPR. Here, we give a brief introduction 
on pulsed EPR as it is also preliminary step to pulsed ENDOR and will be used frequently 
for pulsed ENDOR based quantum computation. 

Remind that NMR most generally is performed with employing a pulse FT technique. 
The reasons are simple as NMR is slow and mostly it contains very narrow lines spread 
over a wide frequency range, as compared to the linewidth. Then NMR in its pulsed 
performance because that it reduces the measurement time by moving to a pulse FT 
technique gives essentially very nice advantages. However, EPR spectra is wide. Then, it 
seems useless to move to pulsed EPR scheme. Although, it is still some more advantages 
through pulsed scheme compared to continuous wave scheme. In a pulsed performance it 
is possible to get more information by measuring the time domain and using the multi- 
dimensional techniques. 

With pulsed EPR e.g. it is possible to measure relaxation times more directly as 
compared to continuous wave EPR. With several pulsed techniques, variety of informa- 
tion can be accessed in a particular experiment. For example, ESSEM, Electron Spin 
Echo Envelope Modulation gives information on interactions of electron spin with the 
surrounding nuclei. 

Pulse EPR spectroscopy is performed by applying a short microwave pulse, less than 
20 ns, but intense, more than 30(W. Then the microwave signals which have been 
generated by the sample is measured and Frourier transformed and gives the experimental 
results. 

Let us assume the laboratory frame as usually is defined with the magnetic field, B , 
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parallel to the z axis, the microwave magnetic field, Bi, parallel to the x axis, and y axis 
to be orthogonal to the x and z axes. Same as the processing which goes on in NMR, 
here also the physical background of the system can be explained simply by assuming a 
spin which is placed in a magnetic field. For EPR, a torque is exerted on the electron spin 
and causes the magnetic moment to precess about the magnetic field with the angular 
frequency named as Larmor frequency defined as follows 

u L = - 7 S , (4.87) 

where ojl is the Larmor frequency. For a free electron which has a ^ /2n value of approx- 
imately —2.8 MHz/Gauss, the Larmor frequency is about 9.75GHz for a field of 3480 
Gauss. Suppose that a large number of electron spins are placed in a strong magnetic 
field, i?o- Then the electron spins are precessing about the magnetic field, however their 
orientation should be completely random in xy plane, since there is not any preferable 
direction in this plane. Now, suppose that a rotating magnetic field, B ± , is applied, let's 
say along the +x direction. B\ is also assumed to be much smaller than B . In order 
to study the effect of B\ on the electron spins, we would be better to move to a frame 
rotating with the microwave frequency. Then, in this frame B\ looks as if it is a station- 
ary field. The assumption here is also that we are working at resonance, meaning that 
uj L = uq, where c^o is the microwave frequency. 

Resulting effect from above is that the magnetization would be precessing about B\ 
at a frequency as follows 

Ul = -7^!, (4.88) 

uj\ is called Rabi frequency. The magnetization will rotate about Bi, +x as assumed 
here, as long as the microwave are applied and the rotation angle is calculated as follows 

a = 7|Si|*p. (4.89) 

The angle a is called the tip angle and t p is the pulse length. It is then clear from the 
bove equation that the tip angle depends on both the magnitude of B\ and the length of 
the pulse, in other words how strong pulse for how long is applied on the sample. 

Now let's assume the original problem. For an on resonance condition and B\ along +x 
direction, the magnetization after a n/2 pulse would have a direction along — y direction. 
In laboratory frame, as introduced above, the magnetization then has a rotating motion 
in xy plane at the Larmor frequency. This rotation actually generates current and voltage 
in the resonator and gives a signal which is called FID, free induction decay. 

The condition above for resonance, which means that the Larmor frequency being 
equal to the microwave frequency, clearly can not be satisfied generally as there are 
many different frequencies in an EPR spectrum. Therefore, most the cases experience 
nonresonance condition. Then, the magnetization after applying the microwave pulse 
would have a rotation in xy plane which might be added up to be canceled partly by the 
rotating field B 1 . The net effect is chosen as follows 

Au = uj - uj . (4.90) 

If Auo > then the magnetization is gaining and will rotate in a counter-clockwise 
fashion. Conversely, if Au < then the magnetization is lagging and will rotate in 
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Figure 4.39: Decoherence time has been measured for different samples. 

a clockwise fashion. Then, different resonance frequencies in the EPR spectrum are 
discriminated as they have different Auj and this is encoded in the FID of the sample. 

The above discussion, however has been realistically true and in real world it can 
not be fulfilled completely. Magnetization after applying any pulse does not remain in 
a particular direction for a long time, for reasons as the spins pertain interactions with 
each other and with the environment. This effect is determined with two constants, T\ 
and T 2 . By definition, T\ is the spin lattice relaxation time and describes how quickly the 
magnetization returns to alignment with the z axis. Also, T 2 is the transverse relaxation 
time and describes how quickly the magnetization in the xy plane disappears. 

Thermal equilibrium is for before applying any pulses. After applying a pulse, then 
there will be no longer the thermal equilibrium. The equation which explains the evolu- 
tion is then explained as follows 



The above equation is for the case that we apply a it pulse and it will be slightly dif- 
ferent for other pulses, e.g. tt/2. therefore, the rate constant at which the magnetization 
recovers to thermal equilibrium is T\ and can be calculated as above after getting the 
signal which has to be taken by the FID of the corresponding experiments. 

We did measurement of T\ for some particular samples, Figure 14.391 The pulse se- 
quence for measurement of T\ is shown in Figure 14.401 The acquired spin lattice relax- 
ation time profiles are given in Table 14.91 The experimental detected saturation recovery 
curves are given in Figure 14.411 T± is temperature dependent. Figure 14.421 shows the 
dependencies for particular sample of our interest. 

Measurement of T 2 also can be understood in almost the same manner. The spins 



m z (t) = tuq x (1 




(4.91) 
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Figure 4.40: Pulses for T\ measurement with saturation recovery. 
Table 4.9: Spin-lattice relaxation times at 10 K. 



Radical 


Ti/mS 


Malonyl radical 


91.5 


DPNO 


392.0 


DPNO-d 10 


42.5 


Alanyl radical 


76.4 


Cystinyl radical 


2.2 




40 

t/ms 



80 



Figure 4.41: Saturation recovery curves at 10 K. Malonyl radical: green, DPNO: yellow, 
DPNO-dio: purple, Alanyl radical: blue and cystinyl radical: red. 
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Figure 4.42: Spin-lattice relaxation times as a function of temperature. Malonyl radical: 
green, DPNO: yellow, DPNO-dio: purple, Alanyl radical: blue and cystinyl radical: red. 



122 




jl/2 7T Electron spin echo 



Figure 4.43: Pulses for T 2 measurement with two-pulse echo decay. 
Table 4.10: Spin-spin relaxation times at 10 K. T 2 for malonyl is given at 20 K. 



Radical 




Malonyl radical 


5.200 


DPNO 


0.777 


DPNO-rfio 


0.489 


Alanyl radical 


1.560 


Cystinyl radical 


1.267 



interact with each other, resulting in mutual and random spin flip-flop. Molecular motion 
can also contribute to this relaxation. These random fluctuations contribute to a faster 
fanning out of the magnetization. The decay of the transverse magnetization or FID from 
this mechanism is in general exponential with relation to the spin-spin relaxation time, 
and then gives a way to measure T 2 , through the following equation. 

m_ y = e {t/T2) . (4.92) 

The pulse sequence to acquire information of T 2 is given in Figure I4.43| with which 
threshold for T 2 of samples are given in Table I4~TU1 Detected experimental results of echo 
decay curves are shown in Figure 14*. 441 and temperature dependencies in Figure H.451 T 2 
is not temperature dependent and this fact is clear from our experimental results, too. 



4.6 Generation of a Pseudo-pure State for ENDOR 
Quantum Computing 

In our experiments, pseudo-pure states are acquired by applying some particular pulse 
sequence and waiting time. For instance, we give an example on how to get a pseudo-pure 
state. Initially, the Larmor frequency of the nuclear spins has been ignored due to that 
it is sma ller than the Larmor frequency of the electron spins. For making pseudo-pure 
state , (jl20[ ). the initial state is the truncated Boltzmann density matrix, p' = — |S Z , 
from the original Boltzmann spin density matrix for high temperature as pe = (1 — 
-Kb)|I4 + Kb ■ Pb, with Kb = Pb-Bo/3A;b^ for g = 2 and p p = 1/4I 4 — 3/2S z . Pseudo- 
pure state is prepared by providing some particular microwave and radiofrequency pulses 
in addition to a waiti ng t ime for making transverse components decay for equalizing 
different populations, (120). 

In a rather similar approach as stated above, we implemented pseudo-pure state in 
an experiment. Figure l4*.47l gives the initial state in comparison to the pseudo-pure state 
after applying the pulse sequence Figure 14.461 for making pseudo-pure state. It is clear 
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Figure 4.44: Echo decay curves at 10 K. Malonyl radical: green, DPNO: yellow, DPNO- 
d\Q. purple, Alanyl radical: blue and cystinyl radical: red. 
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Figure 4.45: Spin-spin relaxation times as a function of temperature. Malonyl radical: 
green, DPNO: yellow, DPNO-^xq: purple, Alanyl radical: blue and cystinyl radical: red. 
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Figure 4.46: The pulse sequence that has been used for generation of pseudo-pure state 
as shown in Figure E. 4 71 
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Figure 4.47: ENDOR measurement before and after applying the pulse sequence for 
making pseudo-pure state. Pseudo-pure state is originated after applying the pulses. 



that the state corresponding to u 34: is suppressed giving us the pseudo-pure state on the 
state corresponding to u>i 2 - We should note that in the pulse sequence Figure [OBI the 
first 109 deg. microwave pulse of co>24 and n/2 radiofrequency pulse of u>u are for pseudo- 
pure state generation and the other pulses are applied for (Davies) ENDOR detection. 
Pseudo-pure state to the other energy levels also would be acquired by applying some 
other form of pulse sequences. More details on how to make different pseudo-pure states, 
by interchanging the rf and microwave pulses, would be clear through our experiments 
on ENDOR quantum computation in the following parts of this chapter. 



4.7 Universal Set of Gates with ENDOR 

Quantum gates are applied by Microwave pulses on electron spin and radiofrequency 
pulses on nuclear spins. For making a universal set of gates, an arbitrary one-qubit gate 
and a non-trivial two-qubit gate are required. One-qubit gate is easy to be performed in 
pulsed ENDOR quantum computation. For instance, NOT-gate is done with a n pulse 
applied in e.g. x direction, Figure l4~4*%l 

Hadamard gate is performed, up to an overall phase, with a ir/2 pulse in y direction 
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Figure 4.48: Not-gate is performed by applying a it pulse. 
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Figure 4.49: The pulse sequence for making a Hadamard gate, up to an overall phase. 

followed by a 7r pulse in x direction, Figure 14.491 

CNOT gate is also performed by applying a 7r radiofrequency pulse in x direction, 
Figure 14*. 501 There would be a phase in this performance for making CNOT gate, however 
as far as we are working in a sequence of gates, for our special interest quantum circuit 
for making entanglement, this would be possible to eliminate the effect of this phase, as 
we did and results are shown in further coming experiments. Also, it is possible to apply 
some more complicated pulse sequences that will give the exact form of the CNOT gate. 
The corresponding part on NMR QC, in Chapter 3, has been useful in this regard. For 
example, the pulse sequence that has been introduced for making CNOT, lr>.12| up to a 
total phase, is essentially applicable for pulsed ENDOR QC, too. 

Experimental evidences on the establishment of different quantum gates will be clear 
in the last part of this chapter that we present the experiments on making entanglement 
with ENDOR and implementation of quantum operations for superdense coding. 
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Figure 4.50: The pulse sequence that we have used for making CNOT in the experi- 
ment for making entanglement. For a general CNOT gate the pulse sequence is more 
complicated. See text. 
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Figure 4.51: A Hahn echo. 

4.8 Readout of Quantum Computation with ENDOR 

Readout in ENDOR experiment is done conventionally through electron spin echo. We 
use the conventional scheme for detection of states in addition to other scheme that will 
be explained in this section. First, we give a brief introduction to the electron spin echo 
detection in ENDOR. Then, we will explain additional detection scheme, TPPI, time 
proportional phase incrementation for quantum computation. 



4.8.1 Electron Spin Echo: 

After that we apply a microwave pulse, it produces a signal that decays away and gives the 
FID. However, if the spectrum is broadend enough it is possible to recover this effect by 
applying another microwave that produces Hahn echo. The processing for echo is simply 
represented in Figure 14*. 511 Here, the transverse relaxation leads to an exponential decay 
in echo height as follows 

-2r 

echo height ~ e T 2 . (4.93) 



Result of the measurement is expressed with the echo intensity that is the signal 
intensity detected after applying the pulse sequence that includes a | followed by a n 
pulse in x direction in addition to a final echo pulse in y direction. Detection through 
electron spin echo is the generally used measurement scheme for any pulsed ENDOR 
experiment. We have used almost the same approach. Information on the nuclear spins 
might exchange and be detected through electron spin echo. However, for some special 
cases, e.g. detection of entanglement that a global information on the state of the two 
spins is required electron spin echo may not be completely adequate and rather more 
elaborated scheme is required that is explained in the following section. 

4.8.2 Time Proportional Phase Incrementation, TPPI: 

TPPI is known in ENDOR for separation of multiple quantum coherences. However, we 
used the scheme in a way that it can give a global information on the state of electron 
and nuclear spins. This is for resolving the restriction of ENDOR spectroscopy that any 
measurement can only be done via electron spin echo. Specially for the case that infor- 
mation on the state of entanglement is required, TPPI is an essential part for detection 
of resolving the entangled state from the simple superposition states. 
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In TPPI, the phase shifts are implemented by incrementing the phase frequency of 
the individual detection pulses in consecutive experiments. What we mean by detection 
pulses will be clear after that we represent our experiments for entanglement with pulsed 
ENDOR. Meanwhile, this term might be understood if we recall the part for entanglement 
witness derived from NMR SDC, see chapter 3. We introduced some unitary back opera- 
tions for transferring the state to the computational state. Hereafter, the corresponding 
unitary back operation would be called detection pulses. 

Firstly, we apply detection pulses with arbitrary phase frequencies, Vj on jth spin. 
Therefore, we have 

Acjj = 2nAv j . (4.94) 
Then consecutive experiments are performed for a time At. We have 

= AujAt. (4.95) 

The phase shifts are detected through the experiment and give information on the 
state. Let us give an example. Suppose that the state in the experiment, is a Bell state 
of the following form 

M = -^(|01)-|10». (4.96) 

Detection pulses that we apply are a microwave n pulse followed by a radiofrequency 
7r/2 pulse, with phases of 0i and 2 , respectively. Then the detection unitary operation 
might be written as 

U d = u M (it/2, 2 )cu 24 (vr, 00. (4.97) 

The first pulse is a microwave n pulse (0^4(71")) and the second pulse is a 7r/2 radiofre- 
quency pulse (a>34(7r/2)). 

Then, measurement is performed with electron spin echo. The intensity would be as 
follows 

= tr{SfU d \i,){i,\Ul] (4.98) 
= ~ cos^- (f) 2 )}. 

Therefore, through detection by TPPI we would get the phase shift as follows 

0i - 02 = A(wi - uj 2 )At. (4.99) 

Here, for the angular frequencies, u>j, we have 

A(cji-cj 2 ) = 2ttA(z/i - z/ 2 ), (4.100) 

Uj is the phase frequency on the jth spin. Detection of a pulse with a phase shift of 0i — 2 
demonstrate that the state has been \ip) and not the separate states of spins. Because, 
if the state under measurement is a separable state of the spins, after applying detection 
unitary operations with v\ and phase frequencies on spins then TPPI measurement 
would give results as only <pi and 2 and not 0i ± 2 . 
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4.9 Pseudo-entanglement with END OR at Q-band 
and Low Temperature 

In the following part we will give report on our experiment as how the phase control has 
been acquired and furthermore the results will show nicely a strong evidence on realizing 
the (pseudo) entanglement. Higher magnetic field and lower temperature and probably 
other experimental properties would give quantum entanglement. 

In our experiment, spins have been selectively addressed by their resonance frequen- 
cies. We did experiment either with pseudo-pure state or without pseudo-pure state. 
Detection has been performed by electron spin echo while phase manipulation has been 
done by TPPI. 

The experiment that is explained here is performed on two different samples. The 
first choice has been malonyl radical. The other sample, has been DPNO, the molecule 
which we have extensively studied and believe that can be a proper electron spin bus 
molecule for quantum computation. For both the cases, we got some reasonably good 
results on implementation of (pseudo-)entangled state. However, we emphasize the spin 
manipulations by means of phase incrementation. This is quite useful point for further 
progress in this field. 

Pseudo-pure state has been made in a similar manner as discussed in the previous 
section and the final state of the system is then |10), after a 109.47° microwave co>24 
followed by a I radio frequency Ui 2 . Slightly different pulse sequences give the other 
choices for the initial state. Then, the initial state changes into a Bell state if pulses are 
applied as follows 

(34), s 1 

|10) -^_(|10> + |11» (4.101) 

( 24 ) I \ 1 

_(|io)-|01». 

This is equivalent to have a Hadamard gate on the nuclear spin followed by a CNOT gate 
on the state of the electron and nuclear spins. 

Then, detection is done by unitary back operations. Entanglement is a global property. 
Then, a single local measurement as opposed by electron spin echo would not generally 
give any information on the status of entanglement. Therefore, we have used the phase 
incrementations on both microwave and radio frequency to extract the required informa- 
tion. We apply phase dependent pulses, 7r microwave u 2 4 with phase 0i followed with a 
7r/2 radio frequency co> 34 with phase 2) that gives the phase dependent echo intensity of 
the following form 

I = -i[l-cos(0i-0 2 )]. (4.102) 

To be noted here, in case that the initial state has been some other state, then the 
other Bell state would be generated after applying the corresponding pulses for making 
entangled states. After the phase dependent pulses for detection, we would get other com- 
bination of the phases. This is shown in Table l4~TT1 for several choices of microwave and 
radiofrequencies for pulse sequence for making entanglement, Figure l4*.53| with resonance 
frequencies shown in Figure E.521 
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Figure 4.52: Energy levels for a system composed of two spin 1/2. The transitions naming 
are in accordance to the ones indicated in the pulse sequence for making entanglement 
Figure 14.531 
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Figure 4.53: The pulse sequence that has been used for generation of entanglement state. 



Table 4.11: Different Bell states derived by different pulse sequences. 
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Also, we would like to emphasize that the scheme which has been explained here is for 
detection of the entangled state from the simple superposition of the state. On the other 
hand, we always get one out of the two possible combinations of the phases and this is not 
sufficient to fully determine the particular Bell state. However, it is sufficient to determine 
the entanglement of the state. Then, the other matter, meaning that determination of 
the particular Bell state would be possible by considering the initial state and the pulses 
which have been applied in order to make the entangled state. 

The artificial phase frequencies in our experiments are Aujj = 27rAz/j, V\ = 1.0 MHz 
and U2 = 5.2 MHZ, as arbitrary values. The resultant spectrum is a 2D spectra, the TPPI 
frequency against time. See Figure l4".54l for an example on TPPI detected 2D spectra for 
different pulse sequences i.e. different Bell states. 

Then the phase interferogram against time is Fourier transformed. See Figure 14.551 
for an example on an interferogram which corresponds to the TPPI detected spectra of 
Figure E. 54( a). The combination of the phases, whether to be addition or subtraction, 
gives evidence on the phase of the Bell state and then proves the entanglement of the 
state. See Figure 14*. 561 for experimental detection of (pseudo-) entanglement for malonyl 
radical. 

The same experiment is done for DPNO. Two-dimensional TPPI frequency against 
time spectra are shown in Figure l4~. 571 Furier transformed experimetal results are given in 
Figure 14.581 The experiment with DPNO involves three qubits. However, because of the 
reasons that we explained before on low power radiofrequency amplifier, nitrogen could 
not manipulated. Therefore, the established entanglement is a bipartite entanglement as 
between proton and electron spin. Nitrogen is a basis for descriminating proton energy 
levels. However, we are planning to use this molecule for more elaborated quantum 
computation after providing a high power radiofrequency amplifier. 

4.10 Implementation of Quantum Operations for Su- 
perdense Coding; X-Band 

Superdense coding introduced by Bennett and Wiesner is a non-local quantum algorithm 
in which two classical bits of information have been transformed from Alice to Bob by 
sending only a single qubit. The scheme, explained in chapter 3, is based on the fact that 
the entangled initial state has been shared between the two involved parties. See chapter 
3 for explanation of SDC and corresponding materials. 

Concepts for ENDOR-based experimental setup for two-qubit superdense coding are 
depicted in Figure I4.59| in which S and I denote the electron-spin part and nuclear spin 
one. The quantum circuit for SDC, Figure 14.591 includes a Hadamard and controlled 
NOT (CNOT) gates. [/, stands for one of the unitary transformation for encoding; out 
of the choices {X, Y, Z, I}. The first Hadamard and CNOT gates generate an entangled 
state between the electron and nuclear spins. Following the unitary transformation Ui, 
the CNOT and Hadamard gates (unitary back operations) the entangled state extract 
the encoded information. Except for a phase factor which dose not affect any signal, 
selective ir/2 and tt pulses are available for the Hadamard gate for single qubit and the 
CNOT gate for two-qubit in pulsed magnetic resonance spectroscopy, respectively. 

Here we give a report on the implementation of the quantum operations for SDC 
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Figure 4.54: 2- dimensional spectra, TPPI frequency against time for detecting the spin 
echo intensity through phase incrementing, in consecutive experiments on malonyl. The 
interferogram of the spectra shown here is Fourier transformed to get the information on 
the individual phases as well as the combination of the phases, in case of existence of any 
entanglement. Particularly, for a) this process is explained in more details in text. Also, 
see Figure 14351 
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Figure 4.55: This figure is to show rather more in details the scheme which has been used 
for detection and analyse of the phases in order to determine demonstration of (pseudo- 
entanglement, a) is the interferogram (one slice) of the TPPI as shown in Figure l4~.54f a). 
b) is the Fourier transform of the spectrum shown in (a). 
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Figure 4.56: Experimental results of demonstration of (pseudo-)entanglement as the 
Fourier transform of the phase interfrograms. Phase frequencies are v\ = 1 MHz and 
v-i = 5.2 MHz. The bottom two correspond to the experiment on malonyl radical, while 
the upper one is due to the experiment on DPNO. 

by ENDOR. We do not call this experiment as an implementation of SDC because that 
the entanglement is actually pseudo-entanglement. If the experimental conditions, as 
discussed to be just in hand, are provided as it will be in a near future, then the same 
experiment with similar processing would be given as an implementation of SDC with 



Here, the main idea is testing the ENDOR system for QIP giving a testing ground for 
QIP and QC to molecule-based ENDOR. In our experiment, states have been prepared 
as pseudo-pure states in a similar approach as explained before in this chapter. 

The pulse sequence that has been used is represented in Figure I4.6UI There are three 
main parts of the sequences, i.e., the preparation of the pseudo-pure states, manipulation 
and finally detection according to the customs in magnetic resonance spectroscopy. 

The outermost left part labeled by "A" is for preparation of the pseudo-pure state. 
Two pulses on electron and nuclear spins, with additional waiting times in order to make 
the off diagonal term of the density matrix vanishing, are required to acquire the pseudo- 
pure state. The first two pulses in the central part of the sequence are for entangling 
and also the last two pulses are for detection of the entanglement. Two phases of <pi and 
02 for the pulses in the detection part are required for clearing up the entangled states 
from the simple superposition states explained in previous sections. In the central part 
between the entangling and detecting the entanglement one of the qubits, nuclear spin 
in our experiment, is encoded by randomly applying one of the four pulses of I, X, Y, Z. 
The necessary pulses for encoding are described in Table 14.121 

Finally, there are pulses for detection by an electron spin echo. For the measurement 
part, the situation has been modified for some detection considerations. This considera- 
tions have been much improved and modified for our following experiment at Q-band and 
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Figure 4.57: 2Dimensional spectra, TPPI frequency against Time for detecting the 
echo intensity through phase incrementing, in consecutive experiments on DPNO. 



135 



1 

1 



+i + H-i-» 




^ = -5.2 MHz 
v z = 1 .0 MHz 



B„ = 1 . 2044T 



4 6 

TPPI Frequency/ M Hz 



10 




4 6 

TPPI Frequency/ M Hz 



B n = 1.2066T 



1 



Figure 4.58: Experimental results of demonstration of (pseudo-)entanglement as the 
Fourier transform of the phase interfrograms of DPNO. 




Figure 4.59: Quantum circuit implementing SDC. 



Table 4.12: The unitary operation and corresponding pulse sequences for encoding. 





Initial State 


Necessary Operation 


Required Pulses 


Encoded State 


U=I 

u=x 

U=Y 
U=Z 


7l(|00) + |H)) 
7l(|00) + |ll)) 
73(100) + |11)) 
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exp(-iirl x ) 
exp(— iirly) 
exp(-iirl x ) 


Pf(^)Pi 2 (vr) 
Pf(2vr)Pf(vr)Pi 2 ( 7 r) 


73(100) + |11)) 
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^(|01>-|10» 
75(100) -|11)) 
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Figure 4.60: Pulse sequence for implementation of superdense coding by molecule-based 
ENDOR. 



Table 4.13: Detection through the angular dependent electron spin echo intensity. 





Angular dependant radio fre- 
quency pulses for encoding 


Detected angular dependant 
echo intensity 
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Jq(— 3cOS# — COS0 + 
4 COS | COS 02 COS [01 — 02]) 


u=z 


mo) 


^(-1 3cos6> + 
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Figure 4.61: Detection through the angular dependent pulses. This scheme for detection 
has been used in case that full analysis based on arbitrary change on microwave and ra- 
diofrequency phases have not been available. Though later the scheme has been improved 
as it is reported in text, following sections. 

low temperature. There, the phases of microwave and radio frequecies could be changed 
as aribitrary. Then detection has been done by TPPI scheme as will be explained later. 
But here we would like to represent the first proposal on detection. In this study, we 
have used the electron spin echo detection. The echo intensities have been detected for 
different angular dependencies of the pulses in the encoding part. As a result, there are 
four sets of angular dependencies for radio frequency pulses which have been used for 
encoding, (see Table l4~T3^) . see Figure HUTl 

Table 14.1 31 shows the angular dependence of the electro-spin-echo intensities in the 
detection part, depending also on experimental conditions for measurements. In the 
SDC experiments, the detected echo intensities incorporate the terms characteristic of 47r 
period, as shown in Table 14.131 We have detected this salient behavior of the intensities, 
exemplifying the case for Ui = X as given in Figure 14.621 In contrast to the 2ir period of 
the population, the observed 47T period originates from the spinor property as the intrinsic 
nature of spins under study with the experimental condition of the selective microwave 
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Figure 4.62: Air period dependence of the electron spin echo detection scheme. Microwave 
frequency: = 9.38725 GHz. Static magnetic field set for ENDOR: B = 335.918 mT. 
Temperature: T = 20 K . 

excitation. The sign difference appearing between Figure 14.621 and Table 14.131 is due to 
the difference between the experimental setup. Therefore, we have introduced pulsed 
ENDOR based approach to make quantum operations, pseudo-pure states and detection 
for SDC. 

4.11 Conclusion 

Our efforts have been made to develop entangling unitary operations in order to check the 
credibility of the molecule-based ENDOR QIP from both the theoretical and experimental 
sides. Angular dependent electron-spin echo intensities in the readout measurements have 
been derived in simple analytical expressions for SDC-ENDOR experiments. TPPI has 
been used for entanglement detection. 

Regarding the concept of entanglement, this work is muc h clo ser to the existence 
of quantum entanglement as compared to the previous work, ([l20h . In addition to the 
experiments at X-band and Q-band, we did also W-band at very low temperature, around 
3.3 K. Under the experimental conditions, the distance to the quantum entanglement is 
much less than the one that has been reported previously. 

The other important fact is the proper sample for quantum computing with ENDOR, 
that can give even quantum entanglement at low temperature and strong magnetic field. 
The previous work has been with malonyl redical. However, we have studied DPNO, 
heavily, and showed that for high magnetic field and low temperature this sample is the 
proper sample. Malonyl radical is not workable at this experimental conditions due to 
very long relaxation time. It does not give any clear ENDOR spectra. However, DPNO 
that has been used for very wide range of experimental conditions has been workable and 
gives very nice ENDOR spectra. 
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Also, about DPNO, it is possible to increase the number of qubits, rather simply by 
isotope labeling and in this regard we did several sample analytical experiments. Nitrogen 
and flourine labeled samples of DPNO have been synthesized, for instances. This sample 
also has been used for liquid state study quantum computing, due to the fact that it gives 
very clear ENDOR spectra for any purposes. 

With ENDOR, we have been throughly engaged in a complete establishment of quan- 
tum computing. This work is still going on, however we can say that the preliminary 
experimental results are well representing the feasibility of the ENDOR quantum infor- 
mation processing. 
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Chapter 5 
Conclusion 



We have studied the concept of quantum entanglement for quantum information process- 
ing and quantum computation. Through the theoretical study, we conclude that specially 
for exponential enhancement, entanglement is a prerequisite quantum property. However, 
entanglement for ensemble quantum information processing and quantum computation is 
not so clear. Firstly we studied implementation of superdense coding by means of nuclear 
magnetic resonance, NMR. We have proved that NMR with low spin polarization, is not 
a quantum system that can give a quantum enhancement over classical processing. For 
instance, it has been proved that NMR implementation of superdense coding is erroneous, 
in the sense that the acquired experimental results are in reach with any classical system. 

We have also introduced a new measure of quantum entanglement, system oriented 
entanglement witness, that is specially applicable for NMR system. This measure of entan- 
glement has been shown to have several advantages over any conventional entanglement 
witness. Firstly, it is possible to detect entanglement in a single run experiment and with 
a simple to implement quantum operations. This measure of entanglement is also useful 
for other ensemble system such as electron nuclear double resonance, ENDOR. 

ENDOR has been introduced for quantum computation and it has been shown that 
the system is more appropriate because of the existence of an electron spin. The electron 
spin polarization, under the same experimental conditions, is three orders of magnitude 
larger than nuclear spin polarizations. This makes the required conditions for existence 
of entanglement much easier as compared with the case that the system only includes 
nuclear spins. 

Quantum computation by means of ENDOR has been studied step by step. Sample 
study is shown to be very critical step as to detect sample with long decoherence time 
and well characterized magnetic tensors. As a result, molecular-spin bus systems for 
ENDOR QC are found. This molecule has been shown to be appropriate for elaborated 
experiments that may come later. 

In the other step for a full quantum computation by means of ENDOR, quantum 
operations for ENDOR QC are demonstrated. Different magnetic field, X-band, Q-band, 
and W-band from the range of 9.5 GHz to 95 GHz have been used for demonstration 
of quantum operations at different temperatures from up to ~ 300K to down to ~ 3K. 
The experimental results are shown to be very much promising on reliability of ENDOR 
physical system for quantum computation and quantum information processing. 

The other very important step, has been measurement for quantum computation. 
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Because of some technological restriction in ENDOR spectroscopy, measurement, at a 
glance may seem to be problematic point for ENDOR quantum computing. However, 
we have modified ENDOR technology in a way that is very well organized for quantum 
computing. For instance, quantum state characterization by TPPI is established and 
has been used widely for investigation of entanglement. This technique is applicable to 
solution ENDOR QC experiments if suitable molecular open-shell entities are prepared. 
W-band ENDOR experiments with the required conditions for entanglement will give 
evidence of establishment of entanglement in a more straightforward manner, disclosing 
nature of entanglement and its new aspects possibly applicable to advanced electron 
magnetic resonance spectroscopy. 

Finally, it should be pointed out that in the context of molecular-spin based EN- 
DOR QC fine structure terms appearing in the spin Hamiltonian are useful for acquiring 
an effective temperature cooling down effect because the zero-field splitting due to the 
fine structure interaction enormously enhances the energy gap between two electronic 
sublevels involving EPR transitions. 
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